Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,870,827 G→A 20.0% intergenic (+70/‑110) 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,8270GA20.0% 218.9 / 43.8 95intergenic (+70/‑110)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (35/41);  new base A (10/9);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 104 total aligned reads displayed.

CAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAA  >  minE/2870759‑2870896
                                                                    |                                                                     
cAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGt                                                                      >  1:1366039/1‑70 (MQ=255)
 aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATgcggcg                                                                       <  1:344591/68‑1 (MQ=255)
 aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATgcggcg                                                                       <  1:3152978/68‑1 (MQ=255)
 aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:566109/70‑1 (MQ=255)
 aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:2423053/70‑1 (MQ=255)
 aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGa                                                                    <  1:973458/71‑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:1140225/1‑71 (MQ=255)
           ttAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                          >  1:910373/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)
                        aagacgcggcaagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAt                                             >  1:2767252‑M2/35‑71 (MQ=255)
                         agacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAAt                                             >  1:1847499‑M2/34‑70 (MQ=255)
                          tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:2285132/45‑1 (MQ=255)
                          tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:3023671/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)
                               tGCTGACAAAATGCGCTTTGTTCATGCCGGATgcggcg                                                                       >  1:3598621/1‑38 (MQ=255)
                               tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   >  1:3700094/1‑42 (MQ=255)
                               tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                <  1:869695/45‑1 (MQ=255)
                               tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                       >  1:1983495/1‑70 (MQ=255)
                               tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                       >  1:2262234/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)
                                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaat                                   >  1:3757004‑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)
                                    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)
                                      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:688005‑M2/52‑4 (MQ=255)
                                       tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc                              <  1:429367‑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:3756211/71‑1 (MQ=255)
                                              cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:714569/71‑1 (MQ=255)
                                              cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:1643017/71‑1 (MQ=255)
                                              cTTTGTTCATGCCGGATGCGGCATGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg                         >  1:394698/1‑69 (MQ=255)
                                               tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                              <  1:3280597/63‑1 (MQ=255)
                                               tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       <  1:3373258/70‑1 (MQ=255)
                                                ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat                           >  1:1006746/1‑65 (MQ=255)
                                                 tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       >  1:138408/1‑68 (MQ=255)
                                                 tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                       >  1:3084404/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:1598197/71‑1 (MQ=255)
                                                       tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg              <  1:789907/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)
                                                                    |                                                                     
CAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAA  >  minE/2870759‑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: