Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,870,829 G→A 23.4% intergenic (+72/‑108) 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
*minE2,870,8290GA23.4% 185.9 / 54.7 94intergenic (+72/‑108)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (30/42);  new base A (12/10);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

Only 100 of 103 total aligned reads displayed.

AATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAA  >  minE/2870760‑2870896
                                                                     |                                                                   
aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:2423053/70‑1 (MQ=255)
aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:566109/70‑1 (MQ=255)
 aTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   <  1:4114383/71‑1 (MQ=255)
  tACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                  >  1:1041946/1‑71 (MQ=255)
  tACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAAc                                                                  >  1:3206774/1‑71 (MQ=255)
     ccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                               <  1:1296359/71‑1 (MQ=255)
      ctctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                               >  1:708021/1‑70 (MQ=255)
         tttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:3900017/61‑1 (MQ=255)
         tttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                             <  1:739766/69‑1 (MQ=255)
          ttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                          >  1:910373/1‑71 (MQ=255)
          ttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                          >  1:1140225/1‑71 (MQ=255)
          gcgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTAt                                                           <  1:1943384‑M2/23‑1 (MQ=255)
          gcgtaggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATc                                                          >  1:3671455‑M2/48‑71 (MQ=255)
           tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:3770162/59‑1 (MQ=255)
           tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                            <  1:1184460/68‑1 (MQ=255)
           tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                         >  1:464967/1‑71 (MQ=255)
           tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATcc                                                         >  1:454735/1‑71 (MQ=255)
           cgtaggcctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTaaa                                                                   >  1:654835‑M2/47‑61 (MQ=255)
            gtaggcctgataagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCg                                                        <  1:1276829‑M2/26‑1 (MQ=255)
              aggcctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTaaa                                                                   >  1:2150327‑M2/44‑58 (MQ=255)
              aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   >  1:446532/1‑58 (MQ=255)
               ggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGAcc                                                     <  1:3522756‑M2/29‑1 (MQ=255)
                gaagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:3589001/54‑1 (MQ=255)
                gaagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                             >  1:2439585/1‑62 (MQ=255)
                gaagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                      <  1:3105312/69‑1 (MQ=255)
                 aagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                        >  1:84484/1‑66 (MQ=255)
                  ctgataagacgcggcaagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTa                                                   >  1:2123991‑M2/40‑70 (MQ=255)
                     ataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTAcaa                                                <  1:3291702‑M2/34‑2 (MQ=255)
                      ggTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATaa                                              <  1:888289/71‑1 (MQ=255)
                       gTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt                                             <  1:3427475/71‑1 (MQ=255)
                        agacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAAt                                             >  1:1847499‑M2/34‑70 (MQ=255)
                         tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:3023671/45‑1 (MQ=255)
                         tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:2285132/45‑1 (MQ=255)
                         tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa                                           >  1:3923109/1‑71 (MQ=255)
                         gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                            >  1:2399343‑M2/33‑70 (MQ=255)
                         gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                            <  1:3374261‑M2/38‑1 (MQ=255)
                          acgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCAc                                          >  1:675455‑M2/32‑71 (MQ=255)
                          aGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   >  1:3544240/1‑46 (MQ=255)
                            aTTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                            <  1:3381010/51‑1 (MQ=255)
                              tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   >  1:3700094/1‑42 (MQ=255)
                              tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                <  1:869695/45‑1 (MQ=255)
                              tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                       >  1:2262234/1‑70 (MQ=255)
                              tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                       >  1:1983495/1‑70 (MQ=255)
                              gttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGc                                        >  1:2794429‑M2/28‑69 (MQ=255)
                                cTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaaa                                     <  1:2552033/70‑1 (MQ=255)
                                 aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaat                                   >  1:529918‑M2/25‑69 (MQ=255)
                                 aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaat                                   >  1:1651814‑M2/25‑69 (MQ=255)
                                 aagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACgctaaat                                   >  1:1748216‑M2/25‑69 (MQ=255)
                                  gACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAtt                                   >  1:3104871/1‑70 (MQ=255)
                                   gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaa                                    <  1:3022600‑M2/46‑2 (MQ=255)
                                   gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaa                                    <  1:2349645‑M2/46‑2 (MQ=255)
                                   aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAtt                                   <  1:2459684/69‑1 (MQ=255)
                                   aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                  <  1:780167/70‑1 (MQ=255)
                                   aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCa                                 <  1:386602/71‑1 (MQ=255)
                                    cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                            <  1:923694‑M2/38‑1 (MQ=255)
                                     gtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATccgac                                                      >  1:405859‑M2/21‑46 (MQ=255)
                                     aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCGTTATCCGGCCTACATaa                                              >  1:2977763/1‑56 (MQ=255)
                                     aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                            <  1:1856755/42‑1 (MQ=38)
                                     aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATa                                               >  1:2457055/1‑55 (MQ=40)
                                     aaaaTGCGCTTTGTTCATGCCGGAAGCGGCGTGAACGCCTTa                                                            <  1:1851950/42‑1 (MQ=255)
                                      tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc                              <  1:429367‑M2/52‑4 (MQ=255)
                                      tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc                              <  1:688005‑M2/52‑4 (MQ=255)
                                      aaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                              <  1:3025956/71‑1 (MQ=255)
                                       aaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa                                           >  1:3154997/1‑57 (MQ=255)
                                       aaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                  <  1:918536/66‑1 (MQ=255)
                                         tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                              <  1:3970162/68‑1 (MQ=255)
                                             cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:1643017/71‑1 (MQ=255)
                                             cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:3756211/71‑1 (MQ=255)
                                             cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:714569/71‑1 (MQ=255)
                                             cTTTGTTCATGCCGGATGCGGCATGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg                         >  1:394698/1‑69 (MQ=255)
                                              tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                              <  1:3280597/63‑1 (MQ=255)
                                              tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:3334544/70‑1 (MQ=255)
                                              tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:3373258/70‑1 (MQ=255)
                                               ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat                           >  1:1006746/1‑65 (MQ=255)
                                                tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       >  1:3084404/1‑68 (MQ=255)
                                                tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       >  1:138408/1‑68 (MQ=255)
                                                 gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATAtt                          <  1:2707301/64‑1 (MQ=255)
                                                  ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       >  1:3830489/1‑66 (MQ=255)
                                                  ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAgaga                   <  1:371141/70‑1 (MQ=255)
                                                  ctCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAt                  >  1:2779281/2‑70 (MQ=2)
                                                      tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg              <  1:789907/71‑1 (MQ=255)
                                                      tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg              <  1:1598197/71‑1 (MQ=255)
                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg              >  1:1312621/1‑70 (MQ=255)
                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       >  1:1453972/1‑60 (MQ=255)
                                                        ccGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       >  1:3671550/1‑60 (MQ=255)
                                                         cGGATGCGGCGTGAACGCCTGATCCGGCCTACATAATCACGCAAATTCaa                                <  1:3465948/50‑1 (MQ=255)
                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCa                                 <  1:61650/47‑1 (MQ=255)
                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc         >  1:3537357/1‑71 (MQ=255)
                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc         <  1:2253416/71‑1 (MQ=255)
                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc         >  1:3081884/1‑71 (MQ=255)
                                                             tGCCGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                              <  1:775932/48‑1 (MQ=255)
                                                                ggcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatat                             <  1:333423/46‑1 (MQ=255)
                                                                 gcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                <  1:1412442/42‑1 (MQ=255)
                                                                  cgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaa  <  1:3122007/71‑1 (MQ=255)
                                                                  cgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaa  <  1:920698/71‑1 (MQ=255)
                                                                  cgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGCCCGGATaa  <  1:3544543/71‑1 (MQ=255)
                                                                    tGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                  <  1:2420706/37‑1 (MQ=255)
                                                                    tGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatat                             <  1:459681/42‑1 (MQ=255)
                                                                    tGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:343781/48‑1 (MQ=255)
                                                                    tGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAgaga                   <  1:2491526/52‑1 (MQ=255)
                                                                     |                                                                   
AATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAA  >  minE/2870760‑2870896

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: