Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,794	T→A	17.3%	T154S (ACT→TCT) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,794	0	T	A	17.3%	233.9 / 33.5	97	T154S (ACT→TCT) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (38/42);  new base A (13/4);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 3.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCC  >  minE/609727‑609864                                                                    |                                                                       ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgtgt                                                                      >  1:2180042/1‑70 (MQ=255) ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgt                                                                        >  1:2812195/1‑68 (MQ=255)  cacggccaTATGCATCTGTCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                                    <  1:5535/71‑1 (MQ=255)  cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                                                    <  1:3006292/71‑1 (MQ=255)   acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                                   >  1:1205796/1‑71 (MQ=255)    cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                                     <  1:1759126/68‑1 (MQ=255)    cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                                     >  1:2006983/1‑68 (MQ=255)    cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                                   >  1:3015164/1‑70 (MQ=255)    cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                                   >  1:2683014/1‑70 (MQ=255)    cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                                   >  1:2393813/1‑70 (MQ=255)    cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                                   >  1:2389073/1‑70 (MQ=255)     tgccaGATGCATATGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                                 >  1:2894211/2‑71 (MQ=255)     ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                                 >  1:96175/1‑71 (MQ=255)     ggccaTATGCATCTGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                                 <  1:1154638/71‑1 (MQ=255)       ccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAcat                                                                >  1:3254109/1‑68 (MQ=255)            tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                                   >  1:89121/1‑62 (MQ=255)            tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                          <  1:3863764/71‑1 (MQ=255)            tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                          <  1:3713503/71‑1 (MQ=255)            tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                          <  1:1557305/71‑1 (MQ=255)              cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACg                                                        <  1:3294142/71‑1 (MQ=255)              cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                                    >  1:3389602/1‑59 (MQ=255)              cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                                    >  1:2722572/1‑59 (MQ=255)               aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATcaa                                                             <  1:2573065/65‑3 (MQ=255)               aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGc                                                       >  1:934271/1‑71 (MQ=255)               aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                             >  1:622401/1‑65 (MQ=255)               aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                             >  1:3472631/1‑65 (MQ=255)                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                                   >  1:75521/1‑57 (MQ=255)                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                                   >  1:430826/1‑57 (MQ=255)                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                     >  1:1200076/1‑71 (MQ=255)                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                     >  1:683712/1‑71 (MQ=255)                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                     >  1:3124160/1‑71 (MQ=255)                 cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                     >  1:3202735/1‑71 (MQ=255)                   ggCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                             <  1:3232434/61‑1 (MQ=255)                    gCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                             <  1:2694772/60‑1 (MQ=255)                      aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                                     <  1:1952353/50‑1 (MQ=255)                       aaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGtta                                               >  1:3413941/2‑70 (MQ=255)                        aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                            <  1:540759/57‑1 (MQ=255)                        aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                              >  1:1962760/1‑71 (MQ=255)                        aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                              <  1:794064/71‑1 (MQ=255)                         ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                                                                  >  1:3332771/1‑50 (MQ=255)                         ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg                                             <  1:1080095/71‑1 (MQ=255)                          cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                            <  1:472274/71‑1 (MQ=255)                           ggTCAGAACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCACTGGTtatggg                                           >  1:2884392/1‑69 (MQ=255)                            gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                            >  1:2894459/1‑69 (MQ=255)                             tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGGCGAACACGCCCTGGTTGTGAGtt                                         <  1:1754112/71‑1 (MQ=255)                              cAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                        >  1:10591/1‑71 (MQ=255)                               aGAACTCCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                                       <  1:2585645/71‑1 (MQ=255)                               aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTTTAAACGTCGAACACGCCCTGGGTGTGAGtttt                                       <  1:3339867/71‑1 (MQ=255)                               aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                              <  1:3287562/48‑1 (MQ=255)                               aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                        <  1:2490538/70‑1 (MQ=255)                               aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                        <  1:858471/70‑1 (MQ=255)                               aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                        <  1:186373/70‑1 (MQ=255)                                 aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                      <  1:816445/70‑1 (MQ=255)                                  acgccAGATTTACGGCATCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                    >  1:1815777/4‑71 (MQ=255)                                   cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                                   >  1:1129451/1‑71 (MQ=255)                                   cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                                   >  1:437384/1‑71 (MQ=255)                                    gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGttt                                        <  1:1631026/64‑1 (MQ=255)                                    gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                    >  1:2061190/2‑69 (MQ=255)                                      cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                                >  1:1531336/1‑71 (MQ=255)                                       aGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAc                                                                  >  1:3579797/1‑36 (MQ=12)                                       aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                                >  1:1452631/1‑54 (MQ=255)                                        gATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                    >  1:1326924/1‑65 (MQ=255)                                        gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACg                                                        >  1:1028313/1‑45 (MQ=255)                                        gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGt                                  >  1:1099139/1‑67 (MQ=255)                                        gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                               >  1:593872/1‑70 (MQ=255)                                        gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                              <  1:1668607/71‑1 (MQ=255)                                           ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                                   >  1:2056657/1‑47 (MQ=255)                                           ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                           >  1:3031529/1‑71 (MQ=255)                                           ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                           >  1:1708068/1‑71 (MQ=255)                                            tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                              <  1:627472/35‑1 (MQ=255)                                            tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                           <  1:2080340/70‑1 (MQ=255)                                             aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                          <  1:3194530/70‑1 (MQ=255)                                             aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                          <  1:2314118/70‑1 (MQ=255)                                              cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGa                         <  1:1607697/70‑1 (MQ=255)                                                  aaCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtt                    >  1:1050678/1‑71 (MQ=255)                                                   aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                     >  1:649819/1‑37 (MQ=255)                                                   aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                     >  1:1502410/1‑37 (MQ=255)                                                   aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGa                      >  1:3117557/1‑68 (MQ=255)                                                     gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAg                                           <  1:82200/45‑1 (MQ=255)                                                     gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                 <  1:1930890/71‑1 (MQ=255)                                                      cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                >  1:1831887/1‑71 (MQ=255)                                                       aGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttttt               >  1:1797344/1‑71 (MQ=255)                                                        ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg              <  1:1372368/71‑1 (MQ=255)                                                         gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa             <  1:2937464/71‑1 (MQ=255)                                                         gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa             <  1:2008984/71‑1 (MQ=255)                                                         gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa             <  1:2074140/71‑1 (MQ=255)                                                          aTGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat            >  1:3690858/1‑71 (MQ=255)                                                           tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatc           <  1:1569361/71‑1 (MQ=255)                                                           tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat            >  1:2696638/1‑70 (MQ=255)                                                              ccGGAGTGTAAACGTCTAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                           <  1:1123984/52‑1 (MQ=255)                                                               cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcatc        <  1:3007321/70‑1 (MQ=255)                                                               cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcatc        <  1:392250/70‑1 (MQ=255)                                                               cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg       >  1:622611/1‑71 (MQ=255)                                                               cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg       <  1:418031/71‑1 (MQ=255)                                                                ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg      <  1:3602155/71‑1 (MQ=255)                                                                ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg      <  1:2932642/71‑1 (MQ=255)                                                                ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg      <  1:351944/71‑1 (MQ=255)                                                                 gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg                            <  1:1069834/48‑1 (MQ=255)                                                                 gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGa     >  1:985011/1‑71 (MQ=255)                                                                    tgtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATcc  <  1:2847058/71‑1 (MQ=255)                                                                    |                                                                       CCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCC  >  minE/609727‑609864

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: ‑