Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,870,827 G→A 8.5% 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,8270GA 8.5% 144.8 / 2.2 47intergenic (+70/‑110)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (26/17);  new base A (2/2);  total (28/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 = 1.00e+00

ATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAA  >  minE/2870757‑2870896
                                                                      |                                                                     
aTCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATgcggcg                                                                       <  1:927122/71‑1 (MQ=255)
    aTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   >  1:1418681/1‑71 (MQ=255)
     tACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   >  1:1289810/1‑70 (MQ=255)
     tACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGa                                                                    <  1:3888894/69‑1 (MQ=255)
      aCCCTCTTGAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                 >  1:411009/1‑71 (MQ=255)
       cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                 <  1:2692973/70‑1 (MQ=255)
       cccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                 <  1:3878218/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)
               aaaaGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   >  1:1852901/1‑60 (MQ=255)
                 aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                      <  1:2062498/71‑1 (MQ=255)
                     agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                      >  1:3865793/1‑67 (MQ=255)
                     agagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                   >  1:3267578/1‑70 (MQ=255)
                          gTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:1402682/47‑1 (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)
                                      aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTg                                                                     <  1:2652678/35‑1 (MQ=255)
                                      aCAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTc                                  <  1:2085475/70‑1 (MQ=255)
                                       cAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                   >  1:1911741/1‑36 (MQ=40)
                                        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:594597/1‑68 (MQ=255)
                                                 tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg                         >  1:1111689/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)
                                                                      |                                                                     
ATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAA  >  minE/2870757‑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: