Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,555,429:1 | +C | 13.3% | intergenic (+25/+49) | sbcB → / ← PP_1366 | exodeoxyribonuclease I/transcriptional regulator MvaT |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,555,429 | 1 | . | C | 13.3% | 143.9 / 21.0 | 45 | intergenic (+25/+49) | sbcB/PP_1366 | exodeoxyribonuclease I/transcriptional regulator MvaT |
Reads supporting (aligned to +/- strand): ref base . (21/18); new base C (4/2); total (25/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01 |
CTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATC > NC_002947/1555281‑1555574 | cTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCa‑a > 2:899565/1‑150 (MQ=255) cgcCAGCAAGCCTGGGAGGGCGCCTACGAGGCAGGTAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGc < 2:192971/150‑1 (MQ=255) cgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCCGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACGTGTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGc < 1:73663/150‑1 (MQ=255) cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGACGCAAGAgg > 2:100159/1‑150 (MQ=255) cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg < 1:363789/150‑1 (MQ=255) gacgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATa > 1:830468/1‑150 (MQ=255) acgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa < 1:702716/150‑1 (MQ=255) aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa < 2:830468/150‑1 (MQ=255) gcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGt > 1:975515/1‑150 (MQ=255) caggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTc < 1:485733/150‑1 (MQ=255) ggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCACTCGGCTGACACCCAATAAACAAAAACGCCGGTA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCca < 1:134681/150‑1 (MQ=255) aggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagca < 2:975515/150‑1 (MQ=255) tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACAT‑TGCCGGCGttt > 2:406437/1‑87 (MQ=38) tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑CATTGCCGGCGttt < 1:406437/87‑1 (MQ=38) tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa < 1:1054519/150‑1 (MQ=255) tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa < 1:1054390/150‑1 (MQ=255) tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc > 2:621453/1‑150 (MQ=255) aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCaa > 2:299231/1‑150 (MQ=255) aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCaa > 1:833204/1‑150 (MQ=255) aaCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCa < 1:212088/150‑1 (MQ=255) cTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCActct > 2:914783/1‑150 (MQ=255) ttCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGCCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTc < 1:131576/150‑1 (MQ=255) aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg > 1:403643/1‑150 (MQ=255) aaTCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg < 1:645416/150‑1 (MQ=255) aTCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAggg > 2:937422/1‑150 (MQ=255) aTCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCAATCTTTCatgt > 2:630957/1‑147 (MQ=255) tGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGAGTTGACCAGCTTACAACCACTTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgt > 1:1098123/1‑150 (MQ=255) aTAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccgc > 2:205847/1‑150 (MQ=255) aaCAAAAACGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCtt > 1:348650/1‑150 (MQ=255) aCAAAAACGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt > 2:784777/1‑150 (MQ=255) aCAAAAACGCCGGCA‑CATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt > 2:787042/1‑150 (MQ=255) aCAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt > 1:414793/1‑150 (MQ=255) caacaaCCCCCGCA‑ACGCGCCGGCGTTTTTGTTGTCCGCAGCAAGAGGGCGTATGCGGATAATCAGGCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAACTGGCTTTCCTCTCATTCAGGGTCTTGTGGTTGCCGCCTTTg < 2:414793/150‑1 (MQ=255) cAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTg > 2:624401/1‑150 (MQ=255) aaaaaCGCCGGCA‑CATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTgg > 1:837055/1‑150 (MQ=255) aaaaCGCCGGCA‑ATGTGACACGGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGt < 1:899565/150‑1 (MQ=255) aaaCGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTc > 2:359459/1‑95 (MQ=38) aaaCGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTc < 1:359459/95‑1 (MQ=38) aaCGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt < 2:403643/150‑1 (MQ=255) aCGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAg < 2:887196/61‑1 (MQ=25) aCGCCGGCA‑CATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAg > 1:887196/1‑61 (MQ=25) cGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCg > 2:654041/1‑150 (MQ=255) gCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa > 2:605409/1‑150 (MQ=255) gCCGGCA‑ATGTGCCGGCGTTTTTATTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGTTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa < 1:860269/150‑1 (MQ=255) cGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATc > 1:14636/1‑150 (MQ=255) cGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGAATGATGATATTCTGCCCAGCGGGGTTACCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTATGGTTGCCGCCTTTGGTTTCGATc > 1:22324/1‑150 (MQ=255) | CTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATC > NC_002947/1555281‑1555574 |
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: ‑ |