| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,870,831 | A→G | 8.7% | intergenic (+74/‑106) | 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,831 | 0 | A | G | 8.7% | 116.7 / 3.5 | 46 | intergenic (+74/‑106) | rnr/rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (25/17); new base G (2/2); total (27/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01 | |||||||||||
ATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTA > minE/2870761‑2870901 | aTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa > 1:1418681/1‑71 (MQ=255) tACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa > 1:1289810/1‑70 (MQ=255) aCCCTCTTGAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg > 1:411009/1‑71 (MQ=255) cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg < 1:3878218/70‑1 (MQ=255) cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg < 1:2692973/70‑1 (MQ=255) cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc < 1:2308556/71‑1 (MQ=255) cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc < 1:2924083/71‑1 (MQ=255) tctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt > 1:1742040/1‑70 (MQ=255) ctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa > 1:1639825/1‑71 (MQ=255) aaaaGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa > 1:1852901/1‑60 (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) cAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa > 1:1911741/1‑36 (MQ=40) aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc > 1:2913063/1‑38 (MQ=40) 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:2739239/1‑71 (MQ=255) tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg > 1:1513446/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) | ATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTA > minE/2870761‑2870901 |
| 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: ‑ |