| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,001,911 | 0 | G | A | 52.2% | 20.9 / 56.4 | 46 | 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 (17/5); new base A (0/24); total (17/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.51e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-07 | |||||||||||
CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGC > NC_000913/2001768‑2002054 | cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGaaaaaa < 2:58069/150‑1 (MQ=255)cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGaaaaaa > 1:468859/1‑150 (MQ=255) aaTTCACGAAAAACAGCCCTGGGTTATATGATTTATCGGCATGATTATCCGTTTCTGAAGGGTTTTTAATCGAACAATTAGGGGGTAAAAGGAGGGGTTATTGGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCgg < 2:76831/151‑1 (MQ=255) cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg > 2:378047/1‑149 (MQ=255) aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc > 2:326310/1‑151 (MQ=255) tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCGGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg > 2:226621/1‑151 (MQ=255) aaaaCAGCCCTGGGTTATATGAGTTATCGGCATGATTATCGGTTTCTGCAGGTTTTTTAAGGGGACGATTAGTGGGTGAAATGGGGGGTTATTGGGGGT‑TTACGGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACg < 1:274367/150‑1 (MQ=255) aaaCGGCCTTGCGTTATTTGTTTTATGGGCAGCGTTATCCTTTTCTGCAGTGTTGTTATTAGGATAATTGGTGGGTGAATTG‑GGGGTTTTTTGGGGGCTTACAGGTAAATTCCGTGCAAAAAGAACCCCCGCCGGTGGCGGGGAAGCACg < 2:506768/150‑1 (MQ=255) aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt > 1:543081/1‑151 (MQ=255) gCCCTGCTTTATATGATTTATCGGCATGATTATCCGTTTCGGCAGGGTTTTTAACGGGACGATTAGGGGGTGAAATGGGGGGTTATTGGGGGG‑TTACAGGTAATTTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTg < 2:278717/151‑1 (MQ=255) gCCCTGCTTTATATGAGTTACCGGCATAATTACCCGTTTCTGCAGGTTTTTAAATCGGACGTTTAGGGGTTAAAATGAGGGTTTATTTGGGGT‑TTACAGGTAATTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCt < 1:597514/150‑1 (MQ=255) cccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg > 1:181060/1‑150 (MQ=255) ggtttATATGATTTAGCGGCATGTTTTTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGGTAAATGAGGGTTTTTTTGGGGT‑TTAAAGTAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa < 1:248831/144‑1 (MQ=255) gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa > 1:126341/1‑150 (MQ=255) ttATCGGCATGATTATCCTTTTCTGCAGGTTTTTTAAGCGGACGATTAGGGGGTAAAAGGGGGGTTTATTTGGGGT‑TTACAGGTAATTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt < 2:578667/150‑1 (MQ=255) cGGCATGCTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg > 1:624091/1‑150 (MQ=255) ggtaTGATTATCCGTTTTTGCGGGTTTTTTAATCGGACTATTAGTGGTTGAAAGGAGGGTTTTTTGGGGGG‑TTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc < 1:22942/148‑1 (MQ=255) tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt > 2:262073/1‑151 (MQ=255) gATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg > 1:149157/1‑151 (MQ=255) ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc > 2:374933/1‑151 (MQ=255) tGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc > 1:148289/1‑151 (MQ=255) gcggggTTTTTAAGGGGACGATTGGGGGGTAAAAGGGGGGGTTTTTGGGGGG‑TTAAAGAAAATTCCCAGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt < 2:399354/148‑1 (MQ=255) ggttttttAATCGGACGATTAGTGGGTAAAATGAGGGTTTTTTTGGGGT‑TAACAGGTAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt < 1:143519/148‑1 (MQ=255) ttAATCGGAGGATTAGTGGGTGAAAGGGGGGTTTTTTTGGGGG‑TTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc < 2:567528/150‑1 (MQ=38) cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa > 2:366407/1‑150 (MQ=255) gaaagATAAGTGGGTAAAAGGAGGGGTTATTGGGGGG‑TTACAGTTAAATCCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa < 2:460571/147‑1 (MQ=37) aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa > 1:192283/1‑150 (MQ=255) aCGATTAGGGGGTAAAATGGGGGTTTATTGGGGGT‑TTACAGGAAAATCCCAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa < 1:536707/151‑1 (MQ=37) gTGGGTAAAAGGGGGGGTTTTTTGGGGT‑TTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa < 1:326317/150‑1 (MQ=37) aaaagggggggTTTTTTGGGGG‑TTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt < 2:192278/141‑1 (MQ=38) aaaGGAGGGTTTATTGGGGGT‑TTACAGTTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCaa < 1:619318/150‑1 (MQ=37) aTGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa > 2:514137/1‑150 (MQ=255) ggtttttttgggggg‑TTAAAGATAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt < 2:277516/139‑1 (MQ=25) ggggTTTTTTGGGGT‑TTACAGATAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt < 1:321706/150‑1 (MQ=255) gttttttgggggg‑TTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa < 1:183971/143‑1 (MQ=37) ggTTTTGTGGGGA‑TCACAGATAAATCCCGGCGAAAAAAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCTGATGCGGCGTAAACCCCTTATCCGCCCTACAAAAATGTGCCAATTCAATAAATTGCa < 1:119212/151‑1 (MQ=255) ggggg‑TTACAGTTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa < 2:183747/149‑1 (MQ=255) ggggaatataggctAATTTCAAGGGAAAACAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTGATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc < 2:443129/135‑1 (MQ=21) gggg‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt > 2:609093/1‑151 (MQ=255) aCAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg < 1:222067/149‑1 (MQ=255) aGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg < 2:551702/148‑1 (MQ=255) ggaaaaTGCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg < 1:526717/148‑1 (MQ=25) aaaTCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa < 1:287004/149‑1 (MQ=255) aaTTCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa < 2:534005/151‑1 (MQ=255) ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc > 2:129941/1‑151 (MQ=255) ccAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg < 1:249062/150‑1 (MQ=255) | CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGC > NC_000913/2001768‑2002054 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |