Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,310	A→T	9.1%	S648S (TCT→TCA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,310	0	A	T	9.1%	149.6 / 7.2	55	S648S (TCT→TCA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (30/20);  new base T (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 = 7.92e-01

GAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGC  >  minE/608240‑608378                                                                       |                                                                     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:1218867/1‑71 (MQ=255)                                         aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                              >  1:1485182/1‑71 (MQ=255)                                         aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                              >  1:3070467/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:3418642/62‑1 (MQ=255)                                                     tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                           <  1:3423367/62‑1 (MQ=255)                                                       tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                                >  1:2472862/1‑55 (MQ=255)                                                       tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                 >  1:2777947/1‑70 (MQ=255)                                                       tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                 >  1:1308062/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)                                                             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)                                                                       |                                                                     GAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGC  >  minE/608240‑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: ‑