Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,433,215 C→G 14.0% A78G (GCC→GGC)  rhlB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,433,2150CG14.0% 217.2 / 21.5 85A78G (GCC→GGC) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (32/41);  new base G (7/5);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAGATCCA  >  minE/2433145‑2433284
                                                                      |                                                                     
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                       >  1:632232/1‑71 (MQ=255)
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                       >  1:2783935/1‑71 (MQ=255)
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                       >  1:3782836/1‑71 (MQ=255)
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                       >  1:2169582/1‑71 (MQ=255)
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                       >  1:4043013/1‑71 (MQ=255)
   ggTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGcc                                                                      >  1:838174/1‑69 (MQ=255)
     tACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGcc                                                                      >  1:2306952/1‑67 (MQ=255)
       ccGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCg                                                                 <  1:2416401/70‑1 (MQ=255)
       ccGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCg                                                                 <  1:3332141/70‑1 (MQ=255)
       ccGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCg                                                                 <  1:1983101/70‑1 (MQ=255)
             aaaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTg                                                          <  1:1148967/71‑1 (MQ=255)
              aaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:1119492/64‑1 (MQ=255)
              aaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:1674833/64‑1 (MQ=255)
              aaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:2029330/64‑1 (MQ=255)
              aaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:2175360/71‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCGTGCGATTGCCg                                                                     <  1:3030862/58‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCg                                                                     <  1:4123974/58‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCg                                                                     <  1:3323117/58‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCg                                                                     <  1:620355/58‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:1928972/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:1937436/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:3142846/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:3914165/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:2237003/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:4112268/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:2740291/63‑1 (MQ=255)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:3432198/62‑1 (MQ=255)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAAt                                                       >  1:365919/1‑71 (MQ=255)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAAt                                                       >  1:2385536/1‑71 (MQ=255)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAAt                                                       >  1:3887193/1‑71 (MQ=255)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:526147/64‑2 (MQ=38)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:503757/64‑2 (MQ=38)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:1618824/64‑2 (MQ=38)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:2176003/64‑2 (MQ=38)
                 gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                                   >  1:2542076/1‑58 (MQ=255)
                 gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                                   >  1:4069336/1‑58 (MQ=255)
                 gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                                   >  1:3869404/1‑58 (MQ=255)
                 gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:2292525/63‑2 (MQ=38)
                 gATGGAGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                                   >  1:1677702/1‑58 (MQ=255)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:3372416/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:2122581/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:1344618/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:51362/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:1326826/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:3425183/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:797047/1‑61 (MQ=37)
                  aTGGCGTCTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGt                                                           <  1:3736961/65‑1 (MQ=255)
                     gCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGc                                                   <  1:2662854/70‑1 (MQ=255)
                     gCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGc                                                   <  1:13277/70‑1 (MQ=255)
                      cGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGc                                                   <  1:3777667/69‑1 (MQ=255)
                          tCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTg                                             >  1:1922770/1‑71 (MQ=255)
                          tCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTg                                             >  1:2263602/1‑71 (MQ=255)
                            ttACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:156789/57‑1 (MQ=255)
                            ttACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:4004773/57‑1 (MQ=255)
                            ttACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:2180067/57‑1 (MQ=255)
                            tcaCGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:288735/55‑1 (MQ=255)
                               cGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTa                                        <  1:2023230/71‑1 (MQ=255)
                                 tCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCcgcg                                               <  1:1051499/62‑1 (MQ=255)
                                 tCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCcgcg                                               <  1:2135128/62‑1 (MQ=255)
                                 tCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCcgcg                                               <  1:2248576/62‑1 (MQ=255)
                                   aaCGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                       >  1:2187611/1‑36 (MQ=255)
                                   aaCGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTa                                    >  1:2907261/1‑71 (MQ=255)
                                    aCGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTAt                                   <  1:1194000/71‑1 (MQ=255)
                                    aCGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTAt                                   >  1:3175961/1‑71 (MQ=255)
                                     cGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATg                                  >  1:570603/1‑71 (MQ=255)
                                      gTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATgg                                 <  1:2265768/71‑1 (MQ=255)
                                      gTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATgg                                 <  1:709049/71‑1 (MQ=255)
                                       tttCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:749032/39‑1 (MQ=255)
                                       tttCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATgg                                 >  1:1534492/1‑70 (MQ=255)
                                       tttCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATgg                                 >  1:1655707/1‑70 (MQ=255)
                                         tCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCAc                              >  1:3100541/1‑71 (MQ=255)
                                         tCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCAc                              >  1:1460130/1‑71 (MQ=255)
                                           aTTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAAt                                      <  1:321543/61‑1 (MQ=255)
                                               tcttctCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGAcgc                        >  1:4125199/1‑71 (MQ=255)
                                                 ttctCTCTCATCGTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGt                      >  1:1199112/1‑71 (MQ=255)
                                                 ttctCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGt                      >  1:184090/1‑71 (MQ=255)
                                                 ttctCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGt                      >  1:1590942/1‑71 (MQ=255)
                                                 ttctCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGt                      >  1:1394298/1‑71 (MQ=255)
                                                   ctctctcATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTg                     <  1:2772054/70‑1 (MQ=255)
                                                        tcATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCtt                                         <  1:2875846/45‑1 (MQ=255)
                                                        tcATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTg               >  1:2113479/1‑71 (MQ=255)
                                                         cATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGc              >  1:1241256/1‑71 (MQ=255)
                                                             cTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGa                           >  1:976008/1‑54 (MQ=255)
                                                             cTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGa                           >  1:3333248/1‑54 (MQ=255)
                                                                cGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTg          >  1:1903904/1‑68 (MQ=255)
                                                                cGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTg          >  1:1269337/1‑68 (MQ=255)
                                                                cGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGc         <  1:3896796/69‑1 (MQ=255)
                                                                    tGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTg                     <  1:2389163/53‑1 (MQ=255)
                                                                     gCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCa                               >  1:4030929/1‑42 (MQ=255)
                                                                     gCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCa                               >  1:1797614/1‑42 (MQ=255)
                                                                     gCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCa                               >  1:2950956/1‑42 (MQ=255)
                                                                     gCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCa                               >  1:2787724/1‑42 (MQ=255)
                                                                     gCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAGATCCa  >  1:2984549/1‑71 (MQ=255)
                                                                      ccGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCTTGCAg       <  1:1248963/65‑1 (MQ=255)
                                                                      ccGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAg       <  1:1798319/65‑1 (MQ=255)
                                                                      ccGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAg       <  1:374099/65‑1 (MQ=255)
                                                                      ccGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAg       <  1:3532566/65‑1 (MQ=255)
                                                                      |                                                                     
ACCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAGATCCA  >  minE/2433145‑2433284

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: