| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,634,661 | 0 | C | G | 60.0% | ‑5.4 / 16.0 | 15 | T75T (ACG→ACC) | tfaQ | putative tail fiber assembly protein TfaQ |
| Reads supporting (aligned to +/- strand): ref base C (3/3); new base G (3/6); total (6/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.22e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTT > NC_000913/1634427‑1634866 | cAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCgatggt > 1:50650/1‑249 (MQ=11) cgcAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGACTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACt > 1:72693/1‑251 (MQ=14) gCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGCTCCGGAACAGTTTTTCTGCTTCCGTATACTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGCGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCCCGCTGACTCATCACTg > 1:148278/1‑251 (MQ=9) ccTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGcc < 1:63897/251‑1 (MQ=11) tnttcNGNNTNNCGANTCCGGAACAGTTTTTCGGCTTNNGTGTCCTTCNCCCAGTCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGGANAATTTTTCGNNAANGGACCNANNNCNNNAANAAATAACGCGTCGCCGGAANCCNNGTCATANACGGTTTTACCCCGATGGTCTTCAAc < 2:35106/191‑1 (MQ=2) ccGCCTCCCGGTGCCGGAACATTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGTTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCGCGATGACTCATCTCTGTTGAAAACAGCTACGAAGCCAGCCGGAATATCTGGCGGt < 1:26805/251‑1 (MQ=11) cAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTGCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTGTACCCCGATGATCTTCAACGAGATGCCACGATGACGCAGCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATAGCGGTACTGTTTg < 1:24563/251‑1 (MQ=9) gTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGAGCCGTTACACTTCTGTCATTCCCCACTCGGAGCTAAACAGGAAAAATTTTCCGGTAACCGCCCGCGTCCAGTAATAAATAACCCGCCGCCGGAAGCCAAGACATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCCCGAAGCCTCATCACTGTTGAAAACAGCCAACAAGCCCGCCGGAATATATGGCGGTGCACTAACGGAACTGTTTGc > 1:78663/1‑250 (MQ=2) tCTGCTCCCGTATCCTTCTCCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACTGACCGAGTTCAGAAATAAATAACGCGTCTCCGGAAGCCACGTCATAGTCGGTTTTCCCCCGATGTTCTGCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAg < 1:113354/196‑1 (MQ=14) ggCTGTGCCGTTCCACTGCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAATCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACTATATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTTCTGTTTGCTGTCAGACCTGTATGAGGCGGAATATATg < 1:84781/251‑1 (MQ=255) gtgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATg > 1:57332/1‑247 (MQ=31) ggACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTgg < 1:130392/251‑1 (MQ=11) gAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGgttgtt > 1:108179/1‑251 (MQ=11) naGANATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCgt < 2:105493/137‑2 (MQ=11) anntaaCGCNNCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCAcaaa < 2:10540/99‑4 (MQ=11) | CAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTT > NC_000913/1634427‑1634866 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |