Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,001,911 | 0 | G | A | 43.8% | 52.0 / 27.5 | 48 | intergenic (‑122/+199) | fliA/fliC | RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin) |
Reads supporting (aligned to +/- strand): ref base G (26/1); new base A (0/21); total (26/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.86e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC > NC_000913/2001770‑2002059 | gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc > 2:169736/1‑151 (MQ=255) agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg > 2:185269/1‑151 (MQ=255) agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg > 1:699147/1‑151 (MQ=255) aaTTCACGATAAACAGCCCTGCGTTATATGATTTATCGGCATGATTATCCTTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGTTTATTG‑GGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCgg < 1:73108/151‑1 (MQ=255) cGATAAACAGCCCTGCTTTATATGAGTTACCGGCATGATTATCCTTTTCTGCAGGTTTTTAAACCGGACGATTAGTGGGTAAAATGGGGGGTTATTG‑GGGGGTTACAGGTAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAg < 1:189267/150‑1 (MQ=255) aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc > 1:825444/1‑151 (MQ=255) aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc > 1:466848/1‑151 (MQ=255) aaaCAGCCCTGCGTTATAGGATTTAGCGGAATGTTTTTCGGTTTCTGGGGGGTTTTAAATCGGACAATTGGGGGGTAAAAGGGGGGGTTTTTT‑GGGGGTTACAGTTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc < 2:152219/149‑1 (MQ=255) aaCAGCCCTGCTTTATATGAGTTATCGGCATGTTTATCCTTTTCTGCGGGGTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGGTTTTTG‑GGGGTTTACGGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt < 1:540590/151‑1 (MQ=255) aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCt > 2:499673/1‑151 (MQ=255) cccTGGGTTTTATGATTTTGCGGCTTGATTACCCGTTTCTGGGGGGTTTTTAATCTGACGATTGGGGGGAAAAATGAGGGTTTTTTG‑GGGGTTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa < 1:1088925/151‑1 (MQ=255) cccTGCTTTATATGATTTATCGGCATGTTTATCCTTTTCTGCAGGTTTTTAAATCGGACGATTAGGGGGTAAAATGAGGGTTTTTTT‑GGGGTTTACAGTTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTg < 2:13325/150‑1 (MQ=255) ttATATGAGTTATCGGCATATTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTTTTG‑GGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa < 1:393355/148‑1 (MQ=255) ttacgggcttaaTTTTCCGTTTCTGCGGTTTTTTAAATCGAACGATGGGGGGGTGAAAGGGGGGTTTTTTG‑GGGGGTAAAAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt < 2:1037275/137‑1 (MQ=255) tagcGGCATGATTATCCGTTTCTGCAGGTTTTTAAATCGGACGATTAGGGGGTAAAATGAGGGGTTTTTT‑GGGGGTGACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg < 2:896836/147‑1 (MQ=255) aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg > 1:1023072/1‑150 (MQ=255) aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGAAAATTCCAGGCAGAAAAAAACACCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg > 1:153407/1‑150 (MQ=255) gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg > 2:78352/1‑151 (MQ=255) ttATCCTTTTCTGCGGGTTTTTTAATCGGACGATTGGGGGGTAAAAGGGGGGTTTATTT‑GGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc < 2:389777/150‑1 (MQ=255) ttATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAAGGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc < 2:258275/150‑1 (MQ=255) aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc > 1:705440/1‑151 (MQ=255) aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc > 1:550314/1‑151 (MQ=255) ccttttCTGCGGGTTTTTTAATGGGACGTTTGGGGGGTAAAAGGAGGGTTTTTTG‑GGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg < 1:609406/147‑1 (MQ=255) ttccggaggggTTTGTAAGCGGTCAATTAGTCGTAGAACAATGGGTTTATTG‑GGGGCTTCCAGGGAAGTCCCATCCAAAAAAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg < 1:244070/137‑1 (MQ=255) ttCTGCGGGGTTTTTAATCGGACAATTAGTGGGTGAAAGGGGGGGTTTTTG‑GGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa < 1:977828/150‑1 (MQ=255) tGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc > 1:518978/1‑151 (MQ=255) gggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt > 2:392704/1‑151 (MQ=255) tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc > 1:355208/1‑151 (MQ=255) aaTCGGACGATTAGGGGGTAAAATGGGGGTTTTTTT‑GGGGGTTAAAGGTAAATCCCAGG‑AAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc < 2:260220/151‑1 (MQ=37) acgggaggtttggTGGGGAAAATGGGGGTTTTTTT‑GGGGGTTCCAGAAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc < 2:836687/139‑1 (MQ=25) ggACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa > 1:1052019/1‑150 (MQ=255) ggACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa > 1:321577/1‑150 (MQ=255) ggACGATTAGGGGGTGAAAGGGGGGTTTTTTT‑GGGGGTTAAAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa < 2:368889/150‑1 (MQ=37) aTTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg > 1:194312/1‑150 (MQ=255) ttAGGGGGTGAAAGGAGGGGTTATTG‑GGGGGTTACAGTTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt < 1:39751/150‑1 (MQ=38) gggTAAAAGGGGGGTTTTTTT‑GGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa < 2:343042/148‑1 (MQ=37) gTGAAATGAGGGGTTATTT‑GGGGATTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc > 2:34419/1‑151 (MQ=255) aaaTGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt > 2:859223/1‑151 (MQ=255) aTGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa > 2:355415/1‑150 (MQ=255) gggggTTTTTTGGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa < 1:884312/148‑1 (MQ=37) ttt‑GGGGGTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt < 1:327745/149‑1 (MQ=255) gggtTTACAGTAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt < 1:78345/151‑1 (MQ=37) cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc > 1:193294/1‑151 (MQ=255) ggTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg > 1:921776/1‑151 (MQ=255) aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa > 2:301037/1‑150 (MQ=255) aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa > 1:864936/1‑150 (MQ=255) aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc > 2:193067/1‑149 (MQ=255) aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCATTTCAACTTGTAGGCCTGATAagcgcagcgc > 1:69087/1‑150 (MQ=255) | GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC > NC_000913/2001770‑2002059 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |