Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,870,831 A→G 15.1% 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
*Exported2,870,8310AG15.1% 126.9 / 13.4 53intergenic (+74/‑106)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/17);  new base G (2/6);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

ACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACG  >  Exported/2870763‑2870899
                                                                    |                                                                    
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)
      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)
           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: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)
                         aTTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                      >  1:1586626/1‑44 (MQ=255)
                         aTTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                      >  1:3604841/1‑44 (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)
                                  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)
                                               ctCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                     <  1:238707/69‑1 (MQ=2)
                                               ctCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                    <  1:769773/70‑1 (MQ=2)
                                                  aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg                                            >  1:4954284/1‑45 (MQ=255)
                                                  aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt                  >  1:3231049/1‑71 (MQ=255)
                                                  aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                 <  1:3836467/71‑1 (MQ=2)
                                                  aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                 <  1:3973050/71‑1 (MQ=2)
                                                   tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                 <  1:4806689/70‑1 (MQ=2)
                                                   tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                 <  1:705961/70‑1 (MQ=2)
                                                   tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                 >  1:3543351/1‑70 (MQ=2)
                                                    gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatata                               >  1:1040324/1‑56 (MQ=255)
                                                    gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt                  <  1:66238/69‑1 (MQ=255)
                                                    gCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                 >  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)
                                                                 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)
                                                                   aaCGCCTTATCCGGCCTACATAATCACGCAAATTAATATATTGCAGAGATCATGTa               >  1:2381709/1‑56 (MQ=255)
                                                                    |                                                                    
ACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACG  >  Exported/2870763‑2870899

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: