Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,870,855 A→G 15.9% intergenic (+98/‑82) 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,8550AG15.9% 159.8 / 20.8 69intergenic (+98/‑82)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (43/15);  new base G (2/9);  total (45/24)
Fisher's exact test for biased strand distribution p-value = 7.59e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

TAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAAC  >  Exported/2870785‑2870923
                                                                      |                                                                    
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:1370511/71‑1 (MQ=255)
                tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat                                                      <  1:4211964/71‑1 (MQ=255)
                tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat                                                      >  1:4884821/1‑71 (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: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:2024884/1‑37 (MQ=255)
                                             aaCGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                           >  1:2510881/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)
                                                gCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGt                      >  1:5519933/1‑71 (MQ=255)
                                                   ttATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTa                                       >  1:4179547/1‑51 (MQ=255)
                                                     aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt                  >  1:645033/1‑70 (MQ=255)
                                                     aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt                  >  1:3675356/1‑70 (MQ=255)
                                                     aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt                  >  1:3344522/1‑70 (MQ=255)
                                                     aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt                  >  1:3119327/1‑70 (MQ=255)
                                                       ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg               >  1:2819967/1‑71 (MQ=255)
                                                       ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg               >  1:285360/1‑71 (MQ=255)
                                                        cGGCCTACATAATCACGCAAATTCAATATATTGCa                                                  >  1:3672851/1‑35 (MQ=255)
                                                        cGGCCTACATAATCACGCAAATTCAATATATTGCa                                                  >  1:3558496/1‑35 (MQ=255)
                                                        cGGCCTACATAATCACGCAAATTCAATATATTGCa                                                  >  1:2983924/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:3272388/71‑1 (MQ=255)
                                                               cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACg       <  1:3719669/71‑1 (MQ=255)
                                                                aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                          >  1:2442516/1‑51 (MQ=255)
                                                                aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                          >  1:376366/1‑51 (MQ=255)
                                                                aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                          >  1:4076443/1‑51 (MQ=255)
                                                                aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                          >  1:2996772/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:3272032/1‑44 (MQ=255)
                                                                      aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAAc                           >  1:1370270/1‑44 (MQ=255)
                                                                      aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTa     <  1:4401021/66‑1 (MQ=255)
                                                                      aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTa     <  1:5433957/66‑1 (MQ=255)
                                                                      aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTa     <  1:3929890/66‑1 (MQ=255)
                                                                      |                                                                    
TAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAAC  >  Exported/2870785‑2870923

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: