Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,153 | G→T | 15.6% | R113S (CGT→AGT) | rsmG ← | 16S RNA methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,153 | 0 | G | T | 15.6% | 95.2 / 3.2 | 31 | R113S (CGT→AGT) | rsmG | 16S RNA methyltransferase |
Reads supporting (aligned to +/- strand): ref base G (18/8); new base T (2/3); total (20/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGTTCTCCATGCTGCTGAAGGCGCGGGAGATGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAACATCCAGCCAGTTATCACGATCGTTGTGGATAAAAGA > NC_002947/2021‑2282 | aGTTCTCCATGCTGCTGAAGGCGCGGGAGATGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACgggg < 2:357210/136‑2 (MQ=255) ttCTCCATGCTGCTGAAGGCGCGGGAGATGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGtct < 2:468081/136‑1 (MQ=255) gctgctGAAGGCGCGGGAGATGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCa > 2:384153/1‑136 (MQ=255) aaGGCGCGGGAGATGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCTGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGt > 1:419110/1‑136 (MQ=255) aaGGCGCGGGAGATGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGt > 1:129667/1‑136 (MQ=255) agaTGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCTGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACTGGTCTTCTTGCCATTGGCGTCCAGCACGGt > 1:122131/1‑136 (MQ=255) agaTGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCTGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACTGGTCTTCTTGCCATTGGCGTCCAGCACGGt > 1:455372/1‑136 (MQ=255) agaTGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGt > 2:464238/1‑136 (MQ=255) agaTGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGt > 2:461907/1‑136 (MQ=255) gTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGt < 2:334550/99‑1 (MQ=255) gTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGt > 1:334550/1‑99 (MQ=255) cGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTtgtg > 1:126799/1‑136 (MQ=255) aaCGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCa < 2:289739/115‑1 (MQ=255) aaCGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCa > 1:289739/1‑115 (MQ=255) aaCGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTtgtgtg < 1:137381/136‑1 (MQ=255) gtgtgCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCTTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACTGGTCTTCTTGACATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGc < 1:63293/135‑1 (MQ=255) ccGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGAt > 1:428214/1‑136 (MQ=255) aaGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCg > 1:237892/1‑136 (MQ=255) gACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCg > 2:462034/1‑136 (MQ=255) gtgGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCa > 1:457030/1‑136 (MQ=255) gATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCAccg < 2:319564/136‑1 (MQ=255) aaCCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACATGAGTCAGGAAACGGGTCTTCTTGACATTGGCGTCCAGCACGGTGACGCGCTTGGGTGGATGCAGGATAGCCAGCTGGATGCCGGGCATGCCAccgccc > 1:310586/1‑135 (MQ=255) gTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGc > 2:443471/1‑136 (MQ=255) ggTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAAc > 1:379626/1‑136 (MQ=255) ggTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAAc > 1:205058/1‑136 (MQ=255) ggTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAAc > 2:400252/1‑136 (MQ=255) cAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAACATccagc > 1:134385/1‑136 (MQ=255) cTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAACATCCAGCCAGTTATCACGATCGTTGTgg < 2:457030/136‑1 (MQ=255) aGGAAACTGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAACATCCAGCCAGTTATCACGATCGTTGTGGATAAAAg < 2:455372/136‑1 (MQ=255) aGGAAACTGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAACATCCAGCCAGTTATCACGATCGTTGTGGATAAAAg < 2:122131/136‑1 (MQ=255) ggAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAACATCCAGCCAGTTATCACGATCGTTGTGGATAAAAGa < 1:464238/136‑1 (MQ=255) ggAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAACATCCAGCCAGTTATCACGATCGTTGTGGATAAAAGa < 1:461907/136‑1 (MQ=255) | AGTTCTCCATGCTGCTGAAGGCGCGGGAGATGATTCCGTCGAACGGTTGTGCCGGCTGGAAGGCCTCGACCCGGCTGTGGATAACCGTGAGGTTGTCCAGCTTGAGTTCCATTTTCACCTGAGTCAGGAAACGGGTCTTCTTGCCATTGGCGTCCAGCACGGTGACGCGCTTGTGTGGATGCAGGATAGCCAGCGGGATGCCGGGCATGCCACCGCCGCTGCCAACATCCAGCCAGTTATCACGATCGTTGTGGATAAAAGA > NC_002947/2021‑2282 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |