Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,684,749 | 0 | A | C | 57.7% | 26.7 / 39.2 | 26 | D164A (GAC→GCC) | tus | inhibitor of replication at Ter, DNA‑binding protein |
Reads supporting (aligned to +/- strand): ref base A (11/0); new base C (0/15); total (11/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-06 |
AAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTG > NC_000913/1684613‑1684881 | aaaCTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTAc > 2:341211/1‑151 (MQ=255) aGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttgg > 2:66241/1‑151 (MQ=255) cGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaat > 1:659066/1‑151 (MQ=255) tatCGCCACGGTTAAATCAGACCCCCCCCCCGCGGCACGTTTTGAAGGGGTGCATCGTCTTTGGCGGGGGCTGATCCCCCTTAATGCTTCCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCtaata < 1:400413/146‑1 (MQ=255) cgcgtttaaaTCAGACCCCCCCCCCGGGGCAGGTTTTGAAGGGGTCCTTCTTCATTTGCCGGGGCGGACCGCCCTTAATGCTTCGGGCGCGCTAACCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACat < 2:594536/137‑1 (MQ=255) cACGGTTGAATCAAACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTTACCGCCCGCCCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATc < 1:614244/148‑1 (MQ=255) tttaaTAAGAACTCCCCCCCGCGCCCCGTTTTGAAGGGGTGCATCGTAATTTGCCGGGGCTGATCCCCCTTAATTTTTACCGCCCGCCCACCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGa < 1:502147/147‑1 (MQ=255) taaaccaaaaCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCCTTAATGCTCCCCCCCCCCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGaa < 2:781753/142‑1 (MQ=255) aaTCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAttt > 1:237274/1‑151 (MQ=255) agacccccccccCGCGTCCCGTTTTGAATGGGTGCACCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTCCCGCCCGCTCACCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACAt < 1:263884/145‑1 (MQ=255) ccgccccccCGCGCCCTGTTTTGATTGGGTGCATCTTCTTTTGCCGGGGCTGACCGCCCTTATTGTTTACCGCGCGCCCACCGTTCTTCACACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGt < 2:177314/148‑1 (MQ=255) ccccccccCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCTTTAATGCTTACCGCACGCTCACCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGt < 1:169868/148‑1 (MQ=255) ccgtttttAAAGGGGTGCATGGTCTTTTGCCGGGGCTAACCCCCCTTATTGCTTCCCGCGCGCTCCCCTTCCGGCACACCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGc < 1:750618/147‑1 (MQ=255) gTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGcac > 1:306098/1‑150 (MQ=255) ttttGAATGGGCGAATCGTCTTTGGCCGGGGCTGATCCCCTTTAATGTTTACCGCCCGCCCACCGTTCTCCACGCCACCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGca < 2:237329/148‑1 (MQ=255) aaTGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTgg > 2:179120/1‑150 (MQ=255) ttGCCGGGGCTGACCCCCCTTAATGCTTCCCGCCCGCTCACCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATca < 1:788637/150‑1 (MQ=255) ggggCTGCCCCCCTTTAATGCTTCCCCCCCCCTCCCCGTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGc < 2:328387/151‑1 (MQ=255) ggggCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGc > 1:666928/1‑151 (MQ=255) gggCTGATCCCCTTTAATGCTTACCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCa < 1:792057/151‑1 (MQ=255) tGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgtgt > 1:217541/1‑151 (MQ=255) cTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcg > 1:517737/1‑151 (MQ=255) aCCGCACGCTCACCGTTCTGCACGACCCCGCCCCTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg > 1:775099/1‑150 (MQ=255) ccGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg > 1:244971/1‑149 (MQ=255) ccccGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCgagg < 2:196224/141‑1 (MQ=255) cacgtttttccACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTg < 2:59360/137‑1 (MQ=255) | AAACTCAAGACCACGTTCGAGCATATCGTCACGGTTGAATCAGAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATGCTTACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTG > NC_000913/1684613‑1684881 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |