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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,687,1590GA44.0% 40.5 / 120.6 91L374L (CTG→TTG) pckphosphoenolpyruvate carboxykinase
Reads supporting (aligned to +/- strand):  ref base G (13/38);  new base A (22/18);  total (35/56)
Fisher's exact test for biased strand distribution p-value = 5.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGCCC  >  minE/2687094‑2687226
                                                                 |                                                                   
cTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:561974/71‑1 (MQ=255)
cTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:750865/71‑1 (MQ=255)
  cAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgaga                                                               >  1:1133984/1‑70 (MQ=255)
    gTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGaaa                                                             >  1:1111586/1‑70 (MQ=255)
    gTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGaaa                                                             >  1:2895554/1‑70 (MQ=255)
    gTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAAc                                                            >  1:3754960/1‑71 (MQ=255)
       ccGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcg                                                          <  1:973381/70‑1 (MQ=255)
       ccGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcg                                                          <  1:3942382/70‑1 (MQ=255)
       ccGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcg                                                          <  1:3804580/70‑1 (MQ=255)
       ccGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACcgg                                                         >  1:2833800/1‑71 (MQ=255)
       ccGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACcgg                                                         >  1:2251693/1‑71 (MQ=255)
       ccGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACcgg                                                         >  1:1469361/1‑71 (MQ=255)
       ccGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACcgg                                                         >  1:1591338/1‑71 (MQ=255)
       ccGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACcgg                                                         >  1:1635870/1‑71 (MQ=255)
           ccAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCggcggc                                                     <  1:3340194/71‑1 (MQ=255)
           ccAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCggcggc                                                     <  1:2523970/71‑1 (MQ=255)
           ccAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCggcggc                                                     <  1:242808/71‑1 (MQ=255)
                 tGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACg                                               >  1:1219867/1‑71 (MQ=255)
                 tGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAacac                                                <  1:2419294/70‑1 (MQ=255)
                       tGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcg                                                          <  1:3648060/54‑1 (MQ=255)
                        gAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAg                                        >  1:1772570/1‑71 (MQ=255)
                        gAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAg                                        >  1:2412177/1‑71 (MQ=255)
                        gAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAg                                        >  1:3020903/1‑71 (MQ=255)
                         aaGCCAGAGAGGATGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAg                                        <  1:3406057/70‑1 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGaaa                                         >  1:3194362/1‑69 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAg                                        >  1:3713969/1‑70 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAg                                        >  1:2816823/1‑70 (MQ=255)
                         aaGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAg                                        >  1:224650/1‑70 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:907111/44‑1 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:3048064/44‑1 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:2298792/44‑1 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:3093941/44‑1 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:3881240/44‑1 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:3732947/44‑1 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:371440/44‑1 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCgag                                                                <  1:1172932/44‑1 (MQ=255)
                           gCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCAt                                     <  1:2962810/71‑1 (MQ=255)
                             cAGAGAGGAAGTGATACTGTGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCAt                                     >  1:1416880/1‑69 (MQ=255)
                             cAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCAt                                     >  1:3175741/1‑69 (MQ=255)
                               gagagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                 <  1:2433534/71‑1 (MQ=255)
                               gagagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                 <  1:2438668/71‑1 (MQ=255)
                               gagagGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                 <  1:3483622/71‑1 (MQ=255)
                               gagagGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                 <  1:1123489/71‑1 (MQ=255)
                               gagagGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                 <  1:2563812/71‑1 (MQ=255)
                                agagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACcggcgg                                                      <  1:2002297/49‑1 (MQ=255)
                                agagGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcag                                  >  1:171/1‑69 (MQ=255)
                                agagGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCaa                                                   >  1:2524481/1‑52 (MQ=255)
                                agagGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCaa                                                   >  1:1008228/1‑52 (MQ=255)
                                agagGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCaa                                                   >  1:2098748/1‑52 (MQ=255)
                                  agGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGt                              <  1:3034616/71‑1 (MQ=255)
                                  agGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGt                              <  1:797265/71‑1 (MQ=255)
                                  agGAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGt                              <  1:1774483/71‑1 (MQ=255)
                                    gAAGTGATACTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGa                                           <  1:2034478/56‑1 (MQ=255)
                                       gTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGa                         <  1:1754654/71‑1 (MQ=255)
                                         gATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAgg                          >  1:464242/1‑68 (MQ=255)
                                         gATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAgg                          >  1:1524469/1‑68 (MQ=255)
                                         gATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAgg                          >  1:1377421/1‑68 (MQ=255)
                                         gATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAgg                          >  1:2097059/1‑68 (MQ=255)
                                         gATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAgg                          >  1:2075953/1‑68 (MQ=255)
                                            aCTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGc                                       <  1:1029724/52‑1 (MQ=255)
                                            aCTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGc                                       <  1:3828780/52‑1 (MQ=255)
                                            aCTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGc                                       <  1:1898899/52‑1 (MQ=255)
                                            aCTGGGTTTGATCGGCAGTCAGGCGAGAAAACGGCGGCAACACGCCGAAAGc                                       <  1:1054245/52‑1 (MQ=255)
                                             cTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGc                                       <  1:2394248/51‑1 (MQ=255)
                                             cTGGGTTTGATCGGCAGTCAAGCGTGAAACCGGCGGCAACACGCCg                                            <  1:1652683/46‑1 (MQ=255)
                                             cTGGGTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCg                                            <  1:2569633/46‑1 (MQ=255)
                                              tGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATc                                    <  1:3492057/53‑1 (MQ=255)
                                               gggTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAAcc                 >  1:416183/1‑71 (MQ=255)
                                                ggTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCt                >  1:3226000/1‑71 (MQ=255)
                                                ggTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCt                >  1:2725084/1‑71 (MQ=255)
                                                ggTTTGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCt                >  1:4204229/1‑71 (MQ=255)
                                                  tttGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagc                                 <  1:3006845/52‑1 (MQ=255)
                                                  tttGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:3937958/71‑1 (MQ=255)
                                                  tttGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:2967244/71‑1 (MQ=255)
                                                  tttGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:1759179/71‑1 (MQ=255)
                                                  tttGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCGGCAGTCAGGAAGATAACCTTa              <  1:1607480/71‑1 (MQ=255)
                                                  tttGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:4027535/71‑1 (MQ=255)
                                                  tttGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:2657527/71‑1 (MQ=255)
                                                  tttGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:2830531/71‑1 (MQ=255)
                                                  tttGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:1118304/71‑1 (MQ=255)
                                                  tttGATCGGCAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:1015387/71‑1 (MQ=255)
                                                  tttGATCGACAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              <  1:231093/71‑1 (MQ=255)
                                                     gATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGc                                       <  1:719586/43‑1 (MQ=255)
                                                      aTCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagca                                <  1:4241061/49‑1 (MQ=255)
                                                      aTCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagca                                <  1:1783234/49‑1 (MQ=255)
                                                      aTCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATcagca                                <  1:1321340/49‑1 (MQ=255)
                                                           cAGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTgg     <  1:735127/71‑1 (MQ=255)
                                                            aGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTa              >  1:3458700/1‑61 (MQ=255)
                                                            aGTCAAGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGc    >  1:2984594/1‑71 (MQ=255)
                                                              tCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAAc                  <  1:752522/55‑1 (MQ=255)
                                                               cAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGccc  >  1:3177900/1‑70 (MQ=255)
                                                                 |                                                                   
CTCAGTACCGGCCAGTTTGGCGGTGAAGCCAGAGAGGAAGTGATACTGGGTTTGATCGGCAGTCAGGCGAGAAACCGGCGGCAACACGCCGAAAGCATCAGCAGTCAGGAAGATAACCTTAGTCGCGTGGCCC  >  minE/2687094‑2687226

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: