Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,687,159 G→A 26.2% L374L (CTG→TTG)  pck ← phosphoenolpyruvate carboxykinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,687,1590GA26.2% 152.3 / 57.1 84L374L (CTG→TTG) pckphosphoenolpyruvate carboxykinase
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (12/10);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGCCCGC  >  minE/2687094‑2687228
                                                                 |                                                                     
cTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAgg                                                                      >  1:36607/1‑67 (MQ=255)
cTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGc                                                                     >  1:365804/1‑68 (MQ=255)
cTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGc                                                                     >  1:3306874/1‑68 (MQ=255)
cTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGc                                                                     >  1:2713568/1‑68 (MQ=255)
cTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGc                                                                     >  1:232029/1‑68 (MQ=255)
cTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGc                                                                     >  1:1385561/1‑68 (MQ=255)
   aGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGaa                                                                <  1:2293445/70‑1 (MQ=255)
   aGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGaa                                                                <  1:2288659/70‑1 (MQ=255)
     tACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAAcc                                                             <  1:549053/71‑1 (MQ=255)
     tACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAAcc                                                             <  1:185585/71‑1 (MQ=255)
      aCCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcg                                                            >  1:3498336/1‑71 (MQ=255)
      aCCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcg                                                            >  1:3673676/1‑71 (MQ=255)
             aGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCaa                                                     <  1:124997/71‑1 (MQ=255)
             aGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCaa                                                     <  1:2987925/71‑1 (MQ=255)
             aGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCaa                                                     <  1:3052028/71‑1 (MQ=255)
              gTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAac                                                    <  1:116867/71‑1 (MQ=255)
              gTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAac                                                    <  1:2154011/71‑1 (MQ=255)
              gTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAac                                                    <  1:3083181/71‑1 (MQ=255)
                ttGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAacac                                                  >  1:2625619/1‑71 (MQ=255)
                 tGGCGGTGAAGCGAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACg                                                 >  1:386554/1‑71 (MQ=255)
                 tGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACg                                                 >  1:3084067/1‑71 (MQ=255)
                 tGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACg                                                 >  1:1598359/1‑71 (MQ=255)
                 tGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACg                                                 >  1:3048923/1‑71 (MQ=255)
                   gCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAac                                                    <  1:2241583/66‑1 (MQ=255)
                     ggTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGa                                             <  1:3964244/71‑1 (MQ=255)
                     ggTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGa                                             <  1:1324371/71‑1 (MQ=255)
                      gTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGa                                             <  1:1556651/70‑1 (MQ=255)
                      gTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGaa                                            <  1:3898577/71‑1 (MQ=255)
                      gTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGaa                                            <  1:2556903/71‑1 (MQ=255)
                        gAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAg                                          >  1:1539793/1‑71 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcgg                                                           >  1:2713765/1‑53 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcgg                                                           >  1:1088133/1‑53 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcgg                                                           >  1:2313494/1‑53 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcgg                                                           >  1:1739259/1‑53 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcgg                                                           >  1:1926969/1‑53 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAAcc                                                             <  1:93029/51‑1 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAAcc                                                             <  1:204673/51‑1 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAAcc                                                             <  1:4113646/51‑1 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTCGGTTTGATCGGCAGTCAAGCGAGAAAcc                                                             <  1:1866825/51‑1 (MQ=255)
                          aGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCa                                        >  1:2746050/1‑71 (MQ=255)
                          aGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCa                                        >  1:3998399/1‑71 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCa                                        >  1:964634/1‑70 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGaaa                                                               >  1:2822131/1‑47 (MQ=255)
                               gagagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcag                                    >  1:2505281/1‑70 (MQ=255)
                               gagagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcag                                    >  1:433317/1‑70 (MQ=255)
                                 gagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgaga                                                                 <  1:1706823/39‑1 (MQ=255)
                                 gagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgaga                                                                 <  1:2311391/39‑1 (MQ=255)
                                 gagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgaga                                                                 <  1:1323970/39‑1 (MQ=255)
                                 gagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgaga                                                                 <  1:1752432/39‑1 (MQ=255)
                                 gagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCcaga                                                                 <  1:3085877/39‑1 (MQ=255)
                                      aGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                   <  1:2054197/64‑1 (MQ=255)
                                      aGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                   <  1:2285906/64‑1 (MQ=255)
                                      aGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                   <  1:2276499/64‑1 (MQ=255)
                                      aGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAg                             >  1:384883/1‑70 (MQ=255)
                                       gTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAac                                                    <  1:2097657/46‑1 (MQ=255)
                                       gTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAac                                                    <  1:2409918/46‑1 (MQ=255)
                                       gTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAac                                                    <  1:3190513/46‑1 (MQ=255)
                                        tGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAgg                            >  1:3962615/1‑69 (MQ=255)
                                        tGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAgg                            >  1:299025/1‑69 (MQ=255)
                                        tGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGaa                          <  1:1311052/71‑1 (MQ=255)
                                         gATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAgg                            >  1:123239/1‑68 (MQ=255)
                                          aTACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGa                        <  1:233922/71‑1 (MQ=255)
                                             cTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAaca                                                   >  1:2929947/1‑41 (MQ=255)
                                             cTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAaca                                                   >  1:2628215/1‑41 (MQ=255)
                                             cTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAaca                                                   >  1:1296084/1‑41 (MQ=255)
                                              tGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAAc                    >  1:576758/1‑71 (MQ=255)
                                               gggTTTGATCGGCAGTCAAGCGAGAAACCGGAGGCAACACGCCGAAAGCATAAGCAGTCAGGAAGATAAcc                   >  1:2299513/1‑71 (MQ=255)
                                                ggTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAAcc                   <  1:775589/70‑1 (MQ=255)
                                                ggTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCt                  <  1:2502153/71‑1 (MQ=255)
                                                 gTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCTACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCtt                 <  1:3221188/71‑1 (MQ=255)
                                                 gTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAAcc                   <  1:2051769/69‑1 (MQ=255)
                                                 gTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAAcc                   <  1:191082/69‑1 (MQ=255)
                                                  tttGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCt                  >  1:798649/1‑69 (MQ=255)
                                                  tttGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa                <  1:1546122/71‑1 (MQ=255)
                                                   ttGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCtt                 >  1:2857125/1‑69 (MQ=255)
                                                    tGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGt              >  1:2947200/1‑71 (MQ=255)
                                                          gCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGa                        <  1:965060/55‑1 (MQ=255)
                                                              tCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGaa                          <  1:2501467/49‑1 (MQ=255)
                                                              tCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGaa                          <  1:2100417/49‑1 (MQ=255)
                                                              tCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGaa                          <  1:1179141/49‑1 (MQ=255)
                                                                aGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAATATAACCTTAGTCGCGTGGCCCGc  >  1:3276145/1‑71 (MQ=255)
                                                                aGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGCCCg   >  1:1060305/1‑70 (MQ=255)
                                                                aGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGCCCGc  >  1:3340338/1‑71 (MQ=255)
                                                                aGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGCCCGc  >  1:1904152/1‑71 (MQ=255)
                                                                aGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGCCCGc  >  1:1259279/1‑71 (MQ=255)
                                                                 |                                                                     
CTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGCCCGC  >  minE/2687094‑2687228

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