Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 558,477 | A→C | 6.9% | E21D (GAA→GAC) | rplR → | 50S ribosomal protein L18 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 558,477 | 0 | A | C | 6.9% | 177.3 / 4.5 | 58 | E21D (GAA→GAC) | rplR | 50S ribosomal protein L18 |
Reads supporting (aligned to +/- strand): ref base A (28/26); new base C (2/2); total (30/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01 |
CAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGC > NC_002947/558331‑558623 | cAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACt > 5:64890/1‑149 (MQ=255) gCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGaa < 3:191513/149‑1 (MQ=255) cTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGt > 6:157823/1‑149 (MQ=255) ttACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGagttc > 1:310794/1‑145 (MQ=255) ttACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAg < 8:287363/147‑1 (MQ=255) ttACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGta > 7:287363/1‑148 (MQ=255) cAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGt < 3:248738/148‑1 (MQ=255) aaggtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgca < 2:310794/149‑2 (MQ=255) aggtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgcgc < 6:64890/149‑1 (MQ=255) aggtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgcgc < 4:139864/149‑1 (MQ=255) aaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGACGCACGAACTCGAAGTCGTGCGCCtgt > 3:257748/1‑149 (MQ=255) aaacgtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGGCTCAAGATGCACGAACTCGAAGTCGTGCGCCtgt > 3:257752/1‑149 (MQ=255) aggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGc < 4:261057/149‑1 (MQ=255) cgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcctc < 6:333038/149‑1 (MQ=255) gAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcag > 5:127863/1‑148 (MQ=255) tcgtAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcatc < 6:89457/114‑4 (MQ=255) tAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGc < 2:320776/116‑1 (MQ=255) tAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGc > 1:320776/1‑116 (MQ=255) tAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTc > 5:69041/1‑122 (MQ=255) tAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTc < 6:69041/122‑1 (MQ=255) agaagCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGt > 1:50353/1‑148 (MQ=255) acgaagaagTCGGGCCTAGTACATGCCCGACAAAAAAGTTATGCGACTGCGGCGCGCTCGCAACGCACGTCGCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg < 4:257752/147‑1 (MQ=255) aagaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGACGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg < 4:257748/149‑1 (MQ=255) tAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcagca < 8:333410/117‑1 (MQ=255) tAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcagca > 7:333410/1‑117 (MQ=255) aGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgcag > 7:137167/1‑149 (MQ=255) aGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAg < 5:166938/122‑1 (MQ=255) aGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAg > 6:166938/1‑122 (MQ=255) gTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAATGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCt > 1:244842/1‑113 (MQ=255) gTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCt < 2:244842/113‑1 (MQ=255) aaaTGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCt < 2:50353/148‑1 (MQ=255) aaTGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg < 6:127863/148‑1 (MQ=255) gACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCATGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGc > 5:28446/1‑147 (MQ=255) aaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcct < 5:184425/85‑1 (MQ=255) aaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcct > 6:184425/1‑85 (MQ=255) tATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTg < 8:137167/148‑1 (MQ=255) aTTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCAttt > 6:319544/1‑106 (MQ=255) aTTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCAttt < 5:319544/106‑1 (MQ=255) ttCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGAc < 5:157823/149‑1 (MQ=255) gACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCa > 5:237466/1‑99 (MQ=255) gACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCa < 6:237466/99‑1 (MQ=255) gCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGAcc > 6:289183/1‑149 (MQ=255) cgcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGt > 2:305228/1‑148 (MQ=255) tCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGAc < 3:204440/133‑1 (MQ=255) tCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGAc > 4:204440/1‑133 (MQ=255) gCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCAc > 8:93785/1‑149 (MQ=255) cACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCATGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGCCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACt < 6:28446/149‑1 (MQ=255) aCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTg > 4:22428/1‑149 (MQ=255) aCGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTg > 4:187099/1‑149 (MQ=255) tctcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGc > 4:86443/1‑148 (MQ=255) tctcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCa < 3:86443/149‑1 (MQ=255) tcAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCaa > 5:56679/1‑148 (MQ=255) ggcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCg > 4:290369/4‑149 (MQ=255) aTGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGa < 3:22428/149‑1 (MQ=255) gCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTg > 2:105942/1‑121 (MQ=255) gCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTg < 1:105942/121‑1 (MQ=255) aCGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGTTGCCAGTGGCGCCATcacgcag > 5:154222/1‑143 (MQ=255) gAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCGCGACCGTGGGCAAAGACCTGCGTGATGGCGCCACTGGCAACCTCGACGCGGc > 4:331895/1‑149 (MQ=255) | CAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGC > NC_002947/558331‑558623 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |