Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,298	A→C	9.1%	V652V (GTT→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,298	0	A	C	9.1%	157.6 / 13.4	66	V652V (GTT→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base C (3/3);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACG  >  minE/608232‑608368                                                                   |                                                                       tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                    <  1:542023/71‑1 (MQ=255) tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                    <  1:491343/71‑1 (MQ=255)   cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                    >  1:2842881/1‑69 (MQ=255)       cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTc                                                               <  1:778552/70‑1 (MQ=255)        ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                             <  1:787453/71‑1 (MQ=255)        ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                             <  1:3146303/71‑1 (MQ=255)        ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                             <  1:1559083/71‑1 (MQ=255)        ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                              <  1:1496844/70‑1 (MQ=255)         gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                            <  1:3485917/71‑1 (MQ=255)         gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                            <  1:2443807/71‑1 (MQ=255)           aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                            >  1:2310104/1‑67 (MQ=255)           aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                          <  1:3595411/71‑1 (MQ=255)           aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                          <  1:1848774/71‑1 (MQ=255)            cGATAGAGAAGGTGTAGGAGATACCCTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagagag                                                         >  1:3732737/1‑71 (MQ=255)            cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagagag                                                         >  1:2751276/1‑71 (MQ=255)              aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                            >  1:3789780/1‑66 (MQ=255)              aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                       >  1:1160041/1‑71 (MQ=255)                   gaAGGAGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                     >  1:116959/1‑68 (MQ=255)                    aaGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGc                                                 <  1:2397246/71‑1 (MQ=255)                     aGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCa                                                >  1:66239/1‑71 (MQ=255)                        tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                                    <  1:1991040/48‑1 (MQ=255)                          tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCttt                                           >  1:1331584/1‑71 (MQ=255)                             gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTTGCAACGGa                                                                   <  1:2092376/44‑1 (MQ=255)                               gaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTc                                                               <  1:3062461/46‑1 (MQ=255)                                 tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                            <  1:3026541/47‑1 (MQ=255)                                 tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                      <  1:433874/69‑1 (MQ=255)                                  aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                              <  1:1565750/44‑1 (MQ=255)                                  aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                              <  1:2312224/44‑1 (MQ=255)                                  aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                   >  1:3711575/1‑71 (MQ=255)                                  aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                   >  1:1997276/1‑71 (MQ=255)                                  aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                   >  1:2276328/1‑71 (MQ=255)                                     aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagag                                                           <  1:3950012/44‑1 (MQ=255)                                      tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCt                                                      <  1:612356/48‑1 (MQ=255)                                      tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                         <  1:931244/61‑1 (MQ=255)                                      tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                         <  1:57827/61‑1 (MQ=255)                                      tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                       <  1:2276715/63‑1 (MQ=255)                                      tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                               >  1:2824860/1‑71 (MQ=255)                                        tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                             >  1:86378/1‑71 (MQ=255)                                          tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                            >  1:2871084/1‑70 (MQ=255)                                          tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                            >  1:162429/1‑70 (MQ=255)                                          tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                            >  1:2274216/1‑70 (MQ=255)                                          tAGCGTAAGCAAACCGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                            >  1:3170732/1‑70 (MQ=255)                                            gCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                  <  1:2989302/46‑1 (MQ=255)                                              gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                       <  1:3083121/71‑1 (MQ=255)                                              gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                       >  1:3206440/1‑71 (MQ=255)                                               tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                       <  1:1532937/54‑1 (MQ=255)                                                 aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                    >  1:3070467/1‑71 (MQ=255)                                                 aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                    >  1:1218867/1‑71 (MQ=255)                                                 aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                    >  1:1485182/1‑71 (MQ=255)                                                     aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                   <  1:51290/68‑1 (MQ=255)                                                      aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                     >  1:1588358/1‑49 (MQ=255)                                                       cGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAc              >  1:2351881/1‑71 (MQ=255)                                                        ggCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc             <  1:1757453/71‑1 (MQ=255)                                                          cAGTTTCGCCACCGAGGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCg           >  1:244307/1‑71 (MQ=255)                                                            gTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCggc         >  1:442481/1‑70 (MQ=255)                                                             tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                 <  1:3423367/62‑1 (MQ=255)                                                             tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                 <  1:3418642/62‑1 (MQ=255)                                                               tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                      >  1:2472862/1‑55 (MQ=255)                                                               tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc       >  1:1308062/1‑70 (MQ=255)                                                               tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc       >  1:2777947/1‑70 (MQ=255)                                                               tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc       >  1:2686875/1‑70 (MQ=255)                                                               tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc       >  1:816383/1‑70 (MQ=255)                                                               tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc       >  1:249141/1‑70 (MQ=255)                                                                aGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                   <  1:3357585/57‑1 (MQ=255)                                                                 gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGa     <  1:2898274/70‑1 (MQ=255)                                                                  cAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCGCGACCGTGCATCgg                      <  1:2989599/52‑1 (MQ=255)                                                                   aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                      >  1:1076369/1‑51 (MQ=255)                                                                   aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg  >  1:3648859/1‑71 (MQ=255)                                                                   aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg  <  1:1243177/71‑1 (MQ=255)                                                                   |                                                                       TGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACG  >  minE/608232‑608368

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