Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,617,993 | C→T | 7.0% | Q340* (CAG→TAG) | PP_3188 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,617,993 | 0 | C | T | 7.0% | 164.2 / 2.1 | 57 | Q340* (CAG→TAG) | PP_3188 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (25/28); new base T (2/2); total (27/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01 |
GGTATTGTCGACCGTTCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCATGAGCGAAGCCCACGGCAAGGA > NC_002947/3617849‑3618141 | ggTATTGTCGACCGTTCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCg > 4:219738/1‑149 (MQ=255) tATTGTCGACCGTTCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGTTCATTGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTc < 7:226585/149‑1 (MQ=255) aTTGTCGACCGTTCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCt > 4:82651/1‑149 (MQ=255) gTCGACCGTTCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTgg < 3:219738/148‑1 (MQ=255) gACCGTTCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCATACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATAGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCt > 1:267277/1‑149 (MQ=255) cGTTCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCa > 8:305673/1‑149 (MQ=255) ttCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCAcc > 1:1025/1‑149 (MQ=255) cTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGCCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCaa > 5:325368/1‑149 (MQ=255) gTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGACCCTTGACGCAAACACGCCCATGCTTTCGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGTCCGCACACCGTCTGGGGATGCACCACAacac > 3:55536/1‑149 (MQ=255) gatgatTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCACCGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGat > 7:330705/1‑149 (MQ=255) atgatTTCGTCAGCAGCACCCTCGGCCCGCTGCTCGCCCGTGGACGCAAACACGCCCATCCTTTGGTGCATACCCTCGTGGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCAGCGCAACACGatc < 4:247862/149‑1 (MQ=255) tgatTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGatca < 3:82651/149‑1 (MQ=255) gTCAGCCGCACCCTCGGCCCGCTGTTCGCCCCTGGACGCAAACACGGCCATGCTTTGGTGCATACGCTCGACGCCTTGCTCATGGTATCTGGCAACGGGCGCAAGGCCGCACAGCGTCTGGGGCTGCACCg < 7:23103/131‑1 (MQ=255) gTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATcagc < 5:85458/148‑1 (MQ=255) gTCAGCCGCACCCTCGCCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGACTTGCTCATGGAAAATGGCAACGGGATCAAGGCCGCACAGCGTCTGGGGCTGCACCg > 8:23103/1‑131 (MQ=255) cGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGTAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATcagcgcatcca < 3:234352/149‑1 (MQ=255) cGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATcagcgcatc > 6:275084/1‑147 (MQ=255) cACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGAAACACGATCAATcagcgcatccagc > 3:85791/1‑149 (MQ=255) cctatgCCCTCTGATAGCCCCTGGACGCAAACACGCCCATGCTTTGGTTCATACCCTAGACGCCTTTCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCTTCTGGGGCTGCACCGCAACAAGAACAATcagcgcatccagcgca < 4:85791/144‑1 (MQ=255) ccTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATcagcgcatccagcgc > 4:107214/1‑148 (MQ=255) ccTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATcagcgcatccagcgc < 3:107214/148‑1 (MQ=255) cGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGaacaa > 5:292632/1‑149 (MQ=255) cGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGaacaa < 2:1025/149‑1 (MQ=255) cTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGaacaac < 8:285834/148‑1 (MQ=255) cTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGaacaac > 7:285834/1‑148 (MQ=255) aTCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAg < 3:5788/149‑1 (MQ=255) ccccTGGACGCATACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATAGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCggg < 2:267277/149‑1 (MQ=255) ccTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCggg > 4:222857/1‑147 (MQ=255) aaCACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGatcaat > 1:1974/1‑104 (MQ=255) aaCACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGatcaat < 2:1974/104‑1 (MQ=255) acGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCATAGCGTCTgggg < 5:293252/80‑1 (MQ=255) acGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCATAGCGTCTgggg > 6:293252/1‑80 (MQ=255) tttGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGat < 7:165151/87‑1 (MQ=255) tttGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGat > 8:165151/1‑87 (MQ=255) tGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTc < 5:191284/147‑1 (MQ=255) gTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGGATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCa < 6:292632/148‑1 (MQ=255) gTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCAt < 5:275084/149‑1 (MQ=255) ccTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCATAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCtt > 2:10992/1‑149 (MQ=255) gACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCg < 3:222857/147‑1 (MQ=255) cGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCg < 7:305673/149‑1 (MQ=255) aTGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCgg > 7:69431/1‑136 (MQ=255) aTGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCgg < 8:69431/136‑1 (MQ=255) aTGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTgg > 2:265858/1‑148 (MQ=255) aaaTGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGgc < 1:265858/149‑1 (MQ=255) gcatctggCTCAAGGGTGCACAGCGTCTGGGGCTGCACCAGTTCAGGATCTATGAGCGCATCCAGCGCATCGAACATCTGATCGGGCTGTGGATGGACGACCCGCTTTTTCGAATGAATGCTTCGGTCGCCATTCTGGTATGGCGGAt < 4:55536/142‑1 (MQ=255) gCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCAt > 1:335750/1‑148 (MQ=255) aaCGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGaacaac < 5:327694/76‑1 (MQ=255) aaCGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGaacaac > 6:327694/1‑76 (MQ=255) aCGGGCTCAAGGCCGCATAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCATGAGc < 1:10992/149‑1 (MQ=255) ggCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCATGAGCGaa < 6:325368/149‑1 (MQ=255) gCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCtt < 6:306876/100‑1 (MQ=255) gCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCtt > 5:306876/1‑100 (MQ=255) tCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGgcgc < 6:176947/137‑1 (MQ=255) tCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGgcgc > 5:176947/1‑137 (MQ=255) aaGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGAGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCATGAGCGAAGCCCa > 5:251468/1‑149 (MQ=255) ggCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCATGAGCGAAGCCCACg < 2:335750/149‑1 (MQ=255) cGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCATGAGCGAAGCCCACGGc > 8:36043/1‑148 (MQ=255) caGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCATGAGCGAAGCCCACGGCAAGGa < 7:36043/149‑1 (MQ=255) | GGTATTGTCGACCGTTCGCTGGTCGATGATTTCGTCAGCCGCACCCTCGGCCCGCTGATCGCCCCTGGACGCAAACACGCCCATGCTTTGGTGCATACCCTCGACGCCTTGCTCATGGAAAATGGCAACGGGCTCAAGGCCGCACAGCGTCTGGGGCTGCACCGCAACACGATCAATCAGCGCATCCAGCGCATCGAACAACTGAGCGGGCAGTCGCTGGACGACCCGCTTTTTCGCATGAATGCTTCGGTCGCCATGCTGGTATGGCGCATGAGCGAAGCCCACGGCAAGGA > NC_002947/3617849‑3618141 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |