| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,634,724 | 0 | G | A | 57.7% | ‑5.6 / 25.3 | 26 | F54F (TTC→TTT) | tfaQ | putative tail fiber assembly protein TfaQ |
| Reads supporting (aligned to +/- strand): ref base G (7/4); new base A (7/8); total (14/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.53e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATTA > NC_000913/1634478‑1634887 | cagtcttttttttGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGGGCCGTTCCACTTCTGATATTCCCCTCCCGGCGAGAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAag < 2:65943/243‑1 (MQ=11)caggcttttttttGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCACCACTGTTGAAAACA‑GCCACAAagc > 1:122160/8‑251 (MQ=11)cAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGCTGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAA‑CACGCCACGAagc > 2:276596/1‑251 (MQ=9) ttttttGTTACTTCCGCCTCCCCGACCCGGAACAGTTTTTCTGCTTCCGTCTCCGTCACCCAGGCTTTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCGCGTCATAGACGGTTTTACCCCGATGGGCTTCAACGAGATGCCACGATGCCTCATCACTGTTTAAATCA‑GCCCCAAAGCCAGCCg < 2:112562/250‑1 (MQ=11) tttttGCTTCTTCCGACTCCCGGATCCGGAACAGTTTTTCTGCATCCGTAACCTTCACCCAGGCTGTGCCGTTCCACTTCTGCTATTCCCCTCCCGGCGATAACCAGGCAAAATTTTCCGGCAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCAACGTCATAGACGGTCTTACCCCGATGGTCTTCAACGAGATGCCACCATGCCTCATCACTGTTGAAAACA‑CCCACACAGCCAGCNggc > 2:262200/1‑250 (MQ=11) cTCCCGGATCCGGAACAGTTTTTCTGCTCCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGACACGATGCCTCATCACTGTTGAACACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCa > 2:132663/1‑250 (MQ=11) cccGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAat > 2:53523/1‑250 (MQ=31) ggATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACACAGGCTGTGCCGTTACACTTATGAAACTACCCTTCCGGCGATAACCAGGTAACCTTTTCCGGTCACGGACCGAGATCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACAGTTTTAGCCCGATGATCTTCAACGAGATGACACAATGCCTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCAGGCGTTGCAATCTCg > 1:154346/1‑251 (MQ=11) gAACAGTTTTTCTGCTTCCGTATACTTCACCCAGGCTGTGCAGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGATCAGAAATAAATAACGCGACGCAGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GACACGAAGCCAGCCGGAATATCAGGCGGTGCAATATCGGTACTGt > 2:275872/1‑251 (MQ=11) gTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCCCTGTTTAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGc < 1:167460/250‑1 (MQ=11) ttttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATAGGTACTTTTTGcagg > 1:35696/1‑248 (MQ=11) tCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTg < 2:49556/250‑1 (MQ=31) tCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCgg < 2:201881/250‑1 (MQ=31) cACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCTCGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGaa < 1:152209/251‑1 (MQ=11) gCTGTGCCGTTCCACTTCCGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCGGTTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCAGAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGAGGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGTCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGc < 1:269444/251‑1 (MQ=11) gtgCCGTTCCACTTCTGAAACTCCCCTTCCGTCGATAACCAGGTAACATTTTCCGTTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGTAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCa < 2:196572/249‑1 (MQ=14) tCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTc > 1:110437/1‑250 (MQ=11) ttCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGAACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAAt > 2:181379/1‑251 (MQ=14) tGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATaa > 2:77342/1‑250 (MQ=31) attCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTc > 1:177559/3‑250 (MQ=11) ggTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATtat < 2:195014/251‑1 (MQ=11) aTTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAAttt > 1:11385/1‑251 (MQ=11) tttCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAAttttt > 2:103370/1‑251 (MQ=31) cGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTCCCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATgg < 2:6980/251‑1 (MQ=11) tCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACA‑GCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCAc < 2:35696/250‑1 (MQ=11) gTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATTa < 2:74590/251‑1 (MQ=31) | CAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACA‑GCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATTA > NC_000913/1634478‑1634887 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |