Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,312,732 | T→A | 9.7% | intergenic (+30/‑132) | PP_2912 → / → hemB | hypothetical protein/5‑aminolevulinate dehydratase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,312,732 | 0 | T | A | 9.7% | 156.1 / 8.7 | 61 | intergenic (+30/‑132) | PP_2912/hemB | hypothetical protein/5‑aminolevulinate dehydratase |
Reads supporting (aligned to +/- strand): ref base T (24/30); new base A (2/4); total (26/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.89e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01 |
GGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAA > NC_002947/3312585‑3312871 | ggTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑c < 1:1004138/150‑1 (MQ=255) tGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑ccc > 2:110753/1‑150 (MQ=255) cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGaa < 2:1581941/150‑1 (MQ=255) taagaTCGACGCCGCCCGGATGGACGACGCTCGTCACCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAaga > 1:1621391/3‑150 (MQ=255) tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAaga < 1:133859/150‑1 (MQ=255) aTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGcc > 2:542147/1‑150 (MQ=255) gccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAAGAGGcggtacag < 2:1405354/150‑8 (MQ=38) ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAg < 1:711189/150‑1 (MQ=255) ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAg < 2:810655/150‑1 (MQ=255) ccGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAAGAGGcg > 2:926532/1‑138 (MQ=255) ccGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGC‑GGGTAACCCCGCGAAGAGGcg < 1:926532/139‑2 (MQ=255) cGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGc > 1:230383/1‑137 (MQ=255) cGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGc < 2:230383/137‑1 (MQ=255) aTGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCa < 1:110753/150‑1 (MQ=255) acgacgCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGc > 2:1122104/1‑150 (MQ=255) tcgtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTTCTGTACCGCCTCTTCGCGGGGGTA‑CCCGCGAAGACGCCAACACAGTCAAACCACGGCCGTTGcc > 1:300729/1‑150 (MQ=255) cgtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCa < 1:50080/150‑1 (MQ=255) gCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGa < 2:336468/150‑1 (MQ=255) gcagcaCGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCaa < 1:434205/150‑1 (MQ=255) acgcCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACAt > 2:1127527/2‑150 (MQ=255) ggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTc > 2:255340/1‑150 (MQ=255) cgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCAtt > 1:1280551/1‑150 (MQ=255) gggCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtt < 1:711747/150‑1 (MQ=255) gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGa < 1:93669/115‑1 (MQ=255) gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGa > 2:93669/1‑115 (MQ=255) ctaaGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTg < 1:108204/148‑1 (MQ=255) ccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGCAGAGGCCAACACGGGCAAACAACGGCCGTTGCTAGTGACGCAAAGCACAATCTCGTCATTTAGCCCTTTCACTg > 1:275448/1‑150 (MQ=255) cAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTCTATTAAAGCTGGGCATCTACACGAGGTTA‑CACGGTAAGAGGCCAACACAGGCAAACCACGGCCGTTGACATTGACGCAAGACACATTCTCGTCAGTTAGCCTTTTCACTgg > 1:444612/1‑150 (MQ=255) aaGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGTGATAA‑GGCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGgc < 1:1349576/150‑1 (MQ=255) agagGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGC‑GGGTAACCCCGCGAAGAGGcg > 2:1501693/1‑75 (MQ=25) agagGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGC‑GGGTAACCCCGCGAAGAGGcg < 1:1501693/76‑2 (MQ=25) ggaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa < 1:542147/150‑1 (MQ=255) aaggTAACGCCGCCTCGGCTGAAGTAGTTTATTCCTGTTCCGCCTCTTCTCGGGGTTA‑CCTGCGAAGAGGCCAACACTGGCCAACCACGGCCCTTTACATTGACGCGTGATACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGCCAg < 2:275448/149‑1 (MQ=255) ccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt < 1:255340/150‑1 (MQ=255) cgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGta < 1:844811/150‑1 (MQ=255) gccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGtat < 1:752561/150‑1 (MQ=255) ctagCTGCAGTAGGTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACg > 2:1227411/4‑150 (MQ=255) gCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTc < 2:1240230/150‑1 (MQ=255) tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGac > 2:338152/1‑87 (MQ=255) tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGac < 1:338152/87‑1 (MQ=255) tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAAt > 2:174919/1‑150 (MQ=255) tATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCAt > 1:1017990/1‑71 (MQ=255) tATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCAt < 2:1017990/71‑1 (MQ=255) aTTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTg > 2:52551/1‑150 (MQ=255) tCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa < 2:1293336/150‑1 (MQ=255) cTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAg < 1:539190/116‑1 (MQ=255) cTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAg > 2:539190/1‑116 (MQ=255) tGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt > 2:80473/1‑150 (MQ=255) aCCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTc < 1:732999/150‑1 (MQ=255) ccGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg > 2:568540/1‑118 (MQ=255) ccGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg < 1:568540/118‑1 (MQ=255) cGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa > 2:1390349/1‑142 (MQ=255) cGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa < 1:1390349/142‑1 (MQ=255) cGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCAt < 1:80473/150‑1 (MQ=255) ccTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGc < 2:423079/110‑1 (MQ=255) ccTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGc > 1:423079/1‑110 (MQ=255) ctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCaa > 2:678496/1‑127 (MQ=255) ctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCaa < 1:678496/127‑1 (MQ=255) ttCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa < 2:1280551/150‑1 (MQ=255) ggggTTA‑CCCTCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa > 2:23227/1‑145 (MQ=255) ggggTGA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa < 1:23227/145‑1 (MQ=255) | GGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAA > NC_002947/3312585‑3312871 |
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: ‑ |