Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,087,376 | C→G | 23.3% | intergenic (+146/+772) | PP_1865 → / ← aroH | transposase/3‑deoxy‑7‑phosphoheptulonate synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,087,376 | 0 | C | G | 23.3% | 59.9 / 16.5 | 30 | intergenic (+146/+772) | PP_1865/aroH | transposase/3‑deoxy‑7‑phosphoheptulonate synthase |
Reads supporting (aligned to +/- strand): ref base C (7/16); new base G (4/3); total (11/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.72e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01 |
TTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATC > NC_002947/2087227‑2087523 | ttGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCccatc < 1:449079/150‑1 (MQ=255) gAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCAAGCGCTGCAGGCCccatcca > 1:44387/1‑150 (MQ=255) aGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATgg < 1:702738/150‑1 (MQ=255) ggATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCCGCGCTCCGGGGTCCCCTCCCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATgggg < 1:725189/150‑1 (MQ=255) gCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTG‑AGGCCCCATGGATGGGGCCtgca < 1:120880/150‑4 (MQ=17) ttACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAgcgc < 2:177714/150‑1 (MQ=255) cggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACg < 2:696758/150‑1 (MQ=255) gTCGTGGCAAAAAGGAACCAAGAACCGTCGGGTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGccc < 2:155087/150‑1 (MQ=255) gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTTCCCTCGCTCCGGCCTTGCGCCCG‑GGAGCACGCGCTGCAGGCCCCACCCATGGGGCCT‑CAGCGCTTGACGGGCATGGATGCCCGTc > 1:391396/1‑150 (MQ=17) gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATGGATGCCCGTc < 1:92497/150‑1 (MQ=255) gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCA‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTc < 1:598450/150‑1 (MQ=255) tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCt < 1:659707/121‑1 (MQ=255) tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCt > 2:659707/1‑121 (MQ=255) gCAAAAAGGAACCAAACACCGGCGGCTCCCATCATCCGGCCCCTAGGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATGGATGCCCGTCaag < 2:391396/150‑3 (MQ=17) gCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATgg < 2:461294/114‑1 (MQ=255) gCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATgg > 1:461294/1‑114 (MQ=255) aaaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATg > 2:449048/1‑150 (MQ=255) aaaCCGTCGGCTCCCATCATCCGGCACCTACGCGTCTCTCCGGCGCCCCCTCCCTCCGGCC‑ACCTGCCGTCGGGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGcc < 1:444697/150‑1 (MQ=17) catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc < 2:44387/150‑1 (MQ=255) atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATGGATGGGGCCtgc < 2:85168/91‑3 (MQ=18) atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATGGATGGGGCCtgc > 1:85168/1‑89 (MQ=17) cccTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCaa > 1:380095/1‑150 (MQ=255) cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGc < 2:569815/150‑1 (MQ=255) ggggTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGAc < 2:635516/150‑1 (MQ=255) cTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGACGTCGCAATCTGCGTCGCCTGGACTACcgcgc > 1:90132/1‑150 (MQ=255) cG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTg > 1:493870/1‑120 (MQ=255) cG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTg < 2:493870/120‑1 (MQ=255) cgcTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGc < 1:396095/150‑1 (MQ=255) gctgaGGCCCCATGGATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGa > 2:706758/5‑150 (MQ=255) tgaGGCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGaa > 2:109861/3‑87 (MQ=18) tgaGGCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGaa < 1:109861/85‑1 (MQ=18) tgaGGCCCCATGGATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAt > 1:90050/3‑150 (MQ=255) ctccatccatGGGACCG‑TATTGCTTAACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGt < 2:563888/148‑1 (MQ=21) catccatGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGag > 1:563888/1‑149 (MQ=255) atccatGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATc > 1:531889/1‑150 (MQ=255) | TTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCG‑GGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATC > NC_002947/2087227‑2087523 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |