Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,870,843 G→A 21.5% intergenic (+86/‑94) 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,8430GA21.5% 177.8 / 40.7 84intergenic (+86/‑94)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (28/37);  new base A (10/8);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGG  >  minE/2870775‑2870911
                                                                    |                                                                    
ggcctgataagacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGAcc                                                                    <  1:3522756‑M2/29‑1 (MQ=255)
 gaagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                     <  1:3105312/69‑1 (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)
          tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa                                                          >  1:3923109/1‑71 (MQ=255)
          gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                                           <  1:3374261‑M2/38‑1 (MQ=255)
          gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                                           >  1:2399343‑M2/33‑70 (MQ=255)
           acgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCAc                                                         >  1:675455‑M2/32‑71 (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:1651814‑M2/25‑69 (MQ=255)
                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaat                                                  >  1:3757004‑M2/25‑69 (MQ=255)
                  aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaat                                                  >  1:529918‑M2/25‑69 (MQ=255)
                  aagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACgctaaat                                                  >  1:1748216‑M2/25‑69 (MQ=255)
                   gACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAtt                                                  >  1:3104871/1‑70 (MQ=255)
                    gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaa                                                   <  1:2349645‑M2/46‑2 (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)
                      aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATa                                                              >  1:2457055/1‑55 (MQ=40)
                       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: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: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:789907/71‑1 (MQ=255)
                                       tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg                             <  1:1598197/71‑1 (MQ=255)
                                        gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg                             >  1:1312621/1‑70 (MQ=255)
                                         ccGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                      >  1:3671550/1‑60 (MQ=255)
                                         ccGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                      >  1:1453972/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:3081884/1‑71 (MQ=255)
                                            gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc                        <  1:2253416/71‑1 (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)
                                                       aaCGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                               >  1:3802744/1‑37 (MQ=255)
                                                       aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAatata                                           >  1:686931/1‑41 (MQ=255)
                                                       aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAata                                             <  1:3901085/39‑1 (MQ=255)
                                                       aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAata                                             <  1:193324/39‑1 (MQ=255)
                                                       aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg              <  1:410316/70‑1 (MQ=255)
                                                       aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg              >  1:351361/1‑70 (MQ=255)
                                                        aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg              >  1:1421789/1‑69 (MQ=255)
                                                               aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt                              >  1:470930/1‑46 (MQ=255)
                                                               aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt      >  1:351908/1‑70 (MQ=255)
                                                               aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt      >  1:3766914/1‑70 (MQ=255)
                                                                tCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt      >  1:2383535/1‑69 (MQ=255)
                                                                  cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTa                           >  1:54003/1‑46 (MQ=255)
                                                                  cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTa                           >  1:316145/1‑46 (MQ=255)
                                                                  cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAgg  >  1:118506/1‑71 (MQ=255)
                                                                  cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAgg  <  1:2327750/71‑1 (MQ=255)
                                                                    gCCTACATAATCACGCAAATTCAATATATTGCagag                                   <  1:1811192/36‑1 (MQ=255)
                                                                    |                                                                    
AGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGG  >  minE/2870775‑2870911

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: