Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,309	G→A	9.1%	L649L (CTG→TTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,309	0	G	A	9.1%	166.2 / 6.1	55	L649L (CTG→TTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (30/20);  new base A (2/3);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGC  >  minE/608239‑608378                                                                       |                                                                      ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                       <  1:3146303/71‑1 (MQ=255) ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                       <  1:787453/71‑1 (MQ=255) ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                                       <  1:1559083/71‑1 (MQ=255)  gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                      <  1:2443807/71‑1 (MQ=255)  gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                      <  1:3485917/71‑1 (MQ=255)    aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                                    <  1:1848774/71‑1 (MQ=255)    aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                                    <  1:3595411/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)                   tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCttt                                                     >  1:1331584/1‑71 (MQ=255)                          tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                                      <  1:3026541/47‑1 (MQ=255)                          tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                                <  1:433874/69‑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:57827/61‑1 (MQ=255)                               tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                                   <  1:931244/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:162429/1‑70 (MQ=255)                                   tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                      >  1:2274216/1‑70 (MQ=255)                                   tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                      >  1:2871084/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:816383/1‑70 (MQ=255)                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                 >  1:249141/1‑70 (MQ=255)                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                 >  1:2777947/1‑70 (MQ=255)                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                 >  1:2686875/1‑70 (MQ=255)                                                        tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                 >  1:1308062/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)                                                              cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACTAcc                                         >  1:3876436/1‑40 (MQ=255)                                                                gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc          <  1:18356/69‑1 (MQ=255)                                                                 aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcg       >  1:3561611/1‑71 (MQ=255)                                                                   gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                              <  1:2751199/46‑1 (MQ=255)                                                                   gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGcc     <  1:3416744/71‑1 (MQ=255)                                                                    tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                      <  1:735118/37‑1 (MQ=255)                                                                    tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                      <  1:2555676/37‑1 (MQ=255)                                                                    tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAg                           <  1:2509218/48‑1 (MQ=255)                                                                     cAGAGAGGCTACTGCACCTTTCTGGTCACGAACGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg   >  1:867294/1‑71 (MQ=255)                                                                      agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa              >  1:3127725/1‑59 (MQ=255)                                                                      agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg   <  1:3095710/70‑1 (MQ=255)                                                                      agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg   <  1:2643037/70‑1 (MQ=255)                                                                      agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGc  <  1:1417372/71‑1 (MQ=255)                                                                       |                                                                      GGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGC  >  minE/608239‑608378

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