Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,286,493 | T→G | 11.8% | N154K (AAT→AAG) | PP_2888 → | ECF family RNA polymerase sigma factor PrtI |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,286,493 | 0 | T | G | 11.8% | 78.2 / 2.5 | 33 | N154K (AAT→AAG) | PP_2888 | ECF family RNA polymerase sigma factor PrtI |
Reads supporting (aligned to +/- strand): ref base T (12/16); new base G (2/2); total (14/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCCGACCTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTGCAG > NC_002947/3286344‑3286633 | gCCGACCTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAAt > 2:385666/1‑150 (MQ=255) ccTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCtt < 2:443149/150‑1 (MQ=255) ttGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGa < 2:589703/150‑1 (MQ=255) ccGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGa > 2:193266/1‑150 (MQ=255) aCGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGcaac < 1:536933/150‑1 (MQ=255) gcggAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGcaacccc < 2:539242/150‑1 (MQ=255) agcGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGcaacccccaacccc < 1:193266/150‑1 (MQ=255) tGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGc < 1:267797/150‑1 (MQ=255) cATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGAcg > 1:50847/1‑150 (MQ=255) gAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGc > 2:347631/1‑150 (MQ=255) aaGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCt > 1:447649/1‑150 (MQ=255) aGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGACTg > 1:606113/1‑150 (MQ=255) cAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgat < 1:347631/150‑1 (MQ=255) aggaggCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgatga > 1:392494/1‑150 (MQ=255) ggCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGACACTGaaa < 1:426606/121‑1 (MQ=255) ggCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGACACTGaaa > 2:426606/1‑121 (MQ=255) gccgccCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAAt < 2:261665/150‑1 (MQ=255) ttGGGCATCCCCATCGGCACTGTGACGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgatga < 1:153714/133‑1 (MQ=255) ttGGGCATCCCCATCGGCACTGTGACGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgatga > 2:153714/1‑133 (MQ=255) tGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTtcgt < 1:456661/150‑1 (MQ=255) aCTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAAc < 2:50847/150‑1 (MQ=255) tgtgATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGgaga > 2:189833/1‑82 (MQ=255) tgtgATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGgaga < 1:189833/82‑1 (MQ=255) tgtgATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCt < 2:447649/150‑1 (MQ=255) tgATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCtgt > 2:250008/1‑150 (MQ=255) gcCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGc > 1:67983/1‑150 (MQ=255) tCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGTGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAgg < 1:385666/150‑1 (MQ=255) tCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAgg > 2:277182/1‑150 (MQ=255) cgcgcgcgcgCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGt < 2:36327/150‑1 (MQ=255) cgcgcgccccCAACGCCTTGCGCGCACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGt < 2:36298/150‑1 (MQ=255) ccgatctcgcAAGGCATTGCGCGAACTGACCGAGGGCAACCCCCAAAACCCGTCCGTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTg < 2:418008/143‑1 (MQ=255) gttcgcgcAAGGCATTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTg < 2:10624/121‑1 (MQ=255) gttcgcgcAAGGCATTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGAACCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTg > 1:10624/4‑124 (MQ=255) ttcgcgcAAGGCATTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCAAGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTGCAg > 2:222973/3‑150 (MQ=255) | GCCGACCTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTGCAG > NC_002947/3286344‑3286633 |
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: ‑ |