| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,870,865 | A→T | 8.6% | intergenic (+108/‑72) | rnr → / → rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,870,865 | 0 | A | T | 8.6% | 154.8 / 6.6 | 58 | intergenic (+108/‑72) | rnr/rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (36/17); new base T (2/3); total (38/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.28e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01 | |||||||||||
AATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATT > minE/2870799‑2870935 | aaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata > 1:3590558/1‑70 (MQ=255) tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg > 1:1111689/1‑68 (MQ=255) tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg > 1:594597/1‑68 (MQ=255) ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:3590846/1‑69 (MQ=255) tactCATGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAgaga < 1:2149275/68‑1 (MQ=0) tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCag < 1:214200/69‑1 (MQ=255) tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:1678999/1‑68 (MQ=255) gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCagag < 1:111531/70‑1 (MQ=255) cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGGAAATTCAatata > 1:915551/1‑59 (MQ=255) cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa > 1:2753095/1‑71 (MQ=255) cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa > 1:1742993/1‑71 (MQ=255) aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt > 1:1381489/1‑71 (MQ=255) aTGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCAt < 1:1497375/70‑1 (MQ=2) aTGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg > 1:3755738/1‑71 (MQ=2) tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg > 1:1513446/1‑71 (MQ=255) tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg > 1:2739239/1‑71 (MQ=255) gCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg > 1:2371438/1‑69 (MQ=2) gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:3613271/1‑57 (MQ=255) aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc < 1:1259442/71‑1 (MQ=255) cggcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGa < 1:2639664/71‑1 (MQ=255) gcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATa < 1:2223766/71‑1 (MQ=255) cgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGAt > 1:1122297/1‑69 (MQ=255) cgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaa < 1:1004566/71‑1 (MQ=255) gTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg > 1:981004/1‑47 (MQ=255) gAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg > 1:2682595/1‑71 (MQ=255) aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg > 1:224174/1‑69 (MQ=255) aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg > 1:2712385/1‑69 (MQ=255) aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt < 1:546197/70‑1 (MQ=255) aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTa < 1:1516463/71‑1 (MQ=255) cGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg > 1:259278/1‑71 (MQ=255) tATCCGGCCTACATAATCACGCAAATTCAATATATTg > 1:2219736/1‑37 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt > 1:805355/1‑43 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt > 1:29821/1‑43 (MQ=255) tCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa < 1:3257019/71‑1 (MQ=255) tCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa < 1:381398/71‑1 (MQ=255) ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg > 1:987823/1‑71 (MQ=255) cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg < 1:1960297/70‑1 (MQ=255) gCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGc > 1:3932473/1‑70 (MQ=255) ccTACATAATCACGCAAATTCAATATATTGCAGAGATc > 1:811427/1‑38 (MQ=255) ccTACATAATCACGCAAATTCAATATATTGCAGAGATc > 1:130741/1‑38 (MQ=255) aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa > 1:2664524/1‑61 (MQ=255) tAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGt > 1:1169558/1‑70 (MQ=255) aaTCGCGCAAATTCTATATATTGCAGAGATCATGTAGGc < 1:1227258/39‑1 (MQ=14) aTCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCgg < 1:2010122/41‑1 (MQ=255) aTCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTaa < 1:3012185/70‑1 (MQ=255) aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCaaa < 1:86551/61‑1 (MQ=255) aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaa > 1:3403788/1‑71 (MQ=255) aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaa < 1:2043393/71‑1 (MQ=255) gCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaac < 1:1817129/70‑1 (MQ=255) cAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGa < 1:3847366/71‑1 (MQ=255) aaaTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGt > 1:3052319/1‑45 (MQ=255) aaaTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAg > 1:2956410/1‑71 (MQ=255) aaTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCaaaa > 1:2875937/1‑57 (MQ=255) aTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg > 1:1052188/1‑42 (MQ=255) aTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg > 1:1047684/1‑42 (MQ=255) aTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTa > 1:2197697/1‑71 (MQ=255) aTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTa > 1:2235854/1‑71 (MQ=255) ttCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTa > 1:612737/1‑70 (MQ=255) aaTATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACAtt < 1:2516928/71‑1 (MQ=255) | AATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATT > minE/2870799‑2870935 |
| 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: ‑ |