Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,276	T→C	6.7%	E326E (GAA→GAG) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,276	0	T	C	6.7%	253.2 / 6.4	89	E326E (GAA→GAG) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (53/30);  new base C (4/2);  total (57/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

CAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATT  >  minE/609207‑609342                                                                      |                                                                   cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGtt                                                                   >  1:4221212/1‑71 (MQ=255) cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                                                    >  1:3660695/1‑70 (MQ=255) cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                                                    >  1:2647042/1‑70 (MQ=255)    ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                                                                <  1:1293650/71‑1 (MQ=255)      tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                                              <  1:3889893/71‑1 (MQ=255)      tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                                              <  1:2079290/71‑1 (MQ=255)          aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                                                          >  1:1698968/1‑71 (MQ=255)          aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAt                                                           >  1:624818/1‑70 (MQ=255)           ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                         >  1:4124624/1‑71 (MQ=255)           ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                         >  1:572776/1‑71 (MQ=255)             gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                       <  1:850734/71‑1 (MQ=255)             gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                       <  1:1691527/71‑1 (MQ=255)              tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                       >  1:4205076/1‑70 (MQ=255)                 aGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                                                    >  1:4514107/1‑70 (MQ=255)                     ccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTGCGGAGTACg                                                >  1:3484677/1‑70 (MQ=255)                       aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                             >  1:4629968/1‑71 (MQ=255)                        tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                            <  1:3820930/71‑1 (MQ=255)                        tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                            <  1:1712763/71‑1 (MQ=255)                          ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                                           <  1:821195/70‑1 (MQ=255)                          ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                                           <  1:3361070/70‑1 (MQ=255)                            accGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                                       >  1:378396/1‑56 (MQ=255)                             ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                         >  1:5616054/1‑53 (MQ=255)                              cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                       <  1:1657633/70‑1 (MQ=255)                               gATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATc                                                               >  1:503729/1‑45 (MQ=255)                                aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTc                                                                  <  1:5490046/41‑1 (MQ=255)                                aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                    >  1:4952182/1‑71 (MQ=255)                                 tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                                         <  1:4114301/49‑1 (MQ=255)                                 tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                      <  1:2988351/68‑1 (MQ=255)                                  aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc                                        >  1:1460356/1‑65 (MQ=255)                                  aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                   >  1:5282517/1‑70 (MQ=255)                                  aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                    >  1:4756260/1‑69 (MQ=255)                                  aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                  >  1:399609/1‑71 (MQ=255)                                   gATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                  <  1:4702818/70‑1 (MQ=255)                                   acttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCa                                  >  1:5649674/3‑70 (MQ=255)                                       cGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACg                              >  1:4007674/1‑70 (MQ=255)                                            gCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGtt                         >  1:1803539/1‑70 (MQ=255)                                             cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                       <  1:3545860/71‑1 (MQ=255)                                             cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                       <  1:5186223/71‑1 (MQ=255)                                              ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                      >  1:3512965/1‑71 (MQ=255)                                               cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                     <  1:3130415/71‑1 (MQ=255)                                               cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                     <  1:2266684/71‑1 (MQ=255)                                               cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                       <  1:42153/53‑1 (MQ=255)                                               cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                       <  1:3613448/53‑1 (MQ=255)                                               cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                       <  1:4897229/53‑1 (MQ=255)                                                aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                     >  1:3455217/1‑70 (MQ=255)                                                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                     >  1:4891136/1‑70 (MQ=255)                                                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                     >  1:123839/1‑70 (MQ=255)                                                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                     >  1:3794498/1‑70 (MQ=255)                                                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                     >  1:3589419/1‑70 (MQ=255)                                                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                     >  1:2193704/1‑70 (MQ=255)                                                aGATCGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                     >  1:4897799/1‑70 (MQ=255)                                                 gATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                   >  1:3407185/1‑71 (MQ=255)                                                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTTCGGAGTACGCAGGAAGCGa                                      >  1:190192/1‑52 (MQ=255)                                                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGc                                               <  1:2008814/43‑1 (MQ=255)                                                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                     >  1:5165173/1‑53 (MQ=255)                                                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                     >  1:4145257/1‑53 (MQ=255)                                                 gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                     >  1:2798119/1‑53 (MQ=255)                                                  aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                     >  1:5447180/1‑36 (MQ=255)                                                  aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                     >  1:440374/1‑36 (MQ=255)                                                  aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                  <  1:3607582/71‑1 (MQ=255)                                                   ttgggTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                   >  1:1024587/3‑69 (MQ=255)                                                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCa                                              >  1:352507/1‑42 (MQ=255)                                                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                      >  1:1816350/1‑50 (MQ=255)                                                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                >  1:3601350/1‑56 (MQ=255)                                                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                >  1:23699/1‑56 (MQ=255)                                                   tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                 <  1:3888221/71‑1 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                   >  1:5637492/1‑52 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                   >  1:5177806/1‑52 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                   >  1:3362506/1‑52 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                   >  1:2447193/1‑52 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                        >  1:722803/1‑63 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                      >  1:5626028/1‑65 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                 <  1:2468755/70‑1 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                 <  1:5527717/70‑1 (MQ=255)                                                    cGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCg                                       >  1:5453611/1‑48 (MQ=255)                                                      ggTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                  <  1:2814762/51‑1 (MQ=255)                                                         cGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtg            >  1:4023666/1‑70 (MQ=255)                                                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                       >  1:3567643/1‑42 (MQ=255)                                                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                  >  1:4516847/1‑63 (MQ=255)                                                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt          >  1:3271539/1‑71 (MQ=255)                                                          gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt          >  1:5303577/1‑71 (MQ=255)                                                            cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa        <  1:5031927/71‑1 (MQ=255)                                                            cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGt                          <  1:228596/53‑1 (MQ=255)                                                            cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa        <  1:278814/71‑1 (MQ=255)                                                             agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAg               >  1:1972603/1‑63 (MQ=255)                                                              gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                <  1:4348556/45‑1 (MQ=255)                                                              gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat       >  1:878281/1‑69 (MQ=255)                                                              gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat       >  1:2931471/1‑69 (MQ=255)                                                                gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa    <  1:1885619/71‑1 (MQ=255)                                                                  aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt  >  1:1719359/1‑71 (MQ=255)                                                                  aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt  >  1:4089161/1‑71 (MQ=255)                                                                   cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt  >  1:2999037/1‑70 (MQ=255)                                                                      ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                        <  1:2220146/45‑1 (MQ=255)                                                                      ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                        <  1:4923374/45‑1 (MQ=255)                                                                      ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                        <  1:64744/45‑1 (MQ=255)                                                                      ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                        <  1:4537443/45‑1 (MQ=255)                                                                      |                                                                   CAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATT  >  minE/609207‑609342

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