Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 783,298 | 0 | G | C | 33.9% | 100.4 / 11.1 | 56 | P141R (CCG→CGG) | greB_1 | Transcription elongation factor GreB |
Reads supporting (aligned to +/- strand): ref base G (5/32); new base C (19/0); total (24/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCTGCGCGTCGACAAACCGCACCCGGC > KT2440/783170‑783430 | ccGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCCGCCGCCGCCCG‑ACCACCCCAGCCCCCCGccccccc > 4:342573/1‑133 (MQ=255) ccGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTACCCCACCCGCCGCCGCCCC‑CCCACCCCAAACCCCCGccccccc > 8:285291/1‑133 (MQ=255) ccGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCGATACCCTGTTCCCCCCACCGCCGCCGCCAC‑CCCACCTCAGCCCCCCGccccccc > 8:128498/1‑133 (MQ=255) ggggCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCCCCTCTATCACCTGTTGCCCCACCGGCCGCCGCCAG‑AACACCAGCTCCCACCgccccgccccccgc > 6:333450/1‑133 (MQ=255) ggggCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGGTCCCCCACCGGCCGCCGCCAC‑ACCACCGCCTCCCCCCGCcccccccccccc > 4:247072/1‑133 (MQ=255) gggCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCa < 2:36465/137‑1 (MQ=255) gCCGCTTTGCGGCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCaaa < 3:29226/137‑1 (MQ=255) tttGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCCCCGGCCGCCGCCCG‑CCCCCCTCATCCCCCCGCCCCCCCCCGCGCCAAccccc > 5:402907/1‑136 (MQ=255) tttGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCCCCGGCCGCCGCCAG‑CCCACCTCATCCCCCCGCCCCGCCCCCCGCAAAgcccc > 5:233613/1‑136 (MQ=255) tGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTGCCCCCCCGGCCGCCGCCCG‑ACCACCGCATCCCCCCGCCCCCCCCCGCGCCAAgcccccc > 2:125899/1‑137 (MQ=255) ttGGGATGCAAAGAAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCg < 7:128498/136‑1 (MQ=255) tGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCCCCCGCCGCCGCCCG‑ACCCCCACAACCCCCCGCCCACCCCCCCGCCAACCCCCCCCCCAAGGCgc > 2:374593/1‑136 (MQ=255) cTGCACAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAAccc < 1:437415/137‑1 (MQ=255) cTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAAccc < 1:269010/137‑1 (MQ=255) aaaGCCGCCCCAGGGCCCTCAACCCACGCCCCGAATCTCAACAACCTCTATGACCTGTTAACCCCCCGGCCGACGACAG‑ACCACCACAACACCCCGCCCAGCCCCCCGCAACGACCCCCCCAAAGCCGAACCCCAgg > 2:338466/1‑136 (MQ=255) aaGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGtt < 2:8380/137‑1 (MQ=255) cgcagcCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTg < 6:121313/136‑1 (MQ=255) ccccAGGGCCCTCACCCCTCGCCCAGAATCTCAACCCCCTCGATCACCTGTTGCCCCACCCGCCGCCGCCAGCACCA‑CACATCCCCCCGCACAGCCCCAAGCCAAGCCCCCCCCAAAGCCGCACCCCACGTGCTCaa > 8:314272/1‑137 (MQ=255) cccAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAAc < 3:421707/137‑1 (MQ=255) cAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCCCCTCATCCCCCGGCCCCGCCCCGCGCCAAGCCCGCCCCCATGGCGCCCCCCAGGTGGACAAccc > 8:36792/1‑137 (MQ=255) cAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCCGCCCCGCCCCGCCCAAAGCCCGCCCCAATGGCGACCCCCCGGTGCTCCAccc > 2:408256/1‑137 (MQ=255) gggCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGt < 7:113497/137‑1 (MQ=255) gggCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGt < 5:92128/137‑1 (MQ=255) cccTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCCCCGGCCGCCGCCAC‑ACCACCTCATCCCCCCGCCCACCCCCCACCCAACCCCGCCCCAAAGGCGAACCCCAGTTGATCAACCCGTGcc > 3:162285/1‑137 (MQ=255) tCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCgg < 5:482305/137‑1 (MQ=255) cAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGc < 7:243969/137‑1 (MQ=255) cTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCg < 1:374593/137‑1 (MQ=255) cTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCg < 4:379505/137‑1 (MQ=255) cccAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑CCCACCTCATACCCCCGCCCCCCCCCGCGCAAAGCCCCCCCCAATGGCGCACCCCACTTGATAAACCCGGGCCCCGCCTCCCCCTc > 6:188700/1‑137 (MQ=255) cccAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCCGCCCCGCCCCCCGCAAATCCCCCCCCAAAGGCGCACCCCACTTGATCAAACCCTGCCCGGCAACCGccac > 3:11175/1‑135 (MQ=255) cAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCCTCCCCCGGCCCCGCCCCCAGCCAACCCCGCCCCCAAGGCGCACCCCAGTTGATCAACCCCTGCCCCGCCCCCGCcacaa > 2:367337/1‑136 (MQ=255) cAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAAAGGCGAACCCCAGTTGATAAACCCGGGCCCCGCATCCGCCTCAt > 2:45971/1‑137 (MQ=255) tctcAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑CCCACCTCATCCCCCCGCCCCCCCCCCAGCCAACCCCGCCCCACTAGCGAAACCCAGTTGATCAACCCGGGCCCCGCATACGCCTCAACATCg > 8:212303/1‑137 (MQ=255) catccaacTCTATCACCTGTTCCCCCACCGTTCTCCGCCAG‑ACCTCCTCATCCCCCGGCCCAGCCCCGAGCAAAGCACGCCCCAATGGCGAACGCCAGTTGATCAACCCGTGCCCTGCATCCGCCTCATCCTCGccc < 8:490287/134‑1 (MQ=255) caccTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCCGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGTACCCCAGTTGATAACACCGTGCCCCGCCTCCCCGGCATCCTCGCccaccc > 8:437768/1‑136 (MQ=255) tctATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCa > 2:469691/1‑54 (MQ=255) tctATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCa < 1:469691/54‑1 (MQ=255) tctATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCCGCCCCGCCCCGCGCAAACCCCCCCCCAAAGGCGCAACCCAATTGATAAACCCCTGCCCCGCAACCGCCACAACAGCGCCCACAAGcag > 8:48574/1‑135 (MQ=255) aTCACCTGTTCCCCCACCGGCCTCCGCCAG‑ACCAACTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGcac < 2:83693/137‑1 (MQ=255) aTCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGcac < 8:192113/137‑1 (MQ=255) tCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATTGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGcacc < 5:460488/137‑1 (MQ=255) tCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGcacc < 1:125899/137‑1 (MQ=255) tCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGcacc < 3:92367/137‑1 (MQ=255) cACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCACGGCCCAGCCCCGAGCAAATCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGcacca < 1:348380/137‑1 (MQ=255) ccTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCaccacc < 1:408256/137‑1 (MQ=255) gTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGc < 1:296287/137‑1 (MQ=255) cGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCt < 7:285291/137‑1 (MQ=255) cGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCCGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCt < 3:342573/137‑1 (MQ=255) cGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCACATCCCCGAGCAACGCCCTCCCCAATGGCGAACCCGAGTTGATCAACCCGTGCCCGGCATCCGACTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCt < 2:390885/137‑1 (MQ=255) ccccgccAG‑ACCACCTCATCCCCCCGCCCCGCCCCCACCCAAGCCCCCCCCAAAGGCCGACCCCCCTTGCACCACCCCAGCACCGCCACCCCCACAAGCTCCCCCCCCCGCGGGCCCCCCGGCTCCTGGGCCTgcgc > 4:257073/3‑137 (MQ=255) ccaccTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCTGCGCGTCGACAAAc < 7:36792/137‑1 (MQ=255) cTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCTGCGCGTCGACAAACCGCa < 3:313413/137‑1 (MQ=255) aaTCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCTGCGCGTCGACAAACCGCAcc > 5:203030/2‑137 (MQ=255) cccccGGCCCAGCCCCGAGCAAAGCCCTCCCCAATGGCGAAACCCAGTTGATGAACACTGTCCCGGCATCCTCCTCATCCTCGGCCTAGATCTGCACCAACTGCTCCTGGTACTGCGCGTCGACATACCTCACCCgg < 1:208445/137‑1 (MQ=255) cccccGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCTGCGCGTCGACAAACCGCACCCgg < 5:333450/137‑1 (MQ=255) ccccGGCCCAGCCACGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCTGCGCGTCGACAAACCGCACCCGGc < 1:259729/137‑1 (MQ=255) | CCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAG‑ACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCCGTGCCCGGCATCCGCCTCATCCTCGCCCACCAGCTGCACCACCTGCTCCTGGTCCTGCGCGTCGACAAACCGCACCCGGC > KT2440/783170‑783430 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |