Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,823,683 | C→T | R428R (CGC→CGT) | spoT → | bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,823,683 | 0 | C | T | 94.3% | 100.8 / ‑0.7 | 35 | R428R (CGC→CGT) | spoT | bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase |
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (22/11); total (23/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01 |
CTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGT > NC_000913/3823548‑3823812 | ctTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAg > 2:107787/1‑139 (MQ=255) ttCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGc > 2:764805/1‑139 (MQ=255) gATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTAcc > 2:37648/1‑139 (MQ=255) agaTTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCt > 1:1085901/1‑139 (MQ=255) tACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGc < 2:982892/139‑1 (MQ=255) aCGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCa > 1:954135/1‑139 (MQ=255) ttCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGc > 1:424552/1‑139 (MQ=255) tGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGtt < 2:935425/139‑1 (MQ=255) tGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCg > 2:417785/1‑139 (MQ=255) cAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCggg > 1:358668/1‑139 (MQ=255) aCGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGgcg < 2:424552/139‑1 (MQ=255) gCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGc > 1:1142926/1‑89 (MQ=255) gCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGc < 2:1142926/89‑1 (MQ=255) aCTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGaa < 1:107787/139‑1 (MQ=255) gCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACTAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCg < 2:1085901/139‑1 (MQ=255) ttATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCt < 1:417785/139‑1 (MQ=255) cATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACt > 1:1215755/1‑139 (MQ=255) aCCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTg > 1:433382/1‑139 (MQ=255) atatCGGTCATGCCTGCGTGGACGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCTGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGt < 2:721654/139‑1 (MQ=255) tatCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt > 2:694768/1‑139 (MQ=255) tatCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGtt > 2:776649/1‑139 (MQ=255) tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCt < 2:465097/97‑1 (MQ=255) tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCt > 1:465097/1‑97 (MQ=255) tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaaa > 1:919316/1‑139 (MQ=255) tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaaa > 2:1074316/1‑139 (MQ=255) tGCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaaa > 1:1212719/1‑139 (MQ=255) gCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCt > 2:261360/1‑116 (MQ=255) gCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCt < 1:261360/116‑1 (MQ=255) gCCTGCGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAg > 1:1052510/1‑139 (MQ=255) cGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgc > 2:1044879/1‑139 (MQ=255) cGTGGGCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTAAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgc > 1:771450/1‑139 (MQ=255) cGTGGGCGCACGCGTTGAACGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGACGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgc > 1:1011002/1‑139 (MQ=255) ggCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCa > 1:1071811/1‑139 (MQ=255) ggCGCACGCGTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCa < 1:1021165/139‑1 (MQ=255) cgTTGACCGTCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGt < 1:764805/139‑1 (MQ=255) | CTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGT > NC_000913/3823548‑3823812 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |