Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,706	G→A	12.5%	D516D (GAC→GAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,706	0	G	A	12.5%	237.0 / 19.3	88	D516D (GAC→GAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (33/44);  new base A (6/5);  total (39/49)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGC  >  minE/608638‑608774                                                                      |                                                                     gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGt                                                                     <  1:882219/70‑1 (MQ=255) gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTc                                                                    >  1:3093960/1‑71 (MQ=255)    gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                                                 >  1:989717/1‑71 (MQ=255)     aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                                                 >  1:1112320/1‑70 (MQ=255)     aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                                                 >  1:3782036/1‑70 (MQ=255)     aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                                                 >  1:827740/1‑70 (MQ=255)         ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAgcg                                                            >  1:2005090/1‑71 (MQ=255)         ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAgcg                                                            <  1:3705642/71‑1 (MQ=255)            gTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCa                                                         <  1:2950475/71‑1 (MQ=255)            gTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCa                                                         <  1:3249682/71‑1 (MQ=255)             tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCa                                                         <  1:1277/70‑1 (MQ=255)             tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                        <  1:609911/71‑1 (MQ=255)             tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                        >  1:3428425/1‑71 (MQ=255)               agcagcAATGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCa                                                      <  1:98677/71‑1 (MQ=255)               agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCa                                                      <  1:2589152/71‑1 (MQ=255)                 cagcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAgcg                                                            <  1:3737942/63‑1 (MQ=255)                 cagcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCaga                                                    <  1:1218174/71‑1 (MQ=255)                  agcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGcc                                                          <  1:895464/64‑1 (MQ=255)                  agcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagag                                                   <  1:3678985/71‑1 (MQ=255)                   gcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCa                                                      >  1:220826/1‑67 (MQ=255)                     aaCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                                                               <  1:2887690/56‑1 (MQ=255)                       ccgACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCt                                              <  1:1898037/69‑1 (MQ=255)                       cGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCt                                              >  1:3421578/1‑71 (MQ=255)                        ggACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCag                                                     >  1:1146727/1‑63 (MQ=255)                        cgACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCtt                                             >  1:150386/2‑71 (MQ=255)                         gACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                                                 >  1:1831780/1‑50 (MQ=255)                          aCAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCg                                           >  1:943171/1‑71 (MQ=255)                          aCAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                           >  1:3372418/1‑71 (MQ=255)                          aCAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                           >  1:3799173/1‑71 (MQ=255)                          aCAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                           >  1:1107471/1‑71 (MQ=255)                           cAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                           >  1:3149032/1‑70 (MQ=255)                           cAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGt                                          <  1:3694058/71‑1 (MQ=255)                            aGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                           >  1:96270/1‑69 (MQ=255)                             gACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCGGAGAGGCTTCGTa                                         <  1:215811/70‑1 (MQ=255)                             gACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTa                                         >  1:3193277/1‑70 (MQ=255)                              aCCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGc                                       <  1:1685785/71‑1 (MQ=255)                                cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACgg                                                                      <  1:2526365/38‑1 (MQ=255)                                cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                                                               <  1:2533402/45‑1 (MQ=255)                                  gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGt                                    >  1:2710663/1‑70 (MQ=255)                                  gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                   >  1:194067/1‑71 (MQ=255)                                  gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                   >  1:1630226/1‑71 (MQ=255)                                   cGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTAc                                  <  1:2629928/71‑1 (MQ=255)                                    gATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTAc                                  <  1:3079912/70‑1 (MQ=255)                                     aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                     <  1:1217645/66‑1 (MQ=255)                                     aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACtt                                <  1:1694886/71‑1 (MQ=255)                                      tACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTg                               >  1:1711142/1‑71 (MQ=255)                                       aCCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                              >  1:411857/1‑71 (MQ=255)                                       aCCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                              >  1:1712925/1‑71 (MQ=255)                                        ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTc                             <  1:3614236/71‑1 (MQ=255)                                         cacaCGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCg                            <  1:256898/71‑1 (MQ=255)                                           catgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGt                           <  1:1649123/67‑1 (MQ=255)                                           cacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                              >  1:1554954/1‑67 (MQ=255)                                             tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa                        >  1:2448498/2‑71 (MQ=255)                                             tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa                        >  1:514557/2‑71 (MQ=255)                                             cGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                   <  1:1576949/60‑1 (MQ=255)                                                 aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATg                      >  1:3881228/1‑69 (MQ=255)                                                 aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTa                    <  1:1979514/71‑1 (MQ=255)                                                 aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTa                    <  1:448813/71‑1 (MQ=255)                                                    gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGt                  <  1:1766988/70‑1 (MQ=255)                                                    gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTg                 <  1:2598374/71‑1 (MQ=255)                                                     tGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                <  1:3700027/71‑1 (MQ=255)                                                     tGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                <  1:693798/71‑1 (MQ=255)                                                       cGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTc                                            >  1:617567/1‑41 (MQ=255)                                                       cGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCTTGCATGTAGTGga               <  1:180471/70‑1 (MQ=255)                                                       cGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga               <  1:3045120/70‑1 (MQ=255)                                                       cGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga               <  1:115108/70‑1 (MQ=255)                                                        ggATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatg             <  1:1051507/71‑1 (MQ=255)                                                        ggATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgat              >  1:1712621/1‑70 (MQ=255)                                                        ggATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgat              >  1:2498510/1‑70 (MQ=255)                                                         gATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga            <  1:1510895/71‑1 (MQ=255)                                                           tA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg          <  1:2111386/71‑1 (MQ=255)                                                            a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt         >  1:1109247/1‑71 (MQ=255)                                                            a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt         >  1:2848019/1‑71 (MQ=255)                                                              acatCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt        <  1:2451649/68‑1 (MQ=255)                                                              aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTc                                            <  1:2784786/35‑1 (MQ=255)                                                              aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt         <  1:284503/70‑1 (MQ=255)                                                                gTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca      <  1:3765323/71‑1 (MQ=255)                                                                atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt        >  1:1375221/2‑69 (MQ=255)                                                                  cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagc    >  1:1683013/1‑70 (MQ=255)                                                                   aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                     <  1:2153764/53‑1 (MQ=255)                                                                   aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                     <  1:3940187/53‑1 (MQ=255)                                                                   aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg   >  1:3690028/1‑71 (MQ=255)                                                                   aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg   <  1:2616968/71‑1 (MQ=255)                                                                   aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg   <  1:583284/71‑1 (MQ=255)                                                                   aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg   >  1:2448796/1‑71 (MQ=255)                                                                   aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg   <  1:3659170/71‑1 (MQ=255)                                                                   aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagt    >  1:1746528/1‑70 (MQ=255)                                                                    cGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCg                            <  1:2605998/45‑1 (MQ=255)                                                                    cGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgc  <  1:1547791/71‑1 (MQ=255)                                                                     ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt        <  1:1271516/64‑1 (MQ=255)                                                                     ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt        <  1:923735/64‑1 (MQ=255)                                                                     ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg   >  1:3745091/1‑69 (MQ=255)                                                                     ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc       >  1:2678580/1‑65 (MQ=255)                                                                      gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGc                         <  1:2052728/46‑1 (MQ=255)                                                                      gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca      <  1:827209/65‑1 (MQ=255)                                                                      |                                                                     GCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGC  >  minE/608638‑608774

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