Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,087,377 | C→G | 8.9% | intergenic (+147/+771) | 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,377 | 0 | C | G | 8.9% | 135.7 / 4.1 | 45 | intergenic (+147/+771) | PP_1865/aroH | transposase/3‑deoxy‑7‑phosphoheptulonate synthase |
Reads supporting (aligned to +/- strand): ref base C (21/20); new base G (2/2); total (23/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 |
AAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGG > NC_002947/2087230‑2087525 | aaaGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccatccat < 1:223838/150‑1 (MQ=255) aGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg > 1:1263901/1‑150 (MQ=255) aGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCAATCCATgg < 2:709529/150‑1 (MQ=255) ccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAgc > 1:326268/1‑150 (MQ=255) gggTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCCCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACgg < 1:1374466/150‑1 (MQ=255) tCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCa < 2:326268/150‑1 (MQ=255) cTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg < 2:267392/126‑1 (MQ=255) cTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg > 1:267392/1‑126 (MQ=255) tGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCa > 1:67673/1‑150 (MQ=255) ggCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATGGATGCCCGTCaa < 1:136496/150‑2 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGc > 1:546835/1‑140 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGc < 2:546835/140‑1 (MQ=255) aaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCN‑CAGCGCNTGACGGGCATCCNTGCCCGTCATCTCCCTCCGCAATGCCt < 1:1379658/150‑1 (MQ=255) aCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTc < 1:1181640/150‑1 (MQ=255) cGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCa > 2:1140041/1‑150 (MQ=255) ggCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCg < 2:462422/150‑1 (MQ=255) tCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcc > 1:1116782/1‑150 (MQ=255) ccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctc > 2:396021/1‑150 (MQ=255) ccATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctc > 1:964990/1‑150 (MQ=255) catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCCTGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc < 1:1355520/150‑1 (MQ=255) catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGcctcc < 2:67673/150‑1 (MQ=255) cGGCCCCTACGCTGGGCTCGGGGGTCCCCTCCCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTc < 2:221086/150‑1 (MQ=255) tACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCAt < 2:752825/97‑1 (MQ=255) tACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCAt > 1:752825/1‑97 (MQ=255) tGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGGCGCAATCTGCGTc > 1:1212693/1‑150 (MQ=255) cgcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGc < 1:1005046/150‑1 (MQ=255) ggTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCaagc > 2:912480/1‑90 (MQ=255) ggTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCaagc < 1:912480/93‑4 (MQ=255) cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgc > 1:596571/1‑150 (MQ=255) cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa > 1:821239/1‑150 (MQ=255) cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGa > 1:483646/1‑150 (MQ=255) ccTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGc < 1:396021/150‑1 (MQ=255) cTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGAt < 2:964990/150‑1 (MQ=255) tCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATc > 1:938851/1‑150 (MQ=255) ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATcc < 1:157118/56‑1 (MQ=38) ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATcc > 2:157118/1‑56 (MQ=38) ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaa > 1:654668/1‑150 (MQ=255) aCGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTACGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCgaga > 2:1023834/1‑150 (MQ=255) gcgcTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAg > 2:985096/1‑150 (MQ=255) cgcTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGc < 2:1263901/150‑1 (MQ=255) gaGGCCCCATGGATGGGACCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTc < 2:248188/99‑1 (MQ=17) gaGGCCCCATGGATGGGACCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTc > 1:248188/2‑100 (MQ=17) ggCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTct < 1:985096/150‑2 (MQ=255) gCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCaa > 1:1241295/1‑90 (MQ=17) gCCCCATGGATGGGGCCTGCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCaa < 2:1241295/90‑1 (MQ=17) cccATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTg > 2:1387055/1‑150 (MQ=255) ccatccatGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGa > 2:799936/1‑150 (MQ=255) ccatccatGCGGCCT‑GAGCTCTTGACGGGCATCCATGCCCGTCATCTCCCCCCTGCATGCCGCCACTCGGCTTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAATGTCAAAAGCGAGAAGCGAATTCGGTGa < 2:225180/150‑1 (MQ=18) atccatGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATc < 1:1140041/150‑1 (MQ=255) atccatGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATGCGGTGATc > 1:1129717/1‑150 (MQ=255) ccatGGGGCCT‑CAGCGCTTGTCGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcgg < 1:384259/150‑1 (MQ=255) | AAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCT‑CAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGG > NC_002947/2087230‑2087525 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |