| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,001,911 | 0 | G | A | 41.0% | 48.4 / 17.9 | 39 | 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 (22/1); new base A (0/16); total (22/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.51e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCG > NC_000913/2001766‑2002053 | ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGaaaa > 1:361229/1‑150 (MQ=255) cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGaaaaaa > 2:762297/1‑150 (MQ=255) agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTg > 2:265581/1‑151 (MQ=255) tccccaaaaacccgCCGTGCTTTATAGGGGTTATCGGCATGATTATCCTTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGGGGGGTTATTTGGGGG‑TTACAGGTAATTTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCggg < 2:409940/137‑1 (MQ=255) aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACg > 1:530385/1‑150 (MQ=255) aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGt > 2:266783/1‑151 (MQ=255) aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCt > 2:179132/1‑151 (MQ=255) gCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg > 1:119547/1‑151 (MQ=255) ggTTATATGAGTGATCGGCATGGTTAAACGTTTCTGCAGGGTTTTT‑AATCGGGCGGTTAGTGGGTGAAATGAGGTGATATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCAAGTTGCTGACAAAt > 1:344893/2‑151 (MQ=255) tataTGAGTTTTCGGCATATTTATCCTTTTCTGCAGGTTTTTT‑AATCGGACGATTAGGGGGAAAAATGGGGGTTTTTTTGGGGG‑TTACAGAAAATTTCCAGCCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc < 2:144741/144‑1 (MQ=255) atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg > 2:505214/1‑151 (MQ=255) tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc > 1:729021/1‑151 (MQ=255) aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt > 1:404663/1‑150 (MQ=255) ggCAGGATTTCCCGTTTCTGCGGGTTTTTT‑AATCGGACGATTAGGGGGTGAAATGAGGGGTTATTGGGGGG‑TTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt < 1:31282/148‑1 (MQ=255) aTGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg > 2:451516/1‑150 (MQ=255) gATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg > 2:339562/1‑151 (MQ=255) ttatcGTTTTTGGGGGGTTTTTAAATGGGAGGTTTAG‑GGGGGAAAAGGGGGGTTTTTTGGGGG‑TTTCGGGGAAATTTCCGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg < 1:584702/145‑1 (MQ=255) tATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg > 1:220035/1‑151 (MQ=255) ttttgaaGGTTTTTT‑AATGGGACGTTTAGGGGGTAAAAGGGGGGTTTATTGGGGGT‑TTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaaa < 2:843201/144‑1 (MQ=255) tCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg > 2:53252/1‑151 (MQ=255) aTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc > 2:655446/1‑151 (MQ=255) ggACGATTGGGGGGTAAAATGGGGGGTTATTTGGGGG‑TTACAGGAAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa < 1:404346/150‑1 (MQ=37) ggACAATTAGGGGGTAAAAGGAGGGTTTTTTGGGGGG‑TTACAGTTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa < 1:34094/151‑1 (MQ=37) gTGGGTAAAAGGAGGGTTTTTTTGGGGG‑TTAAGGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg < 2:796371/148‑1 (MQ=37) gggTAAAATGGGGGTTTATTGGGGGT‑TTACAGGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc < 2:413281/147‑1 (MQ=37) tGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa > 2:168483/1‑150 (MQ=255) tGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa > 2:178962/1‑150 (MQ=255) atgggggTTTTTTGGGGGTTTAAAGTTAAATTCAAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa < 2:698545/148‑1 (MQ=25) gtttttttGGGGGTTTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt < 2:273325/144‑1 (MQ=37) ggTTATTGGGGGT‑TTAAAGATAAATCCCAGGAAAAAAAACA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa < 2:307037/151‑1 (MQ=25) tATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt > 2:597799/1‑151 (MQ=255) ttgggggt‑TTACAGGTAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt < 2:273761/146‑1 (MQ=37) ttGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa > 1:685368/1‑150 (MQ=255) tggggtg‑TTACAGATAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa < 1:810688/148‑1 (MQ=37) ggggg‑TTACAGGAAAATCCCGGGAAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa < 2:876550/149‑1 (MQ=25) gggt‑TTACAGATAAAGCCCAGGAAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt < 2:559467/151‑1 (MQ=255) ggTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg > 1:69763/1‑151 (MQ=255) ttCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag < 2:799489/150‑1 (MQ=255) tCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg > 2:328695/1‑151 (MQ=255) | TTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCG > NC_000913/2001766‑2002053 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |