Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,899,213 | Δ1 bp | 9.4% | coding (1395/1413 nt) | PP_2552 → | DOPA decarboxylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,899,213 | 0 | G | . | 9.4% | 159.4 / 12.1 | 53 | coding (1395/1413 nt) | PP_2552 | DOPA decarboxylase |
Reads supporting (aligned to +/- strand): ref base G (26/22); new base . (2/3); total (28/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01 |
ATACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGC > NC_002947/2899065‑2899362 | aTACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAgg > 2:664518/1‑150 (MQ=255) tACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGa > 2:963837/1‑150 (MQ=255) cAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGt > 2:236208/1‑150 (MQ=255) aaGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCCTGCA‑GACGTg < 2:1270309/150‑1 (MQ=255) aGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCCTGCA‑GACGTgc < 1:899887/150‑2 (MQ=255) gggCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGCGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCCGTGGGCACGT‑CTGCAGGACGTGAt > 2:698874/1‑150 (MQ=255) gggCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAgg > 1:1420775/1‑150 (MQ=255) cGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTg > 2:498077/1‑150 (MQ=255) cGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGc < 2:1261814/150‑1 (MQ=255) aTCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAggg > 2:129539/1‑132 (MQ=255) aTCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAggg < 1:129539/132‑1 (MQ=255) cgcTGCTGATGTGACCGCGGCTACACTGGACGGGCGGTTGATGGTGCGGGTTTCGATTGGTGCGCTGCCGCCCGAGCGGAGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGCGCTGAAGCCTGTACGGTACCTGttt < 1:620446/150‑3 (MQ=255) gctgctTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATg > 2:753505/1‑150 (MQ=255) ttATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCgg < 2:1420775/150‑1 (MQ=255) cTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGTTTCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAg < 2:1045428/150‑1 (MQ=255) tACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAgg < 1:236208/150‑1 (MQ=255) tACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAgg < 1:963837/150‑1 (MQ=255) acTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGc > 2:158034/1‑150 (MQ=255) gACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCa > 1:1299823/1‑150 (MQ=255) gggCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGctttg < 1:158034/150‑5 (MQ=255) ggCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAggg > 2:18878/1‑140 (MQ=255) ggCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAggg < 1:18878/140‑1 (MQ=255) tGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGCCCCttc < 1:498077/150‑3 (MQ=255) tGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCGGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAAt > 1:1265880/1‑150 (MQ=255) cGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCt > 2:1345161/1‑150 (MQ=255) tCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGt < 2:1299823/150‑1 (MQ=255) aTTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCaa < 2:758339/138‑1 (MQ=255) aTTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCaa > 1:758339/1‑138 (MQ=255) ggTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATc < 2:39046/63‑1 (MQ=255) ggTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATc > 1:39046/1‑63 (MQ=255) tGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagcag > 2:86803/1‑124 (MQ=255) tGCGCTGCCGACCGAGCGGGGGGATGGGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAagcag < 1:86803/124‑1 (MQ=255) cgaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGACGGGGCTGGGAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaaa < 2:15391/150‑1 (MQ=255) gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCTTTGCAGccc > 1:646534/1‑118 (MQ=255) gaccgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAg > 2:511724/1‑150 (MQ=255) gaccgaGCGGGGGGATGTGCAGCGGCTGTGCGCAGGC‑CTGCAGGGCGCGGGTTATGTCTAAAGCCTGTAGGGGTGCGGTATGAGCGGGCTTGCCCGCGAAGGGGCGGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAg > 1:26594/1‑150 (MQ=255) gaccgaGCGGGGGGATGTGAAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCAGTGCAATCTATCAGTAGCTCACCACCTGAAAg > 1:164551/1‑150 (MQ=255) cgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTa > 1:283845/1‑150 (MQ=255) cGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCtt > 2:570653/1‑150 (MQ=255) gggggATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTa > 2:1240350/1‑150 (MQ=255) cAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc < 2:1407441/135‑1 (MQ=255) cAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc > 1:1407441/1‑135 (MQ=255) cAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCtt > 1:735881/1‑150 (MQ=255) aGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTc < 1:664518/150‑1 (MQ=255) ggCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATa < 1:698874/150‑1 (MQ=255) ggCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATa < 1:753505/150‑1 (MQ=255) ttgatcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGCGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCt > 2:1320155/6‑150 (MQ=255) gggCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGc < 2:283845/150‑1 (MQ=255) tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc > 1:1219624/2‑145 (MQ=255) tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGc < 2:1219624/144‑1 (MQ=255) gCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCt < 1:873543/114‑1 (MQ=255) gCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCt > 2:873543/1‑114 (MQ=255) gCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCt > 1:910832/1‑150 (MQ=255) cACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCtt < 2:735881/150‑1 (MQ=255) ggACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGc > 2:297783/1‑150 (MQ=255) | ATACCAAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGCCTTTTCCATCGGC > NC_002947/2899065‑2899362 |
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: ‑ |