Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,870,848 C→G 14.8% intergenic (+91/‑89) 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,8480CG14.8% 155.1 / 14.3 54intergenic (+91/‑89)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/10);  new base G (2/6);  total (38/16)
Fisher's exact test for biased strand distribution p-value = 5.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACG  >  Exported/2870778‑2870918
                                                                      |                                                                      
agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                        <  1:2947504/71‑1 (MQ=255)
  agGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                       >  1:2556739/1‑70 (MQ=255)
  agGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAt                                                                      <  1:4607426/71‑1 (MQ=255)
      ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt                                                                   >  1:5389853/1‑70 (MQ=255)
      ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt                                                                   >  1:1966575/1‑70 (MQ=255)
      ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt                                                                   >  1:399353/1‑70 (MQ=255)
      ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt                                                                   >  1:4549262/1‑70 (MQ=255)
       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: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)
                                   aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg                                                               >  1:4954284/1‑45 (MQ=255)
                                   aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt                                     >  1:3231049/1‑71 (MQ=255)
                                   aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                    <  1:3973050/71‑1 (MQ=2)
                                   aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                    <  1:3836467/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:2936130/1‑58 (MQ=255)
                                          aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga                                           >  1:5062382/1‑58 (MQ=255)
                                          aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCaga                                           >  1:5315801/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:3344522/1‑70 (MQ=255)
                                                            aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt             >  1:3119327/1‑70 (MQ=255)
                                                            aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt             >  1:3675356/1‑70 (MQ=255)
                                                              ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg          >  1:2819967/1‑71 (MQ=255)
                                                              ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg          >  1:285360/1‑71 (MQ=255)
                                                               cGGCCTACATAATCACGCAAATTCAATATATTGCa                                             >  1:2983924/1‑35 (MQ=255)
                                                               cGGCCTACATAATCACGCAAATTCAATATATTGCa                                             >  1:3672851/1‑35 (MQ=255)
                                                               cGGCCTACATAATCACGCAAATTCAATATATTGCa                                             >  1:3558496/1‑35 (MQ=255)
                                                                      cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAgg                                >  1:4230332/1‑41 (MQ=255)
                                                                      cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc                              >  1:485966/1‑43 (MQ=255)
                                                                      cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg                  >  1:5458084/1‑55 (MQ=255)
                                                                      cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg                  >  1:648384/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)
                                                                      |                                                                      
AGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACG  >  Exported/2870778‑2870918

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: