Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,870,833 G→A 8.3% intergenic (+76/‑104) 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,8330GA 8.3% 143.9 / 3.0 48intergenic (+76/‑104)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (27/17);  new base A (2/2);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

ACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAG  >  minE/2870763‑2870902
                                                                      |                                                                     
aCCCTCTTGAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                       >  1:411009/1‑71 (MQ=255)
 cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                       <  1:3878218/70‑1 (MQ=255)
 cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                       <  1:2692973/70‑1 (MQ=255)
 cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                      <  1:2924083/71‑1 (MQ=255)
 cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                      <  1:2308556/71‑1 (MQ=255)
    tctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                    >  1:1742040/1‑70 (MQ=255)
     ctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                  >  1:1639825/1‑71 (MQ=255)
           aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                            <  1:2062498/71‑1 (MQ=255)
               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                            >  1:3865793/1‑67 (MQ=255)
               agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                         >  1:3267578/1‑70 (MQ=255)
                     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)
                                  aaaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                      >  1:2913063/1‑38 (MQ=40)
                                   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:2753095/1‑71 (MQ=255)
                                                 cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa                      >  1:1742993/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)
                                                                      |                                                                     
ACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAG  >  minE/2870763‑2870902

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: