Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,433,217 G→A 14.1% D79N (GAT→AAT)  rhlB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,433,2170GA14.1% 234.4 / 25.4 85D79N (GAT→AAT) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (32/41);  new base A (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

CCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAGATCCA  >  minE/2433152‑2433284
                                                                 |                                                                   
ccGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCg                                                                 <  1:1983101/70‑1 (MQ=255)
ccGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCg                                                                 <  1:2416401/70‑1 (MQ=255)
ccGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCg                                                                 <  1:3332141/70‑1 (MQ=255)
      aaaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTg                                                          <  1:1148967/71‑1 (MQ=255)
       aaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:1119492/64‑1 (MQ=255)
       aaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:2029330/64‑1 (MQ=255)
       aaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:1674833/64‑1 (MQ=255)
       aaCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:2175360/71‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCGTGCGATTGCCg                                                                     <  1:3030862/58‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCg                                                                     <  1:620355/58‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCg                                                                     <  1:4123974/58‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCg                                                                     <  1:3323117/58‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:3142846/63‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:1937436/63‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:1928972/63‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:2237003/63‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:2740291/63‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:3914165/63‑1 (MQ=255)
        aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:4112268/63‑1 (MQ=255)
         cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                                <  1:3432198/62‑1 (MQ=255)
         cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAAt                                                       >  1:2385536/1‑71 (MQ=255)
         cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAAt                                                       >  1:365919/1‑71 (MQ=255)
         cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAAt                                                       >  1:3887193/1‑71 (MQ=255)
         cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:503757/64‑2 (MQ=38)
         cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:1618824/64‑2 (MQ=38)
         cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:526147/64‑2 (MQ=38)
         cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:2176003/64‑2 (MQ=38)
          gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                                   >  1:2542076/1‑58 (MQ=255)
          gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                                   >  1:3869404/1‑58 (MQ=255)
          gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                                   >  1:4069336/1‑58 (MQ=255)
          gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                              <  1:2292525/63‑2 (MQ=38)
          gATGGAGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                                   >  1:1677702/1‑58 (MQ=255)
           aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:1344618/1‑61 (MQ=37)
           aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:51362/1‑61 (MQ=37)
           aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:1326826/1‑61 (MQ=37)
           aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:3372416/1‑61 (MQ=37)
           aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:3425183/1‑61 (MQ=37)
           aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:797047/1‑61 (MQ=37)
           aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                            >  1:2122581/1‑61 (MQ=37)
           aTGGCGTCTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGt                                                           <  1:3736961/65‑1 (MQ=255)
              gCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGc                                                   <  1:13277/70‑1 (MQ=255)
              gCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGc                                                   <  1:2662854/70‑1 (MQ=255)
               cGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGc                                                   <  1:3777667/69‑1 (MQ=255)
                   tCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTg                                             >  1:2263602/1‑71 (MQ=255)
                   tCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTg                                             >  1:1922770/1‑71 (MQ=255)
                     ttACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:4004773/57‑1 (MQ=255)
                     ttACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:2180067/57‑1 (MQ=255)
                     ttACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                         <  1:156789/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)
                            aaCGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTa                                    >  1:2907261/1‑71 (MQ=255)
                             aCGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTAt                                   >  1:3175961/1‑71 (MQ=255)
                             aCGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTAt                                   <  1:1194000/71‑1 (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:1655707/1‑70 (MQ=255)
                                tttCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATgg                                 >  1:1534492/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:1590942/1‑71 (MQ=255)
                                          ttctCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGt                      >  1:184090/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:1269337/1‑68 (MQ=255)
                                                         cGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTg          >  1:1903904/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:2950956/1‑42 (MQ=255)
                                                              gCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCa                               >  1:2787724/1‑42 (MQ=255)
                                                              gCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCa                               >  1:1797614/1‑42 (MQ=255)
                                                              gCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAGATCCa  >  1:2984549/1‑71 (MQ=255)
                                                               ccGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCTTGCAg       <  1:1248963/65‑1 (MQ=255)
                                                               ccGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAg       <  1:3532566/65‑1 (MQ=255)
                                                               ccGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAg       <  1:374099/65‑1 (MQ=255)
                                                               ccGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAg       <  1:1798319/65‑1 (MQ=255)
                                                                 gATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATgg                                 <  1:1839539/37‑1 (MQ=255)
                                                                 |                                                                   
CCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGCAGATCCA  >  minE/2433152‑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: