Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 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
*Exported2,870,8650AT13.6% 208.6 / 22.7 81intergenic (+108/‑72)rnr/rlmBexoribonuclease 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  >  Exported/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)
              ctCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                                                         <  1:238707/69‑1 (MQ=2)
              ctCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                        <  1:769773/70‑1 (MQ=2)
                 aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt                                                      >  1:3231049/1‑71 (MQ=255)
                 aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                                     <  1:3973050/71‑1 (MQ=2)
                 aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                                     <  1:3836467/71‑1 (MQ=2)
                  tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                                     >  1:3543351/1‑70 (MQ=2)
                  tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                                     <  1:705961/70‑1 (MQ=2)
                  tGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                                     <  1:4806689/70‑1 (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: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)
                                  aaCGCCTTATCCGGCCTACATAATCACGCAAATTAATATATTGCAGAGATCATGTa                                                   >  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  >  Exported/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: