Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,433,213 T→C 13.8% I77I (ATT→ATC rhlB → ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,433,2130TC13.8% 228.1 / 26.2 87I77I (ATT→ATCrhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base C (7/5);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGC  >  minE/2433145‑2433277
                                                                    |                                                                
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                >  1:632232/1‑71 (MQ=255)
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                >  1:2169582/1‑71 (MQ=255)
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                >  1:2783935/1‑71 (MQ=255)
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                >  1:4043013/1‑71 (MQ=255)
aCCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGc                                                                >  1:3782836/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:1983101/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:1937436/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                         <  1:2740291/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                         <  1:3142846/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                         <  1:4112268/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                         <  1:1928972/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                         <  1:3914165/63‑1 (MQ=255)
               aCGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGc                                                         <  1:2237003/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:1618824/64‑2 (MQ=38)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                       <  1:2176003/64‑2 (MQ=38)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                       <  1:526147/64‑2 (MQ=38)
                cGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                       <  1:503757/64‑2 (MQ=38)
                 gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                            >  1:4069336/1‑58 (MQ=255)
                 gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                            >  1:3869404/1‑58 (MQ=255)
                 gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                            >  1:2542076/1‑58 (MQ=255)
                 gATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCag                                                       <  1:2292525/63‑2 (MQ=38)
                 gATGGAGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGAt                                                            >  1:1677702/1‑58 (MQ=255)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                     >  1:3425183/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                     >  1:3372416/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                     >  1:51362/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                     >  1:1326826/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                     >  1:797047/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                     >  1:2122581/1‑61 (MQ=37)
                  aTGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATCGGCAATCGCagga                                                     >  1:1344618/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:1922770/1‑71 (MQ=255)
                          tCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTg                                      >  1:2263602/1‑71 (MQ=255)
                            ttACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                  <  1:4004773/57‑1 (MQ=255)
                            ttACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGa                                                  <  1:156789/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:2135128/62‑1 (MQ=255)
                                 tCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCcgcg                                        <  1:1051499/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: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:184090/1‑71 (MQ=255)
                                                 ttctCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGt               >  1:1394298/1‑71 (MQ=255)
                                                 ttctCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGt               >  1:1590942/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:3333248/1‑54 (MQ=255)
                                                             cTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGa                    >  1:976008/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)
                                                                    |                                                                
ACCGGTACCGGGAAAACGATGGCGTTTCTTACGTCAACGTTTCATTATCTTCTCTCTCATCCTGCGATTGCCGATCGCAAGGTGAATCAGCCGCGTGCCTTAATTATGGCACCGACGCGTGAACTTGCCGTGC  >  minE/2433145‑2433277

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: