| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,870,848 | C→G | 8.0% | intergenic (+91/‑89) | 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,848 | 0 | C | G | 8.0% | 154.5 / 3.1 | 50 | intergenic (+91/‑89) | rnr/rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (29/17); new base G (2/2); total (31/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.29e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01 | |||||||||||
TTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGC > minE/2870784‑2870912 | ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:2663924/1‑70 (MQ=255)ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATc < 1:1877163/71‑1 (MQ=255) tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa > 1:2724357/1‑71 (MQ=255) tGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAt > 1:1385038/1‑70 (MQ=255) aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc < 1:2085475/70‑1 (MQ=255) aaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAtt < 1:1681453/66‑1 (MQ=255) aaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata > 1:3590558/1‑70 (MQ=255) gcgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc < 1:1032439/63‑1 (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) tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAGTc < 1:3044267/47‑1 (MQ=12) tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATc > 1:735710/1‑47 (MQ=255) 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) aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc < 1:1791198/52‑1 (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:2712385/1‑69 (MQ=255) aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg > 1:224174/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:29821/1‑43 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt > 1:805355/1‑43 (MQ=255) tCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa < 1:381398/71‑1 (MQ=255) tCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa < 1:3257019/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) | TTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGC > minE/2870784‑2870912 |
| 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: ‑ |