Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,316,796 G→A 28.0% L157L (CTG→CTA glpX → fructose 1,6‑bisphosphatase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,316,7960GA28.0% 124.6 / 56.2 82L157L (CTG→CTAglpXfructose 1,6‑bisphosphatase II
Reads supporting (aligned to +/- strand):  ref base G (36/23);  new base A (12/11);  total (48/34)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

TCTGAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTACGC  >  minE/2316727‑2316865
                                                                     |                                                                     
tCTGAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGa                                                                      >  1:3436387/1‑71 (MQ=255)
 cTGAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGAc                                                                     >  1:892967/1‑71 (MQ=255)
 cTGAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGAc                                                                     >  1:2561043/1‑71 (MQ=255)
   gAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACgg                                                                   >  1:2211734/1‑71 (MQ=255)
   gAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACgg                                                                   >  1:1916248/1‑71 (MQ=255)
    aaCCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGt                                                                  <  1:2655145/71‑1 (MQ=255)
     aCCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGt                                                                  >  1:1628046/1‑70 (MQ=255)
           cGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGAtt                                                           >  1:110345/1‑71 (MQ=255)
           cGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGAtt                                                           >  1:2443667/1‑71 (MQ=255)
                gCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAAAGATTCTg                                                        <  1:3063080/69‑1 (MQ=255)
                        gcTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACTGTAACGATTCTGGCTAAAccac                                              <  1:974154/70‑1 (MQ=255)
                         cTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCt                                                         <  1:2286159/59‑1 (MQ=255)
                         cTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacg                                             >  1:2570322/1‑71 (MQ=255)
                         cTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacg                                             >  1:2375150/1‑71 (MQ=255)
                          tGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgc                                            >  1:1277261/1‑71 (MQ=255)
                          tGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgc                                            >  1:1877128/1‑71 (MQ=255)
                           gcgcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTc                                                          <  1:3518826/56‑1 (MQ=255)
                           gcgcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTc                                                          <  1:3093711/56‑1 (MQ=255)
                           gcgcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTc                                                          <  1:2554437/56‑1 (MQ=255)
                           gcgcAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAAccacgc                                            >  1:491338/1‑70 (MQ=255)
                              cAATGTAGCGTCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgccacg                                        <  1:1444172/71‑1 (MQ=255)
                              cAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgccacg                                        <  1:1877771/71‑1 (MQ=255)
                               aaTGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGa                                       <  1:456029/71‑1 (MQ=255)
                                  gTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTa                                                    <  1:1549208/55‑1 (MQ=255)
                                   tAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTg                                                                       <  1:323880/35‑1 (MQ=255)
                                   tAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGcc                                   >  1:2601866/1‑71 (MQ=255)
                                   tAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGcc                                   >  1:344055/1‑71 (MQ=255)
                                   tAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGcc                                   >  1:1420237/1‑71 (MQ=255)
                                     gcggcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGt                                 >  1:319946/1‑71 (MQ=255)
                                       ggcggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTa                               <  1:3055579/71‑1 (MQ=255)
                                          ggcgCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCg                            >  1:3265691/1‑71 (MQ=255)
                                            cgcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                          >  1:1078392/1‑71 (MQ=255)
                                            cgcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                          <  1:2639276/71‑1 (MQ=255)
                                             gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgccac                                         <  1:2265059/55‑1 (MQ=255)
                                             gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTAt                              <  1:3807715/66‑1 (MQ=255)
                                             gcTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                          >  1:1894089/1‑70 (MQ=255)
                                             gcTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                          >  1:2265844/1‑70 (MQ=255)
                                             gcTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                          >  1:1225925/1‑70 (MQ=255)
                                              cTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGa                        <  1:261243/71‑1 (MQ=255)
                                              cTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGa                        <  1:1637096/71‑1 (MQ=255)
                                              cTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCt                          >  1:1215383/1‑69 (MQ=255)
                                              cTCGGCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGa                        <  1:1334058/71‑1 (MQ=255)
                                               tCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGa                        >  1:780747/1‑70 (MQ=255)
                                               tCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGa                        >  1:3832451/1‑70 (MQ=255)
                                                 ggCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTaaa                                                  >  1:1656290/1‑42 (MQ=255)
                                                 ggCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGaaa                      >  1:2788309/1‑70 (MQ=255)
                                                 ggCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGaaa                      >  1:2433399/1‑70 (MQ=255)
                                                 ggCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAAt                     >  1:2734182/1‑71 (MQ=255)
                                                 ggCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAAt                     >  1:1133771/1‑71 (MQ=255)
                                                  gCAAACCGTTTAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTa                               <  1:1404140/60‑1 (MQ=255)
                                                  gCAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                    <  1:38744/71‑1 (MQ=255)
                                                   cAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAAccacgc                                            <  1:2193029/46‑1 (MQ=255)
                                                   cAAACCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAAccacgc                                            <  1:2197482/46‑1 (MQ=255)
                                                    aaaCCGTTGAGCGAACTAACGGTAACGATTCTGGCTaaa                                                  >  1:187352/1‑39 (MQ=255)
                                                    aaaCCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATc                             >  1:3767131/1‑60 (MQ=255)
                                                    aaaCCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATc                             >  1:1884678/1‑60 (MQ=255)
                                                     aaCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgccac                                         >  1:1166062/1‑47 (MQ=255)
                                                     aaCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgccac                                         >  1:3239719/1‑47 (MQ=255)
                                                     aaCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                    >  1:1505109/1‑68 (MQ=255)
                                                     aaCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATg                    >  1:1756934/1‑68 (MQ=255)
                                                      aCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcagc                <  1:2594621/71‑1 (MQ=255)
                                                      aCCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcagc                <  1:3658272/71‑1 (MQ=255)
                                                      aCCGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTa                               >  1:59736/1‑56 (MQ=255)
                                                       ccGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATgcagca               <  1:2073655/71‑1 (MQ=255)
                                                       ccGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCAAGCCACGATGCCGTTATCGCTGAAATgcagc                >  1:3214498/1‑70 (MQ=255)
                                                        ctttGAGCGAACTGACGGTAACGATTCTGGCTAAAccacgccac                                         <  1:1825643/42‑1 (MQ=255)
                                                        cGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCaa              <  1:1914549/71‑1 (MQ=255)
                                                        cGTTGAGCGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGcc                                   <  1:1747793/50‑1 (MQ=255)
                                                         gTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAAc             <  1:87219/71‑1 (MQ=255)
                                                         gTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAAc             >  1:3809165/1‑71 (MQ=255)
                                                           tGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTc           >  1:2099641/1‑71 (MQ=255)
                                                             aGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCg                                  <  1:747231/46‑1 (MQ=255)
                                                             aGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCg                                  <  1:2353133/46‑1 (MQ=255)
                                                               cGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCgg         >  1:3910196/1‑69 (MQ=255)
                                                               cGAACTAACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTc           >  1:193206/1‑67 (MQ=255)
                                                                 aaCTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGt      <  1:536281/70‑1 (MQ=255)
                                                                 aaCTGACGGTAACGATTCTGGCTAAACCACGCCACAATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTa     >  1:410771/1‑71 (MQ=255)
                                                                  aCTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTAc    <  1:2849166/71‑1 (MQ=255)
                                                                  aCTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTAACGCTGAAATGCAGCAACTCg          >  1:3499733/1‑65 (MQ=255)
                                                                  aCTAACGGTAACGATTCTGGCTAAACCACGCCACGATg                                     <  1:1494978/38‑1 (MQ=255)
                                                                   cTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTAcg   >  1:2074378/1‑71 (MQ=255)
                                                                   cTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTAcg   >  1:2693778/1‑71 (MQ=255)
                                                                    tGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTAcgc  >  1:151240/1‑71 (MQ=255)
                                                                     |                                                                     
TCTGAACCTGCCGCTGGCGGATAACCTGCGCAATGTAGCGGCGGCGCTCGGCAAACCGTTGAGCGAACTGACGGTAACGATTCTGGCTAAACCACGCCACGATGCCGTTATCGCTGAAATGCAGCAACTCGGCGTACGC  >  minE/2316727‑2316865

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: