| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,870,836 | T→A | 9.1% | intergenic (+79/‑101) | 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,836 | 0 | T | A | 9.1% | 117.2 / 4.4 | 44 | intergenic (+79/‑101) | rnr/rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base T (25/15); new base A (2/2); total (27/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.34e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01 | |||||||||||
TCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAG > minE/2870767‑2870902 | tctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt > 1:1742040/1‑70 (MQ=255) ctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa > 1:1639825/1‑71 (MQ=255) aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc < 1:2062498/71‑1 (MQ=255) agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc > 1:3865793/1‑67 (MQ=255) agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa > 1:3267578/1‑70 (MQ=255) 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:1742993/1‑71 (MQ=255) cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa > 1:2753095/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) | TCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAG > minE/2870767‑2870902 |
| 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: ‑ |