Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,292,617 A→G 19.7% intergenic (+168/‑253) fsaB → / → gldA fructose‑6‑phosphate aldolase 2/glycerol dehydrogenase, NAD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,292,6170AG19.7% 175.5 / 37.7 81intergenic (+168/‑253)fsaB/gldAfructose‑6‑phosphate aldolase 2/glycerol dehydrogenase, NAD
Reads supporting (aligned to +/- strand):  ref base A (40/25);  new base G (6/10);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 9.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTAGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGTGCTGGC  >  minE/2292548‑2292687
                                                                     |                                                                      
ttattaGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGaaa                                                                       >  1:3648993/1‑71 (MQ=255)
    taGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGaaa                                                                       >  1:683980/1‑67 (MQ=255)
    taGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACgg                                                                    >  1:93785/1‑70 (MQ=255)
    taGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACgg                                                                    >  1:3727111/1‑70 (MQ=255)
          cACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACgg                                                                    <  1:3461354/64‑1 (MQ=255)
            cTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTAc                                                           >  1:2666191/1‑71 (MQ=255)
              ggCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGaaa                                                                       >  1:2832008/1‑57 (MQ=255)
               gCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGt                                                        >  1:2661299/1‑71 (MQ=255)
               gCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGt                                                        >  1:2344065/1‑71 (MQ=255)
                 aGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGtt                                                       >  1:507774/1‑70 (MQ=255)
                 aGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTc                                                      <  1:499337/71‑1 (MQ=255)
                 aGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTc                                                      <  1:481374/71‑1 (MQ=255)
                 aGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGt                                                             <  1:2465507/64‑1 (MQ=255)
                  ggggCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCGAGGGCGGAAACGGCATTCGTACGGTTCa                                                     <  1:24181/71‑1 (MQ=255)
                  ggggCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGaa                                                                        <  1:2307532/52‑1 (MQ=255)
                   gggCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAg                                                    <  1:853203/71‑1 (MQ=255)
                     gCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAgg                                                   <  1:2373379/70‑1 (MQ=255)
                         aaTACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTg                                              >  1:3108216/1‑71 (MQ=255)
                          aTACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGCCGGCATTCGTACGGTTCAGGAGCTg                                              >  1:3852332/1‑70 (MQ=255)
                          aTACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTg                                              >  1:618422/1‑70 (MQ=255)
                           tACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCa                                            <  1:393427/71‑1 (MQ=255)
                           tACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACgg                                                                    >  1:1699033/1‑47 (MQ=255)
                           tACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCa                                                     >  1:819400/1‑62 (MQ=255)
                           tACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCa                                            <  1:1216311/71‑1 (MQ=255)
                             cGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCaaa                                          >  1:2216173/1‑71 (MQ=255)
                             cGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAg                                                    >  1:3135785/1‑61 (MQ=255)
                             cGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCaaa                                          >  1:3341952/1‑71 (MQ=255)
                             cGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCaa                                           >  1:246607/1‑70 (MQ=255)
                              gTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACgg                                                                    <  1:2660297/44‑1 (MQ=255)
                               ttGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACg                                        <  1:3718779/71‑1 (MQ=255)
                                   tCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGt                                    >  1:821448/1‑71 (MQ=255)
                                   tCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGt                                    >  1:889418/1‑71 (MQ=255)
                                    ccGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGt                                    >  1:261353/1‑70 (MQ=255)
                                     cGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCt                                               >  1:361859/1‑58 (MQ=255)
                                     cGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCaaa                                          >  1:3318744/1‑63 (MQ=255)
                                     cGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAAc                                         <  1:2395850/64‑1 (MQ=255)
                                     cGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAAc                                         <  1:2040330/64‑1 (MQ=255)
                                     cGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTa                                  <  1:3295275/71‑1 (MQ=255)
                                       tATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTg                                              >  1:1829074/1‑57 (MQ=255)
                                       tATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGa                                <  1:12702/71‑1 (MQ=255)
                                        aTGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGa                                >  1:2293092/1‑70 (MQ=255)
                                         tGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaaa                              >  1:1798882/1‑71 (MQ=255)
                                         tGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaaa                              >  1:3879934/1‑71 (MQ=255)
                                         tGTTAACCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaa                               <  1:2829435/70‑1 (MQ=255)
                                         tGTTAACCGCGTAGATGCCCAGGGCGGAACCGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaaa                              >  1:3496643/1‑71 (MQ=255)
                                         tGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaaa                              >  1:1255266/1‑71 (MQ=255)
                                         tGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaa                               <  1:717588/70‑1 (MQ=255)
                                           ttAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATg                            >  1:3635975/1‑71 (MQ=255)
                                           ttAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATg                            >  1:1581767/1‑71 (MQ=255)
                                            tAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATg                            >  1:2429958/1‑70 (MQ=255)
                                            tAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATg                            >  1:3764414/1‑70 (MQ=255)
                                            tAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATg                            >  1:3612614/1‑70 (MQ=255)
                                             aaCCGCGTAGATGCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCa                          >  1:1458257/1‑71 (MQ=255)
                                             aaCCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCa                          >  1:1634653/1‑71 (MQ=255)
                                             aaCCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCa                          >  1:2673637/1‑71 (MQ=255)
                                              aCCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCa                          <  1:836508/70‑1 (MQ=255)
                                              aCCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCa                          <  1:1461869/70‑1 (MQ=255)
                                              aCCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCAc                         >  1:2969812/1‑71 (MQ=255)
                                                cgcgTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCt                                               >  1:345261/1‑47 (MQ=255)
                                                   gTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAgg                                                   <  1:2378310/40‑1 (MQ=255)
                                                   gTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGtt                                   >  1:356844/1‑56 (MQ=255)
                                                   gTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGtt                                   >  1:2620806/1‑56 (MQ=255)
                                                       aTGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGa                 >  1:1151709/1‑70 (MQ=255)
                                                        tGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGa                 >  1:79087/1‑69 (MQ=255)
                                                         gCCCAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCa                                            <  1:4044879/41‑1 (MQ=255)
                                                         gCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaaa                              >  1:3647068/1‑55 (MQ=255)
                                                         gCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAg              >  1:160586/1‑71 (MQ=255)
                                                         gCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAg              >  1:3518577/1‑71 (MQ=255)
                                                          cccAGGGCGGAGACGTCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAg              <  1:2300967/70‑1 (MQ=255)
                                                          cccAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAg              <  1:3080600/70‑1 (MQ=255)
                                                          cccAGGGCGGAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAg              <  1:1068591/70‑1 (MQ=255)
                                                          cccAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATg                            <  1:470855/56‑1 (MQ=255)
                                                            cAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGGGCCAGAAAGCa            <  1:2401974/70‑1 (MQ=255)
                                                            cAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCAt           <  1:3347259/71‑1 (MQ=255)
                                                              gggCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATgg         <  1:217435/71‑1 (MQ=255)
                                                              gggCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATgg         >  1:395663/1‑71 (MQ=255)
                                                               ggCGGAAACGGCATTCGTACGGTTCAGGAGCTGCa                                            <  1:3173565/35‑1 (MQ=255)
                                                               ggCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGt        <  1:3288979/71‑1 (MQ=255)
                                                               ggCGGAAACGGCATTCGTACCGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGt        <  1:4009098/71‑1 (MQ=255)
                                                                gCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGa                                <  1:2317586/46‑1 (MQ=255)
                                                                 cGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGaaa                              >  1:2592983/1‑47 (MQ=255)
                                                                  ggAGACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGa                                <  1:4073038/44‑1 (MQ=255)
                                                                   gAAACTGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGa                                <  1:1060699/43‑1 (MQ=255)
                                                                   gAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGa                                <  1:1908272/43‑1 (MQ=255)
                                                                   gAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCACAAAGCATGGTGCTg    >  1:606021/1‑71 (MQ=255)
                                                                     aaCGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGTGCTGGc  <  1:1279549/71‑1 (MQ=255)
                                                                     aaCGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCCCCAGAAAGCATGGTGCTGGc  <  1:917869/71‑1 (MQ=255)
                                                                     |                                                                      
TTATTAGCCGCACTGGCAGGGGCAAAATACGTTGCTCCGTATGTTAACCGCGTAGATGCCCAGGGCGGAAACGGCATTCGTACGGTTCAGGAGCTGCAAACGCTGTTAGAAATGCACGCGCCAGAAAGCATGGTGCTGGC  >  minE/2292548‑2292687

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: