Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,319	T→G	8.2%	A645A (GCA→GCC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,319	0	T	G	8.2%	146.5 / 10.1	61	A645A (GCA→GCC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base G (2/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGT  >  minE/608251‑608389                                                                     |                                                                       aaGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGc                                                                      <  1:2397246/71‑1 (MQ=255)  aGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCa                                                                     >  1:66239/1‑71 (MQ=255)       tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCttt                                                                >  1:1331584/1‑71 (MQ=255)              tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                                           <  1:433874/69‑1 (MQ=255)               aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                        >  1:3711575/1‑71 (MQ=255)               aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                        >  1:2276328/1‑71 (MQ=255)               aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                        >  1:1997276/1‑71 (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: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:3206440/1‑71 (MQ=255)                           gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                                            <  1:3083121/71‑1 (MQ=255)                            tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                                            <  1:1532937/54‑1 (MQ=255)                              aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                         >  1:1218867/1‑71 (MQ=255)                              aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                         >  1:3070467/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:249141/1‑70 (MQ=255)                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                            >  1:2777947/1‑70 (MQ=255)                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                            >  1:816383/1‑70 (MQ=255)                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                            >  1:1308062/1‑70 (MQ=255)                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                            >  1:2686875/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:2555676/37‑1 (MQ=255)                                                        tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                 <  1:735118/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:2643037/70‑1 (MQ=255)                                                          agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg              <  1:3095710/70‑1 (MQ=255)                                                          agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGc             <  1:1417372/71‑1 (MQ=255)                                                             cagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacg          <  1:2164731/70‑1 (MQ=255)                                                              agGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacg          <  1:1962260/70‑1 (MQ=255)                                                                gCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAg                                      <  1:813800/40‑1 (MQ=255)                                                                    cTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg   >  1:2325586/1‑71 (MQ=255)                                                                    cTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg   <  1:360988/71‑1 (MQ=255)                                                                     tGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGt  >  1:1337665/1‑71 (MQ=255)                                                                     tGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGt  >  1:4043511/1‑71 (MQ=255)                                                                     |                                                                       AAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGT  >  minE/608251‑608389

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