| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,870,843 | G→T | 14.2% | intergenic (+86/‑94) | 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,843 | 0 | G | T | 14.2% | 161.9 / 9.8 | 56 | intergenic (+86/‑94) | rnr/rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (37/11); new base T (2/6); total (39/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.01e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01 | |||||||||||
AAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAG > minE/2870774‑2870910 | aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc > 1:164125/1‑71 (MQ=255) agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc < 1:2947504/71‑1 (MQ=255) agGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa > 1:2556739/1‑70 (MQ=255) agGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAt < 1:4607426/71‑1 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:5389853/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:1966575/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:399353/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:4549262/1‑70 (MQ=255) tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa > 1:3208322/1‑71 (MQ=255) gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaa < 1:2078551/70‑1 (MQ=255) aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAtt < 1:4532247/69‑1 (MQ=255) aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc < 1:528715/70‑1 (MQ=255) cAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa < 1:737980/71‑1 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat > 1:4884821/1‑71 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat < 1:4211964/71‑1 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat < 1:1370511/71‑1 (MQ=255) tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:1168042/1‑68 (MQ=255) ctCATGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGAt < 1:238707/69‑1 (MQ=2) ctCATGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATc < 1:769773/70‑1 (MQ=2) aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg > 1:4954284/1‑45 (MQ=255) aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt > 1:3231049/1‑71 (MQ=255) aTGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg < 1:3836467/71‑1 (MQ=2) aTGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg < 1:3973050/71‑1 (MQ=2) tGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg < 1:705961/70‑1 (MQ=2) tGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg > 1:3543351/1‑70 (MQ=2) tGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg < 1:4806689/70‑1 (MQ=2) gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatata > 1:1040324/1‑56 (MQ=255) gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt < 1:66238/69‑1 (MQ=255) gCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg > 1:4622790/1‑69 (MQ=2) gCCAGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatata > 1:215123/1‑56 (MQ=255) aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga > 1:2936130/1‑58 (MQ=255) aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga > 1:5315801/1‑58 (MQ=255) aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga > 1:5062382/1‑58 (MQ=255) aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc < 1:2857564/71‑1 (MQ=255) tGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa > 1:3216614/1‑46 (MQ=255) tGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt > 1:185567/1‑60 (MQ=255) gcggcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa > 1:1504040/1‑45 (MQ=255) tGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATa > 1:4364611/1‑68 (MQ=255) aaCGCCTTATCCGGCCTACATAATCACGCAAATTCaa > 1:2510881/1‑37 (MQ=255) aaCGCCTTATCCGGCCTACATAATCACGCAAATTCaa > 1:2024884/1‑37 (MQ=255) aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAat > 1:2760744/1‑38 (MQ=255) aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAat > 1:5474362/1‑38 (MQ=255) aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTa > 1:454123/1‑57 (MQ=255) aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg < 1:2576102/70‑1 (MQ=255) aaCGCCTTATCCGGCCTACATAATCACGCAAATT‑AATATATTGCAGAGATCATGTa > 1:2381709/1‑56 (MQ=255) gCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGt > 1:5519933/1‑71 (MQ=255) ttATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTa > 1:4179547/1‑51 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt > 1:3344522/1‑70 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt > 1:3119327/1‑70 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt > 1:645033/1‑70 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt > 1:3675356/1‑70 (MQ=255) ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg > 1:285360/1‑71 (MQ=255) ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg > 1:2819967/1‑71 (MQ=255) cGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:3558496/1‑35 (MQ=255) cGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:3672851/1‑35 (MQ=255) cGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:2983924/1‑35 (MQ=255) | AAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAG > minE/2870774‑2870910 |
| 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: ‑ |