| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | Exported | 2,870,833 | G→A | 14.2% | intergenic (+76/‑104) | 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 | ||
| * | Exported | 2,870,833 | 0 | G | A | 14.2% | 157.0 / 11.8 | 56 | intergenic (+76/‑104) | rnr/rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (30/18); new base A (2/6); total (32/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.27e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGT > Exported/2870763‑2870903 | aCCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg > 1:1050702/1‑71 (MQ=255) cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc < 1:4319416/71‑1 (MQ=255) cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc < 1:4787815/71‑1 (MQ=255) ccTCTTTAACAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTGCATGCCGGATGCGGCGTGAACGcc < 1:4653919/71‑1 (MQ=255) tttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt > 1:4168040/1‑71 (MQ=255) tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa < 1:3118162/68‑1 (MQ=255) tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa < 1:5477257/68‑1 (MQ=255) tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt < 1:1328406/69‑1 (MQ=255) 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:4549262/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:399353/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:5389853/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:1966575/1‑70 (MQ=255) tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa > 1:3208322/1‑71 (MQ=255) gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaa < 1:2078551/70‑1 (MQ=255) tGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc > 1:1865884/1‑49 (MQ=255) aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt > 1:1450940/1‑43 (MQ=40) aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAtt < 1:4532247/69‑1 (MQ=255) aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc < 1:528715/70‑1 (MQ=255) cAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa < 1:737980/71‑1 (MQ=255) aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt > 1:2504180/1‑40 (MQ=40) aTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt < 1:3676738/40‑1 (MQ=33) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat < 1:4211964/71‑1 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat < 1:1370511/71‑1 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat > 1:4884821/1‑71 (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: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:2936130/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: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) | ACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGT > Exported/2870763‑2870903 |
| 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: ‑ |