| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,634,578 | 0 | T | C | 57.1% | ‑3.3 / 16.6 | 21 | E103G (GAA→GGA) | tfaQ | putative tail fiber assembly protein TfaQ |
| Reads supporting (aligned to +/- strand): ref base T (6/3); new base C (4/8); total (10/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.98e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGC > NC_000913/1634410‑1634798 | tcctccGTTGCAATTTCCCGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGGAAAGTACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAAc > 1:198074/1‑251 (MQ=9)tccgccGTTGCAAGTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCTCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATaaa < 1:5903/246‑1 (MQ=14) aTTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCg < 1:258884/251‑1 (MQ=11) ttCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGCAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCgat > 2:157286/1‑251 (MQ=14) cTGCAGCATCCTGAAGCTGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCATGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTc > 1:102672/1‑250 (MQ=14) cTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACAGAGTTCAGAATTAAATAACGCGTCGACGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTc > 2:197705/1‑250 (MQ=14) gcatccgtaaGCGGCGCAATATGCTCACTGGCTACCTGCATCATGCTTTTTTTTGTTTCGTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCGTATATTCCCCTCCCGGCGATAACCATGTAAAATGTTCCGGTAACGGACCGAGTGCATAAATAAATAACGCGTCGCCGGAATCCACGTCATAGTCGGTTTTACCCCaatggtcggcaacgtg < 2:151168/243‑13 (MQ=11) ccTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATg < 1:224576/250‑1 (MQ=11) gcAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtc < 2:3987/251‑1 (MQ=31) aCTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGCATC‑AGAACAAGTTTTTCTGATTCCGTCTCATTCACCCAGGATGTGACGTTACACTTCAGAAACTCCACTTACGGCGACAACCAGGCAACATTTTACGGTAACGAACCGAGTTCAGAACTAAACAACGCGTCAAAGGAAGCAACGGCCTAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTg > 1:17951/1‑250 (MQ=2) ggCTAGCTGCATCCGGCTTTTTTATGTTTATTACGCCTCCCGGATCCGGCACA‑GTTTTTCTGCTATCGTATCCTTAACCCAGGCAGTGCCGTTCCCATTCTGATATTCCCCTCCCGGCGATAACCCGGTAAAATTTTCCGGTAACGGACAGGGATCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTATTCAACGAGATGCCACGATGCCTCATAACTGTTGaaa > 2:150297/1‑250 (MQ=255) aTCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCAcaa > 1:216693/1‑248 (MQ=11) ggattttttttGTTTCTTCCTCCTACCGGATCCGGTACA‑GTTTTGCTGCTCCCGTATCCTTCCCCCAGGCTGGGCCGTTCCACTTCTGAATCTCCCATTCCGACGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTTCCCCGACGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagcca < 2:714/246‑1 (MQ=11) tttttttGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAagccagc < 2:192312/248‑1 (MQ=11) ttGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAat < 1:70168/250‑1 (MQ=11) tcCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCg < 1:217080/250‑1 (MQ=11) cA‑GTTTATCTGCTTCCGCATCCTTCACGCAGGCTGAGCCGTTCCACTTCTGTTATTCCCCTCCCGGCGATACCCCGGTAACATTTTCCGGTAACGGACCGACTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGCTGCCACGATGCCTCATCACTGTTGAAAACAACCCCAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTg > 2:104920/1‑251 (MQ=11) gCTTCCGTATCCTTCACCCAGGCTGGTCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAAAAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGAGGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATGTCGGTACTGTTTgcaggcag < 2:255492/248‑1 (MQ=11) cGTATCCTTCACCCAGGCTGTGCCGTTCCCCTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCCCGTCATAGACGGTTTTACCCCGATGGTCTTCACCGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACCAAGCCCGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTAt > 1:74626/1‑251 (MQ=11) ccTTCCCCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCGGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAgg < 2:20559/250‑1 (MQ=11) gCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGc > 1:10989/1‑251 (MQ=31) | TCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACA‑GTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGC > NC_000913/1634410‑1634798 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |