| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 59,338 | C→T | intergenic (+59/+349) | yabQ → / ← rluA | protein YabQ/23S rRNA pseudouridine(746) and tRNA pseudouridine(32) synthase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 59,338 | 0 | C | T | 95.7% | 72.3 / ‑4.0 | 23 | intergenic (+59/+349) | yabQ/rluA | protein YabQ/23S rRNA pseudouridine(746) and tRNA pseudouridine(32) synthase |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (16/6); total (17/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01 | |||||||||||
GAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACC > NC_000913/59264‑59413 | gAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAAt > 1:14320/1‑80 (MQ=255) gAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGa > 2:14319/1‑80 (MQ=255) ccAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGa < 1:14284/81‑1 (MQ=255) ccAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGa > 1:14316/1‑81 (MQ=255) aTAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTg > 2:14327/1‑80 (MQ=255) aGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAg < 1:14285/81‑1 (MQ=255) gAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTa > 1:14314/1‑79 (MQ=255) aaGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTa > 1:14309/1‑78 (MQ=255) aGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTagag > 1:14321/1‑80 (MQ=255) ggggCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTagag > 1:14313/1‑79 (MQ=255) ggCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGagaaga < 1:14286/81‑1 (MQ=255) cTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAgg > 1:14311/1‑81 (MQ=255) cccTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACg > 2:14312/1‑81 (MQ=255) cTTCAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCa < 1:14287/81‑1 (MQ=255) ttCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGc > 2:14323/1‑79 (MQ=255) cAACCTTAGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCaa < 1:14288/79‑1 (MQ=255) aGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGcgtcg < 2:14289/79‑1 (MQ=255) aGCAGTAGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTg > 2:14322/1‑81 (MQ=255) aGCGTGGGATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGa > 1:14330/1‑80 (MQ=255) gggATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGa > 1:14325/1‑81 (MQ=255) gggATGATTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGa > 2:14326/1‑81 (MQ=255) gatTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAcc > 1:14317/1‑81 (MQ=255) atTTTACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAcc > 1:14339/1‑80 (MQ=255) | GAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACC > NC_000913/59264‑59413 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |