| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,429,160 | 0 | C | G | 61.1% | 12.2 / 17.4 | 19 | L38V (CTG→GTG) | stfR | putative prophage side tail fiber protein StfR |
| Reads supporting (aligned to +/- strand): ref base C (3/4); new base G (7/4); total (10/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.30e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGG > NC_000913/1428918‑1429406 | tGAAGGCCCCGGCAGTGGCGACGGGCGCACTGACGGTTTCATCGTGGGTTTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTAGGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTCCAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCtn < 1:133480/250‑2 (MQ=255) gCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACCGTGTTATGACAGCCCGCCGGACCCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTACAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGc > 1:197290/1‑251 (MQ=255) cACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGACGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTTAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTa > 1:119838/1‑250 (MQ=255) aTCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACATCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGTTGAATATGGCAGTAAAAATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAACCAGCTGAAAGCAAAACGTAACAGCACCACTGTGGTGGTGAACACCCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTGACAGCATGGACGTTg < 1:212929/251‑1 (MQ=255) gCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAgtgtcatcctgcagg > 1:85562/1‑237 (MQ=255) ccggaaccggtgggCTTTTTTGTGGGGTGAATATGGCAGGAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACTGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGc < 1:148227/171‑1 (MQ=255) gnaaccggtgggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCCCGGTGGTGGGGAACACGGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGatgtgg > 1:170339/10‑185 (MQ=255) tgGGGCGAATATGGCAGTAAAGATTTCAGGAGTACTGAAAGACGGCACAGGAAAAACGGTACAGAACTGCACCATTCAGATGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATACTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGAtt > 1:65686/1‑251 (MQ=255) ggTGANNATNNNAGNNAAGATTTCAGGTGTACTGAAAGANGGCACAGGAAAACCGGTACAGAACTGCACNANCNAGCTNNAAGNAAAACGTAACNGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGCTGAAGNCNGGCGTTACAGCATGGACGTTGAGT‑ACGGTCANTACAGCGTTATTCTGTTGGTGGAANGATTNncgncnnnnnnnnnn > 2:159747/1‑213 (MQ=255) gCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTAc > 1:109053/1‑144 (MQ=255) gtccTGAAAGACGGCACAGGAAAACCGGTACAGAACGGCCCCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTgggct < 1:121742/94‑4 (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) | TGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGG > NC_000913/1428918‑1429406 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |