Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,870,850 T→A 15.7% 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,8500TA15.7% 264.9 / 38.5 102intergenic (+93/‑87)rnr/rlmBexoribonuclease R, RNase R/23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (41/45);  new base A (10/6);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 107 total aligned reads displayed.

GGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGT  >  minE/2870782‑2870919
                                                                    |                                                                     
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:2399343‑M2/33‑70 (MQ=255)
   gacgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATc                                                                   <  1:3374261‑M2/38‑1 (MQ=255)
    acgcgttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCAc                                                                 >  1:675455‑M2/32‑71 (MQ=255)
        tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                                              >  1:2262234/1‑70 (MQ=255)
        tGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCa                                                              >  1:1983495/1‑70 (MQ=255)
        gttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGc                                                               >  1:2794429‑M2/28‑69 (MQ=255)
          cTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCaaa                                                            <  1:2552033/70‑1 (MQ=255)
           aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaat                                                          >  1:3757004‑M2/25‑69 (MQ=255)
           aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaat                                                          >  1:529918‑M2/25‑69 (MQ=255)
           aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaat                                                          >  1:1651814‑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)
             gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaa                                                           <  1:2349645‑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)
               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:714569/71‑1 (MQ=255)
                       cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                              <  1:1643017/71‑1 (MQ=255)
                       cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                              <  1:3756211/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)
                         ttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAatatat                                                  >  1:1006746/1‑65 (MQ=255)
                          tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                              >  1:3084404/1‑68 (MQ=255)
                          tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                              >  1:138408/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)
                                     gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc                                >  1:3081884/1‑71 (MQ=255)
                                     gATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGc                                >  1:3537357/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:920698/71‑1 (MQ=255)
                                            cgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaa                         <  1:3122007/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:193324/39‑1 (MQ=255)
                                                aaCGCCTTATCCGGCCTACATAATCACGCAAATTCAata                                                     <  1:3901085/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:3766914/1‑70 (MQ=255)
                                                        aTCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCAt              >  1:351908/1‑70 (MQ=255)
                                                           cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTa                                   >  1:316145/1‑46 (MQ=255)
                                                           cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTa                                   >  1:54003/1‑46 (MQ=255)
                                                           cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAgg          <  1:2327750/71‑1 (MQ=255)
                                                           cGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAgg          >  1:118506/1‑71 (MQ=255)
                                                             gCCTACATAATCACGCAAATTCAATATATTGCagag                                           <  1:1811192/36‑1 (MQ=255)
                                                              ccTACATAATCACGCAAATTCAATATATTGCAgaga                                          <  1:2872915/36‑1 (MQ=255)
                                                              ccTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGc         <  1:1989548/69‑1 (MQ=255)
                                                               cTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTa                    >  1:827803/1‑57 (MQ=255)
                                                               cTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCaa       >  1:2539346/1‑70 (MQ=255)
                                                                tACATAATCACGCAAATTCAATATATTGCAGAGAt                                         >  1:2999299/1‑35 (MQ=255)
                                                                tACATAATCACGCAAATTCAATATATTGCAGAGATCa                                       <  1:2603851/37‑1 (MQ=255)
                                                                tACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt                     <  1:666149/55‑1 (MQ=255)
                                                                tACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAg           <  1:2964134/65‑1 (MQ=255)
                                                                  cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGa                            >  1:2113488/1‑46 (MQ=255)
                                                                  cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGa                            >  1:1873332/1‑46 (MQ=255)
                                                                  cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaaa                        >  1:2676823/1‑50 (MQ=255)
                                                                  cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATaaa                        >  1:1361808/1‑50 (MQ=255)
                                                                  cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt                     >  1:183268/1‑53 (MQ=255)
                                                                  cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACg   <  1:703422/71‑1 (MQ=255)
                                                                  cATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACg   >  1:713082/1‑71 (MQ=255)
                                                                   aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACg   >  1:488435/1‑70 (MQ=255)
                                                                   aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACg   >  1:2222457/1‑70 (MQ=255)
                                                                   aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGt  <  1:1593059/71‑1 (MQ=255)
                                                                   aTAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGt  <  1:1588188/71‑1 (MQ=255)
                                                                    tAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACg                      <  1:2184324/50‑1 (MQ=255)
                                                                    tAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCa               >  1:3479401/1‑57 (MQ=255)
                                                                    tAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGt  >  1:293860/1‑70 (MQ=255)
                                                                    tAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGt  >  1:1137408/1‑70 (MQ=255)
                                                                    |                                                                     
GGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGT  >  minE/2870782‑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: