Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,827,811 | G→C | 14.8% | G223A (GGG→GCG) ‡ | mtgA → | monofunctional biosynthetic peptidoglycan transglycosylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,827,811 | 0 | G | C | 14.8% | 69.2 / 6.3 | 27 | G223A (GGG→GCG) ‡ | mtgA | monofunctional biosynthetic peptidoglycan transglycosylase |
Reads supporting (aligned to +/- strand): ref base G (16/7); new base C (2/2); total (18/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.82e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01 |
TCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGG > NC_002947/5827667‑5827960 | tCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCa < 1:343031/150‑1 (MQ=255) ccGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGAt < 2:155146/150‑1 (MQ=255) cAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCGCCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCAGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCt > 1:449875/1‑150 (MQ=255) cagcgcagcTGGCAGCGGTGCTGCCAAGCCCACTGGATGGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGACGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGaa > 1:356207/1‑150 (MQ=255) tgcCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccggc > 1:50346/1‑150 (MQ=255) aaGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCtt > 2:360911/1‑150 (MQ=255) gCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTc > 2:366512/1‑150 (MQ=255) gTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGACCCTTTCGCGGGTAAAc > 1:309624/1‑150 (MQ=255) tGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAAt < 2:483757/122‑1 (MQ=255) tGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAAt > 1:483757/1‑122 (MQ=255) cAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTa > 1:393850/1‑150 (MQ=255) gaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGCTGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGc > 2:86802/1‑150 (MQ=255) gaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCtt < 1:102870/104‑1 (MQ=255) gaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCtt > 2:102870/1‑104 (MQ=255) cgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCACTCCTACAAGCaaa > 2:401062/1‑149 (MQ=255) tGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGc < 2:302721/150‑1 (MQ=255) gCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcct > 1:513219/1‑150 (MQ=255) cgggccggCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTcgc > 1:158593/1‑150 (MQ=255) gCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCaaa > 1:236491/1‑150 (MQ=255) ggATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAAccc > 2:154542/1‑150 (MQ=255) tCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCgcgc < 2:236491/150‑1 (MQ=255) atcaggtaGCTGGGCCCGCCCAGCTGCCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCaa > 1:201147/8‑150 (MQ=38) tcaggtaGCTGGGCCCGCCCAGCTGGCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCcgc < 1:187486/131‑1 (MQ=38) tcaggtaGCTGGGCCCGCCCAGCTGGCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCcgc > 2:187486/7‑137 (MQ=38) tcaggtaGCTGGGCCCGCCCAGCTGCCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAAc < 2:201147/144‑1 (MQ=39) gAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAAcc < 1:154542/150‑1 (MQ=255) tGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGGAAACCCGCGCCTACAGGCAACCGCATcccc > 1:120513/1‑150 (MQ=255) cgggcCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACgg > 2:46821/1‑87 (MQ=255) cgggcCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACgg < 1:46821/87‑1 (MQ=255) ggcCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAgg > 1:229838/1‑150 (MQ=255) | TCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGG > NC_002947/5827667‑5827960 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |