Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,794	T→A	16.7%	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	16.7%	195.4 / 30.8	78	T154S (ACT→TCT) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (40/25);  new base A (8/5);  total (48/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

GCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATC  >  minE/609726‑609863                                                                     |                                                                      gccacggccaTATGCATCTGGCAGACCGGTCAGAAGACCAGATTTACGGCAACGCAGGATGTCCGGAgtg                                                                      >  1:530807/1‑70 (MQ=255)  ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                                    <  1:3791377/71‑1 (MQ=255)   cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                                    <  1:293023/70‑1 (MQ=255)      ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                                <  1:2958091/71‑1 (MQ=255)       gccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                                  >  1:2523725/1‑68 (MQ=255)          ttatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                                    <  1:4005554/62‑1 (MQ=255)          atatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                            <  1:1666750/71‑1 (MQ=255)           tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                           <  1:1281979/71‑1 (MQ=255)             tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                         >  1:2272082/1‑71 (MQ=255)             tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                         >  1:1758809/1‑71 (MQ=255)             tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                         >  1:2561492/1‑71 (MQ=255)             tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                         >  1:1967169/1‑71 (MQ=255)             tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                         >  1:2104948/1‑71 (MQ=255)             gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                           >  1:721075/3‑69 (MQ=255)             gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                           >  1:1757432/3‑69 (MQ=255)             gccATCTGGTAAACAGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                           >  1:1214781/3‑69 (MQ=255)               cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATca                                                             >  1:1686080/1‑64 (MQ=255)                 tCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGcc                                                     >  1:2871227/1‑71 (MQ=255)                 tCTGGCAGACC‑GTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                    >  1:3835619/1‑71 (MQ=255)                  cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                                  >  1:3630091/1‑57 (MQ=255)                   tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                                  >  1:4131320/1‑56 (MQ=255)                   tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                                  >  1:4108083/1‑56 (MQ=255)                   tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                                   >  1:1463254/1‑71 (MQ=255)                   tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                                   >  1:2667979/1‑71 (MQ=255)                      cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGt                                                >  1:3626158/1‑71 (MQ=255)                       aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt                                                >  1:588441/3‑70 (MQ=255)                       aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt                                                >  1:2203540/3‑70 (MQ=255)                       aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                               >  1:2484232/1‑71 (MQ=255)                       aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                               >  1:583371/1‑71 (MQ=255)                        gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                           <  1:482496/58‑1 (MQ=255)                        gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                           <  1:3222016/58‑1 (MQ=255)                        gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtg                                              >  1:2602522/1‑71 (MQ=255)                          ccGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTtatg                                            >  1:916180/1‑71 (MQ=255)                          ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                             >  1:635880/1‑70 (MQ=255)                           cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                             >  1:1938058/1‑69 (MQ=255)                           cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                             >  1:1610704/1‑69 (MQ=255)                              tCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                                                   <  1:1840715/44‑1 (MQ=12)                                 gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                         >  1:2404495/1‑51 (MQ=255)                                 gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                                      <  1:32959/70‑1 (MQ=255)                                  aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                                  >  1:3022398/1‑41 (MQ=255)                                  aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAc                                    >  1:677873/1‑71 (MQ=255)                                   acgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                   <  1:241043/68‑1 (MQ=255)                                   acacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                   <  1:3980117/71‑1 (MQ=255)                                   acacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                   <  1:1709024/71‑1 (MQ=255)                                    cgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                   >  1:952727/3‑70 (MQ=255)                                    cgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                   <  1:1354525/68‑1 (MQ=255)                                      ccAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                           <  1:2034711/44‑1 (MQ=255)                                      ccAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTa                                >  1:2956282/1‑71 (MQ=255)                                       cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAGACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                               >  1:3123276/1‑71 (MQ=255)                                       cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                               >  1:3127472/1‑71 (MQ=255)                                       cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                               >  1:1292000/1‑71 (MQ=255)                                        aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACg                                                       >  1:3724182/1‑46 (MQ=255)                                        aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                              >  1:2109699/1‑71 (MQ=255)                                         gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTgg                                                 <  1:746411/51‑1 (MQ=255)                                          aTTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGttt                                       <  1:3061216/60‑1 (MQ=12)                                          aTTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCa                            <  1:2775210/71‑1 (MQ=255)                                            ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                          >  1:3992948/1‑71 (MQ=255)                                            ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                          >  1:838025/1‑71 (MQ=255)                                               cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg                           <  1:2792215/67‑1 (MQ=255)                                               cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg                           <  1:1061055/67‑1 (MQ=255)                                                 gCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                           <  1:577511/49‑1 (MQ=255)                                                  cAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAt                    <  1:787331/71‑1 (MQ=255)                                                     cGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt                 <  1:180482/71‑1 (MQ=255)                                                      gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                <  1:1141544/71‑1 (MQ=255)                                                      gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                <  1:2791061/71‑1 (MQ=255)                                                      gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt                 >  1:3853816/1‑70 (MQ=255)                                                        aGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                         <  1:3829167/60‑1 (MQ=255)                                                         ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg             <  1:2070941/71‑1 (MQ=255)                                                          gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                         <  1:1436360/58‑1 (MQ=255)                                                           aTGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa            <  1:2790006/70‑1 (MQ=255)                                                            tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat           >  1:3114934/1‑70 (MQ=255)                                                             gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt               >  1:638912/1‑65 (MQ=255)                                                             gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt               >  1:113296/1‑65 (MQ=255)                                                             gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtatttt               >  1:903910/1‑65 (MQ=255)                                                             gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatca         >  1:477762/1‑71 (MQ=255)                                                                cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt               >  1:2178093/1‑62 (MQ=255)                                                                cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcatc       <  1:107495/70‑1 (MQ=255)                                                                cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg      >  1:1080820/1‑71 (MQ=255)                                                                  gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa            >  1:1803670/1‑63 (MQ=255)                                                                    gtgtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATc  <  1:280475/71‑1 (MQ=255)                                                                     tgtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                                  <  1:2589183/38‑1 (MQ=255)                                                                     |                                                                      GCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATC  >  minE/609726‑609863

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