| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,429,207 | 0 | T | G | 52.4% | 10.8 / 14.7 | 21 | V53V (GTT→GTG) | stfR | putative prophage side tail fiber protein StfR |
| Reads supporting (aligned to +/- strand): ref base T (4/6); new base G (8/3); total (12/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.98e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCC > NC_000913/1428958‑1429441 | aTCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACATCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGTTGAATATGGCAGTAAAAATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAACCAGCTGAAAGCAAAACGTAACAGCACCACTGTGGTGGTGAACACCCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTGACAGCATGGACGTTg < 1:212929/251‑1 (MQ=255) gCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAgtgtcatcctgcagg > 1:85562/1‑237 (MQ=255) tgGGGCGAATATGGCAGTAAAGATTTCAGGAGTACTGAAAGACGGCACAGGAAAAACGGTACAGAACTGCACCATTCAGATGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATACTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGAtt > 1:65686/1‑251 (MQ=255) ggTGANNATNNNAGNNAAGATTTCAGGTGTACTGAAAGANGGCACAGGAAAACCGGTACAGAACTGCACNANCNAGCTNNAAGNAAAACGTAACNGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGCTGAAGNCNGGCGTTACAGCATGGACGTTGAGT‑ACGGTCANTACAGCGTTATTCTGTTGGTGGAANGATTNncgncnnnnnnnnnn > 2:159747/1‑213 (MQ=255) aggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCAACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg > 1:15394/2‑251 (MQ=255) gCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGAGGTGAACACGGAGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCa > 1:193195/1‑243 (MQ=255) aaaaCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGGAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGt < 1:188498/122‑1 (MQ=255) aaCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGGGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATGCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGcc < 1:159747/251‑1 (MQ=255) aGAACTGCACAAACCGGCGGAAAGCAAAAGGTAGCAGCACCACGGGGGTGGTGAACGCGGTGCCCTCAGAAAATCCGGATGAATCCGTCTGTTACAGCATGGACGTTG‑GTGAGGGTAATTATGGCGTCACTGGGTTTGTGGAAGGATTCCCGCCGTCACATGCCGGGGCTATTGCCGTGTATGAAGATTCTCAGCCCGGGACGCTGAATGATTTTCTCGGTCCCAGGACGGAGGATGATGCCCGTCCGGAg < 1:198554/251‑1 (MQ=255) caATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTgcgc > 1:169723/1‑251 (MQ=255) gAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTTTCATCCTGCAGGTGGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGc < 1:178861/222‑1 (MQ=255) cacGTTGGCCTAAGAAAATCCGGATGAAGCCGGGCGTTACAGCAGGGAGGTGGAGT‑ACGGTCAGTACAGCGTTATTCTGTGGGTGGAAGGCTTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACGGATGGTGGAAGAGGTGGTGCGTAACGCGTCCGCgg < 1:172378/251‑1 (MQ=255) gcTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTAAAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGCTTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCAATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTATTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGt > 1:150114/2‑237 (MQ=255) cctaaGAAAATCCGGATGAAGCCGGGCATTCCAGCATGGACGTTGAGT‑ACGGTCGGTACAGCGGTATTCTGTTGGTGGATGGTTTCCCGCCGTCACATGCCGTGACCACTACCGTGTATGAAGAGTCCCAACCCGGTACGCTGAATGATTGTCTCGGTGACATGACGGAGGATGATGCCCGTCCGGAGGCCCTGCGCCGTTTTGAACTGATGGTGGACGAGGTGGGGCTTAACGCGTGCGCGGTGGCACAg < 1:16949/247‑1 (MQ=255) gaaTCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGGGGTGGCACAGagtacgg > 1:93871/2‑245 (MQ=255) ggATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGCCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCCCCATCGCACGCCGGGACCATCACCGTGTATGCAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGCAGAGGTGGAGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgacg > 1:29840/1‑248 (MQ=255) gggCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTANNGTGTAtnn > 2:17771/1‑96 (MQ=255) gggCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCAGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTg < 1:17771/222‑1 (MQ=255) gTTACAGCATGGCTGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGCTGACGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGcc > 1:218268/1‑251 (MQ=255) gTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGATGTGCTGCGTCGTCGTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGcc < 1:73144/251‑1 (MQ=255) gTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGcc > 1:99726/1‑251 (MQ=255) | ATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCC > NC_000913/1428958‑1429441 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |