| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,870,865 | A→T | 13.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 | 13.6% | 208.6 / 22.7 | 81 | intergenic (+108/‑72) | rnr/rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (51/19); new base T (2/9); total (53/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.44e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01 | |||||||||||
CAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATT > minE/2870796‑2870935 | cAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa < 1:737980/71‑1 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat < 1:1370511/71‑1 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat > 1:4884821/1‑71 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat < 1:4211964/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) aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt > 1:3231049/1‑71 (MQ=255) aTGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg < 1:3973050/71‑1 (MQ=2) aTGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg < 1:3836467/71‑1 (MQ=2) tGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg > 1:3543351/1‑70 (MQ=2) tGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg < 1:705961/70‑1 (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:5315801/1‑58 (MQ=255) aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga > 1:5062382/1‑58 (MQ=255) aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga > 1:2936130/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:5474362/1‑38 (MQ=255) aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAat > 1:2760744/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:3675356/1‑70 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt > 1:3119327/1‑70 (MQ=255) aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt > 1:645033/1‑70 (MQ=255) ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg > 1:285360/1‑71 (MQ=255) ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg > 1:2819967/1‑71 (MQ=255) cGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:2983924/1‑35 (MQ=255) cGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:3558496/1‑35 (MQ=255) cGGCCTACATAATCACGCAAATTCAATATATTGCa > 1:3672851/1‑35 (MQ=255) cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAgg > 1:4230332/1‑41 (MQ=255) cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc > 1:485966/1‑43 (MQ=255) cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg > 1:648384/1‑55 (MQ=255) cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg > 1:5458084/1‑55 (MQ=255) cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCaaaa > 1:798657/1‑69 (MQ=255) cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAAc < 1:5144743/70‑1 (MQ=255) cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACg < 1:3719669/71‑1 (MQ=255) cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACg < 1:3272388/71‑1 (MQ=255) aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg > 1:376366/1‑51 (MQ=255) aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg > 1:4076443/1‑51 (MQ=255) aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg > 1:2996772/1‑51 (MQ=255) aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg > 1:2442516/1‑51 (MQ=255) aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGt < 1:3400942/71‑1 (MQ=255) aaaTCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGa < 1:5558718/43‑1 (MQ=11) aaaTCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGa < 1:69229/43‑1 (MQ=11) aaaTCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGa < 1:4626623/43‑1 (MQ=11) aaTCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTa > 1:1967699/1‑70 (MQ=255) aTCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATa > 1:1810995/1‑44 (MQ=255) tCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTa > 1:2501356/1‑49 (MQ=255) tCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTaaa > 1:1408647/1‑70 (MQ=255) tCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaac < 1:2019605/71‑1 (MQ=255) cACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaac < 1:585002/70‑1 (MQ=255) aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAAc > 1:1370270/1‑44 (MQ=255) aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAAc > 1:3272032/1‑44 (MQ=255) aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTa < 1:4401021/66‑1 (MQ=255) aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTa < 1:3929890/66‑1 (MQ=255) aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTa < 1:5433957/66‑1 (MQ=255) cGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaac < 1:2686869/71‑1 (MQ=255) gCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaac > 1:2375848/1‑70 (MQ=255) gCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACg < 1:2995769/71‑1 (MQ=255) aaaTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAg < 1:5019282/71‑1 (MQ=255) aaTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGc < 1:2351983/46‑1 (MQ=255) aaTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGc < 1:3103855/46‑1 (MQ=255) aTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg > 1:596739/1‑42 (MQ=255) aTTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTa > 1:1227187/1‑71 (MQ=255) ttCAATATATTGCAGAGATCATGTAGGCCGTATAAACGt > 1:2453387/1‑39 (MQ=255) ttCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt > 1:3029351/1‑39 (MQ=255) ttCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt > 1:357125/1‑39 (MQ=255) ttCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt > 1:2204111/1‑39 (MQ=255) tCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACa > 1:4021441/1‑71 (MQ=255) tCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACa > 1:4461709/1‑71 (MQ=255) cAATATATTGCAGAGATCATGTAGGCCGTATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACAt > 1:3206258/1‑71 (MQ=255) cAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg > 1:631960/1‑47 (MQ=255) aaTATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACAtt < 1:2988855/71‑1 (MQ=255) | CAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATT > minE/2870796‑2870935 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |