Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,870,850 T→A 7.8% intergenic (+93/‑87) 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,8500TA 7.8% 138.9 / 4.2 51intergenic (+93/‑87)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (30/17);  new base A (2/2);  total (32/19)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

TTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGT  >  minE/2870784‑2870919
                                                                  |                                                                     
ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAAt                                                                    >  1:2663924/1‑70 (MQ=255)
ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATc                                                                   <  1:1877163/71‑1 (MQ=255)
 tAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCa                                                                  >  1:2724357/1‑71 (MQ=255)
         tGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAt                                                           >  1:1385038/1‑70 (MQ=255)
           aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                                         <  1:2085475/70‑1 (MQ=255)
              aaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAtt                                                          <  1:1681453/66‑1 (MQ=255)
               aaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                                                     >  1:3590558/1‑70 (MQ=255)
                  gcgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                                         <  1:1032439/63‑1 (MQ=255)
                      tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg                                                >  1:1111689/1‑68 (MQ=255)
                      tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg                                                >  1:594597/1‑68 (MQ=255)
                       ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                              >  1:3590846/1‑69 (MQ=255)
                        tactCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAgaga                                          <  1:2149275/68‑1 (MQ=0)
                        tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAGTc                                                                   <  1:3044267/47‑1 (MQ=12)
                        tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATc                                                                   >  1:735710/1‑47 (MQ=255)
                        tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCag                                             <  1:214200/69‑1 (MQ=255)
                        tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                              >  1:1678999/1‑68 (MQ=255)
                         gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCagag                                           <  1:111531/70‑1 (MQ=255)
                            cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGGAAATTCAatata                                                   >  1:915551/1‑59 (MQ=255)
                            cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa                                       >  1:1742993/1‑71 (MQ=255)
                            cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa                                       >  1:2753095/1‑71 (MQ=255)
                             aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                                         <  1:1791198/52‑1 (MQ=255)
                             aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt                                      >  1:1381489/1‑71 (MQ=255)
                             aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCAt                                      <  1:1497375/70‑1 (MQ=2)
                             aTGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                     >  1:3755738/1‑71 (MQ=2)
                              tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg                                     >  1:1513446/1‑71 (MQ=255)
                              tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATg                                     >  1:2739239/1‑71 (MQ=255)
                               gCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATg                                     >  1:2371438/1‑69 (MQ=2)
                                   gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                              >  1:3613271/1‑57 (MQ=255)
                                    aTGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc                               <  1:1259442/71‑1 (MQ=255)
                                       cggcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGa                            <  1:2639664/71‑1 (MQ=255)
                                         gcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATa                          <  1:2223766/71‑1 (MQ=255)
                                          cgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGAt                           >  1:1122297/1‑69 (MQ=255)
                                          cgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaa                         <  1:1004566/71‑1 (MQ=255)
                                           gTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg                                                >  1:981004/1‑47 (MQ=255)
                                             gAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                      >  1:2682595/1‑71 (MQ=255)
                                               aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                      >  1:2712385/1‑69 (MQ=255)
                                               aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                      >  1:224174/1‑69 (MQ=255)
                                               aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt                     <  1:546197/70‑1 (MQ=255)
                                               aCGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTa                    <  1:1516463/71‑1 (MQ=255)
                                                cGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg                   >  1:259278/1‑71 (MQ=255)
                                                     tATCCGGCCTACATAATCACGCAAATTCAATATATTg                                                >  1:2219736/1‑37 (MQ=255)
                                                      aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt                                         >  1:29821/1‑43 (MQ=255)
                                                      aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGAt                                         >  1:805355/1‑43 (MQ=255)
                                                       tCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa            <  1:3257019/71‑1 (MQ=255)
                                                       tCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa            <  1:381398/71‑1 (MQ=255)
                                                        ccGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg           >  1:987823/1‑71 (MQ=255)
                                                         cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg           <  1:1960297/70‑1 (MQ=255)
                                                           gCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGc         >  1:3932473/1‑70 (MQ=255)
                                                            ccTACATAATCACGCAAATTCAATATATTGCAGAGATc                                        >  1:811427/1‑38 (MQ=255)
                                                            ccTACATAATCACGCAAATTCAATATATTGCAGAGATc                                        >  1:130741/1‑38 (MQ=255)
                                                                 aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa            >  1:2664524/1‑61 (MQ=255)
                                                                  tAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGt  >  1:1169558/1‑70 (MQ=255)
                                                                  |                                                                     
TTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGT  >  minE/2870784‑2870919

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: