| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,542,683 | 0 | T | . | 61.9% | 3.1 / 36.9 | 21 | intergenic (+50/‑432) | fimE/fimA | regulator for fimA/type 1 fimbriae major subunit |
| Reads supporting (aligned to +/- strand): ref base T (3/5); new base . (6/7); total (9/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTA > NC_000913/4542474‑4542930 | ggCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTACCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTAtttttt > 2:366223/1‑251 (MQ=255) gaggggggtgCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGGGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa < 2:55271/245‑1 (MQ=255) aaGATACTCGTTTAATTAAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaaa < 2:179658/250‑1 (MQ=255) aCTCGTTTAATTCAGGATTATCTCGGGAATCAAAATATCCGCCATACTGTGGGTTATTCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCGCATAAACGAAAAATTAAAAAGAGAAGAGGTCTGATTTAACTTAGTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt < 1:230876/250‑1 (MQ=255) ttttAATTCAGAATTACATCGGTCATCAAAATATTCGCCAAATTGGGCGTTATACCGCCAGTAATGCTGCTCGTTTTGACGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAG‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaa < 2:85719/249‑1 (MQ=255) gTTTAATTCAGGATTATCTAGGCCATCGAAATATTCGCCATACTGTGCTTAATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaaa < 2:70536/251‑1 (MQ=255) aTTATCTCGGGCCTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGA‑AATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGc > 1:271313/1‑250 (MQ=255) cGGGCATCGAAATATTCGCCATCCAGAGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATACTATCATAAAAGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAACACAAATAAATACAAGACAAT‑GGGACAAAACTGTCAATATCATACATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa > 2:360335/1‑250 (MQ=255) tatTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGACGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAA‑TGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgt > 1:156293‑M1/1‑162 (MQ=255) tgCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAtt > 1:369602/1‑251 (MQ=255) ccGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAA‑TGGTGCCAttttgactcatagaggaaagcatctcggacaatctttttcagtttatttgttggcttaatattcgatggttatctttatttatagatgtttatattgcatgaggtggttttt < 1:19913‑M1/251‑113 (MQ=255) aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAA‑TGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagagaagaatgaggaagatgcgtcgagccacagaaacgttagctttacata < 2:5628‑M1/251‑166 (MQ=255) ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAA‑TGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagagaagaatgaggaagatgcgtcgagccacagaaacgttagctttacatatagcggaggtgatgtgaaattaa < 1:152448‑M1/251‑189 (MQ=255) ttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAA‑TGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagagaagaatgaggaagatgcgtcgagccacagaaacgttagctttacatatagcggaggtgatgtgaaattaatttacaa < 1:101517‑M1/250‑196 (MQ=255) ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAA‑TGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagagaagaatgaggaagatgcgtcgagccacagaaacgttagctttacatatagcggaggtgatgtgaaattaatttacaatagaa > 1:22700‑M1/1‑50 (MQ=255) tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCa > 1:5760/1‑250 (MQ=255) tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTACAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAACTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCa > 1:323514/1‑250 (MQ=255) aataaAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATaaa < 1:270591/250‑1 (MQ=255) ataaAGTTAAAAAAACAAATAAATACAAGACAAT‑GGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACa < 1:345257/251‑1 (MQ=255) aaaaaaCAAATAAATACAAGACAA‑TGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagagaagaatgaggaagatgcgtcgagccacagaaacgttagctttacatatagcggaggtgatgtgaaattaatttacaatagaaataatttacatatcaaaca < 2:38423‑M1/251‑220 (MQ=255) cAAATAAATACAAGACAA‑TGGGGCCAttttgactcatagaggaaagcatcgcggacaacctttttcagtttatttgttggcttaatattctattgttacctttctttatagatgtttatattgcatgaggtggtttttggagagaagaatgaggaagatgcgtcgagccacagaaacgttagctttacatatagcggaggtgatgtgaaattaatttacaatagaaataatttacatatcaaacagttagc > 1:346151‑M1/1‑26 (MQ=255) caatGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTa < 2:149721/248‑1 (MQ=255) | GGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTA > NC_000913/4542474‑4542930 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |