Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,222	G→C	23.9%	L344L (CTC→CTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,222	0	G	C	23.9%	177.2 / 50.2	92	L344L (CTC→CTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base C (14/8);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCG  >  minE/609154‑609289                                                                     |                                                                    gaCGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAg                                                                   <  1:1688737/71‑1 (MQ=255)    ggccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc                                                                >  1:88760/4‑71 (MQ=255)    ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                                <  1:831753/71‑1 (MQ=255)    ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                                <  1:4683451/71‑1 (MQ=255)    ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                                >  1:4277363/1‑71 (MQ=255)     gccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc                                                                >  1:5585024/3‑70 (MQ=255)     gccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc                                                                >  1:1969608/3‑70 (MQ=255)       cccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCAt                                                             >  1:2876034/1‑71 (MQ=255)        ccATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGa                                                                    >  1:1727957/1‑63 (MQ=255)        ccATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGa                                                                    >  1:3538263/1‑63 (MQ=255)         cATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac                                                           <  1:1498445/71‑1 (MQ=255)         cATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCAtcc                                                           <  1:28585/71‑3 (MQ=255)          aTGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac                                                           >  1:3725049/1‑70 (MQ=255)           tGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATacca                                                         >  1:1355781/1‑71 (MQ=255)             gtgtGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATcccgcc                                                       >  1:3131886/1‑71 (MQ=255)             gtgtGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATcccgc                                                        >  1:4181715/1‑68 (MQ=255)             gtgtACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATaccac                                                        >  1:1664215/1‑70 (MQ=255)                tACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGAt                                                    <  1:2213120/71‑1 (MQ=255)                  cAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                     >  1:1558318/1‑68 (MQ=255)                  cAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATa                                                   <  1:4373535/70‑1 (MQ=255)                    gggTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATAACACCGATAGAt                                                <  1:4219930/71‑1 (MQ=255)                         ttCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                                <  1:1232033/50‑1 (MQ=255)                         ttCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                                <  1:2075826/50‑1 (MQ=255)                         ttCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGt                                           >  1:146097/1‑71 (MQ=255)                           cAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac                                                           >  1:1681998/1‑53 (MQ=255)                            aGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                         <  1:1815458/70‑1 (MQ=255)                            aGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGc                                        <  1:3223435/71‑1 (MQ=255)                                aaCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCAt                                                             >  1:180679/1‑46 (MQ=255)                                aaCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAg                                    <  1:5120824/71‑1 (MQ=255)                                 aCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGa                                   <  1:1618778/71‑1 (MQ=255)                                  cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCg                                                      >  1:254371/1‑51 (MQ=255)                                  cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                     >  1:3229867/1‑52 (MQ=255)                                  cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                     >  1:1655566/1‑52 (MQ=255)                                  cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                       >  1:2858751/1‑66 (MQ=255)                                  cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGAt                                  <  1:1468218/71‑1 (MQ=255)                                   ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAt                                                    >  1:2928662/1‑52 (MQ=16)                                   ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAt                                                    >  1:3848782/1‑52 (MQ=16)                                   ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCAt                                                             >  1:3627170/1‑43 (MQ=255)                                   ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                            >  1:2470007/1‑44 (MQ=255)                                   ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                       >  1:4980736/1‑65 (MQ=255)                                      aGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                               <  1:2476638/70‑1 (MQ=255)                                       ggagTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGt                             <  1:2964678/68‑1 (MQ=255)                                          gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                  >  1:1750265/1‑45 (MQ=16)                                          gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                  >  1:5579429/1‑45 (MQ=16)                                          gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                  >  1:417144/1‑45 (MQ=16)                                          gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                  >  1:2007762/1‑45 (MQ=16)                                          gTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCgg                                            >  1:4654558/1‑53 (MQ=255)                                           tttttGGTCACCAGCGTACGACCGTCCAGCCCCATcccgc                                                        <  1:2577542/40‑3 (MQ=12)                                           tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtc                                                   <  1:3875433/45‑2 (MQ=12)                                           tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTggg                              >  1:4588029/1‑66 (MQ=255)                                           tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                            <  1:3431981/36‑1 (MQ=255)                                           tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                            <  1:1202690/36‑1 (MQ=255)                                           tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                            <  1:3161205/36‑1 (MQ=255)                                           tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                            <  1:1174576/36‑1 (MQ=255)                                           tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                            <  1:1735616/36‑1 (MQ=255)                                           tttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                          <  1:151642/70‑1 (MQ=255)                                            ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                          >  1:1515688/1‑69 (MQ=255)                                            ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                          >  1:2138134/1‑69 (MQ=255)                                             tttGGTAACCAGGGTACGACCGTCGAGGCCCATacca                                                         >  1:1993855/1‑37 (MQ=255)                                             tttGGTAACCAGGGTACGACCGTCGAGGCCCATacc                                                          >  1:266622/1‑36 (MQ=255)                                             tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCg                                >  1:2000489/1‑62 (MQ=255)                                             tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                         >  1:2712682/1‑69 (MQ=255)                                             tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                         >  1:2474570/1‑69 (MQ=255)                                              ttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                       <  1:1269378/70‑3 (MQ=255)                                                ggTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCagag                     >  1:1969865/1‑70 (MQ=255)                                                 gtcaCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGaa                   <  1:4104932/68‑1 (MQ=255)                                                 gTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                    <  1:2512797/70‑1 (MQ=255)                                                    aCCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                >  1:845085/1‑71 (MQ=255)                                                    aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGcc                         <  1:765651/62‑1 (MQ=255)                                                    aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAg                <  1:3508966/71‑1 (MQ=255)                                                      cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGtt              >  1:4221212/1‑71 (MQ=255)                                                      cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt               >  1:2647042/1‑70 (MQ=255)                                                      cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt               >  1:3660695/1‑70 (MQ=255)                                                      cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                       >  1:2796722/1‑60 (MQ=255)                                                         ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt           <  1:1293650/71‑1 (MQ=255)                                                           tACGACCGTCGAGGCCCATACCACCGATAGATTCgg                                            >  1:4678409/1‑36 (MQ=255)                                                           tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg         <  1:3889893/71‑1 (MQ=255)                                                           tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg         <  1:2079290/71‑1 (MQ=255)                                                             cGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                         <  1:3613668/37‑1 (MQ=255)                                                             cGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                         <  1:3607041/37‑1 (MQ=255)                                                               aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt     >  1:1698968/1‑71 (MQ=255)                                                               aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAt      >  1:624818/1‑70 (MQ=255)                                                                ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc    >  1:572776/1‑71 (MQ=255)                                                                ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc    >  1:4124624/1‑71 (MQ=255)                                                                  gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                               <  1:2518763/42‑1 (MQ=255)                                                                  gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg  <  1:1691527/71‑1 (MQ=255)                                                                  gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg  <  1:850734/71‑1 (MQ=255)                                                                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                 <  1:4247526/55‑1 (MQ=255)                                                                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                 <  1:4986670/55‑1 (MQ=255)                                                                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                 <  1:3899602/55‑1 (MQ=255)                                                                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                 <  1:3307660/55‑1 (MQ=255)                                                                   tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg  >  1:4205076/1‑70 (MQ=255)                                                                     |                                                                    GACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCG  >  minE/609154‑609289

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