Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,528,612 | T→A | 18.4% | L129Q (CTG→CAG) | dnaEB → | error‑prone DNA polymerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,528,612 | 0 | T | A | 18.4% | 90.5 / 13.4 | 38 | L129Q (CTG→CAG) | dnaEB | error‑prone DNA polymerase |
Reads supporting (aligned to +/- strand): ref base T (19/12); new base A (2/5); total (21/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.07e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGC > NC_002947/3528467‑3528750 | gACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGcc > 2:171093/1‑150 (MQ=255) ggCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTtgt < 2:465058/150‑1 (MQ=255) cccAAGCTGGTGCTGCTGGTCTAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtg < 2:144562/150‑4 (MQ=255) ccAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtgg < 1:214891/150‑5 (MQ=255) cAAGCTGGTGCTGCTGGTCGAAAACCTGCCTGTTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACTAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTgtgggt < 1:585464/150‑1 (MQ=255) cAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCCACAAGGCCAGCAGGCCTtggtggt < 2:33329/150‑6 (MQ=255) cAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCAGCAGGCCTtggtggt < 1:1014/150‑6 (MQ=255) tGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCg > 1:605487/1‑150 (MQ=255) gctgctGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGa > 2:505901/1‑150 (MQ=255) tgGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTg < 1:427197/150‑1 (MQ=255) gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGAt > 1:546730/1‑150 (MQ=255) cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTgg < 2:605487/150‑1 (MQ=255) cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTgg > 1:22877/1‑150 (MQ=255) gAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCt < 1:171093/150‑1 (MQ=255) aCCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTgc < 1:448828/150‑1 (MQ=255) cTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTgcgc > 2:137225/1‑150 (MQ=255) cTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTgcgc > 1:248388/1‑150 (MQ=255) cTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTc > 1:537580/1‑150 (MQ=255) ccGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTgg < 1:452131/105‑1 (MQ=255) ccGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTgg > 2:452131/1‑105 (MQ=255) cGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAgcgc < 2:22877/150‑1 (MQ=255) cggcggGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGtggctg < 1:137225/150‑1 (MQ=255) atagaAGTGGGCGTACTAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCa > 2:469076/3‑150 (MQ=255) ggCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCg < 2:537580/150‑1 (MQ=255) tGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACg > 1:397617/1‑150 (MQ=255) tttCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGa > 1:528033/1‑107 (MQ=255) tttCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGa < 2:528005/107‑1 (MQ=255) tttCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGa < 2:528033/107‑1 (MQ=255) tttCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGa > 1:528005/1‑107 (MQ=255) ttGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTg > 1:343847/1‑150 (MQ=255) ttGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTg > 2:458966/1‑150 (MQ=255) ctgcACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACgg > 2:460166/1‑150 (MQ=255) tgcACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGt > 2:392036/1‑150 (MQ=255) ccaccaAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCg > 1:227736/1‑150 (MQ=255) ccacaAGGCCAGCAGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCaa > 2:533460/4‑141 (MQ=255) ccacaAGGCCAGCAGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTCCGCATTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCaa < 1:533460/138‑1 (MQ=255) caccaAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGGGCGTGGCTGCGCAGTGTGTTGGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGAGGAAAGAACGATTCCGTGCGCCTGCAACGGTTGCgg > 2:232517/1‑150 (MQ=255) acaAGGCCAGCAGGCCTTGGTGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTCGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGc > 1:402137/2‑150 (MQ=255) | GACGGGCCCAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGTGTTCGCTGAGCGCCTGTGGCTGGCCGTGCACCTGCATCGGGGCAGTGACGATGCCGTGCGCCTGCAACGGTTGCGGGC > NC_002947/3528467‑3528750 |
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: ‑ |