| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | Exported | 2,870,827 | G→A | 13.3% | intergenic (+70/‑110) | 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,827 | 0 | G | A | 13.3% | 169.9 / 13.1 | 60 | intergenic (+70/‑110) | rnr/rlmB | exoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (31/21); new base A (2/6); total (33/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.24e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01 | |||||||||||
TCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCC > Exported/2870758‑2870890 | tCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGCATgcggcg > 1:4703482/1‑70 (MQ=255) cAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGt > 1:1556590/1‑70 (MQ=255) aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGa < 1:3540087/71‑1 (MQ=255) aCCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg > 1:1050702/1‑71 (MQ=255) cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc < 1:4787815/71‑1 (MQ=255) cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc < 1:4319416/71‑1 (MQ=255) ccTCTTTAACAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTGCATGCCGGATGCGGCGTGAACGcc < 1:4653919/71‑1 (MQ=255) tttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGt < 1:3995475/60‑1 (MQ=255) tttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt > 1:4168040/1‑71 (MQ=255) tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa < 1:5477257/68‑1 (MQ=255) tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa < 1:3118162/68‑1 (MQ=255) tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt < 1:1328406/69‑1 (MQ=255) aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAAc > 1:1308489/1‑59 (MQ=255) aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAAc > 1:1382266/1‑59 (MQ=255) aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAAc > 1:413789/1‑59 (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:1966575/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:4549262/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:399353/1‑70 (MQ=255) ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt > 1:5389853/1‑70 (MQ=255) tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa > 1:3208322/1‑71 (MQ=255) aTTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg < 1:829195/42‑1 (MQ=255) aTTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg < 1:2625817/42‑1 (MQ=255) aTTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa > 1:3604841/1‑44 (MQ=255) aTTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa > 1:1586626/1‑44 (MQ=255) gCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaa < 1:2078551/70‑1 (MQ=255) tGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc > 1:1865884/1‑49 (MQ=255) gACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg > 1:3387439/1‑36 (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) aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGaa > 1:1443164/1‑35 (MQ=40) aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt > 1:2504180/1‑40 (MQ=40) aTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt < 1:3676738/40‑1 (MQ=33) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat > 1:4884821/1‑71 (MQ=255) tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat < 1:1370511/71‑1 (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) aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg > 1:4954284/1‑45 (MQ=255) 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:705961/70‑1 (MQ=2) tGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg < 1:4806689/70‑1 (MQ=2) tGCCGGATGCGGCATGAGTACCATATCCTTCCT‑GAAAATCGCGCAAATTCTATATATTGCAGAGATCATg > 1:3543351/1‑70 (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:5062382/1‑58 (MQ=255) aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga > 1:5315801/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) | TCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCC > Exported/2870758‑2870890 |
| 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: ‑ |