Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,094	C→G	8.9%	M720I (ATG→ATC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,094	0	C	G	8.9%	196.4 / 8.2	67	M720I (ATG→ATC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base G (3/3);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

AGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCT  >  minE/608026‑608163                                                                     |                                                                      aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                                                     >  1:2956635/1‑71 (MQ=255)     aGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATcgcg                                                                 <  1:2970336/71‑1 (MQ=255)         gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTc                                                               >  1:3850627/1‑69 (MQ=255)            aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                            >  1:255551/1‑69 (MQ=255)               tACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc                                                                       >  1:2191725/1‑55 (MQ=255)               tACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc                                                                       >  1:2499353/1‑55 (MQ=255)                 cTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAtt                                                     <  1:2941602/71‑1 (MQ=255)                 cTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAtt                                                     >  1:3755243/1‑71 (MQ=255)                 cTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAtt                                                     >  1:833715/1‑71 (MQ=255)                   gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAtt                                                     >  1:789006/1‑69 (MQ=255)                   gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAtt                                                     >  1:1781679/1‑69 (MQ=255)                        gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                               <  1:4059636/70‑1 (MQ=255)                        gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                               <  1:1554197/70‑1 (MQ=255)                        gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                               <  1:2782407/70‑1 (MQ=255)                        gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                               <  1:3759950/70‑1 (MQ=255)                        gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                              >  1:2263107/1‑71 (MQ=255)                         ccAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGtt                                             >  1:678018/1‑71 (MQ=255)                              tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                            >  1:2113819/1‑51 (MQ=255)                              tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                 >  1:2117928/1‑62 (MQ=255)                              tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatca                                        >  1:2134546/1‑71 (MQ=255)                                cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACg                                                <  1:3307487/61‑1 (MQ=255)                                cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                       <  1:954018/70‑1 (MQ=255)                                  cGGATGGTCAGCT‑‑GGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                              >  1:3782096/1‑59 (MQ=38)                                   ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTa                                   <  1:2867888/71‑1 (MQ=255)                                   ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTa                                   <  1:2862875/71‑1 (MQ=255)                                    gATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttga                                  <  1:2025127/71‑3 (MQ=255)                                    gATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa                                  <  1:3732417/71‑1 (MQ=255)                                        gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                         >  1:2968005/1‑60 (MQ=255)                                        gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                              <  1:1895377/71‑1 (MQ=255)                                         tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                             <  1:2209672/71‑1 (MQ=255)                                          cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                              >  1:341219/1‑69 (MQ=255)                                           aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                               >  1:873134/1‑51 (MQ=255)                                           aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                            <  1:428259/70‑1 (MQ=255)                                               taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAttt                                                    >  1:1509493/3‑42 (MQ=12)                                               taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                       >  1:2042885/3‑70 (MQ=255)                                                cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                          <  1:3068170/51‑1 (MQ=255)                                                cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACg                                <  1:1592085/61‑1 (MQ=255)                                                cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                      <  1:2993827/71‑1 (MQ=255)                                                aggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgtt                      >  1:1240154/2‑69 (MQ=255)                                                 ggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTg                     <  1:808871/71‑1 (MQ=12)                                                  gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGtt                                             <  1:3237633/46‑1 (MQ=255)                                                  gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                    >  1:87146/1‑71 (MQ=255)                                                  gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGAtt                                             <  1:2344523/46‑1 (MQ=12)                                                   atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                              <  1:2751272/60‑1 (MQ=255)                                                   atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGac                   >  1:2458871/1‑71 (MQ=255)                                                      tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                          <  1:2823193/45‑1 (MQ=255)                                                        tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                            >  1:1845417/1‑57 (MQ=255)                                                        tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                      >  1:3168975/1‑63 (MQ=255)                                                        tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc               >  1:747842/1‑70 (MQ=255)                                                           ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc            <  1:1257222/70‑1 (MQ=255)                                                           ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc            <  1:2698739/70‑1 (MQ=255)                                                            cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                        <  1:2694289/57‑1 (MQ=255)                                                            cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacc                  <  1:1602284/63‑1 (MQ=255)                                                              ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg        <  1:662534/71‑1 (MQ=255)                                                              ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg        <  1:2205812/71‑1 (MQ=255)                                                              ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg        <  1:757415/71‑1 (MQ=255)                                                               gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                              <  1:3973057/48‑1 (MQ=255)                                                               gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                              <  1:178005/48‑1 (MQ=255)                                                               gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg        >  1:1728911/1‑70 (MQ=255)                                                                ttttCCATCGCGTCGAGCATCATTTCACGGTTCATCACGTTAACGtt                              <  1:3039242/47‑1 (MQ=255)                                                                ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                             <  1:2731694/48‑1 (MQ=255)                                                                 tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                       <  1:1306678/37‑1 (MQ=255)                                                                 tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                       <  1:2931425/37‑1 (MQ=255)                                                                 tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                       <  1:2926582/37‑1 (MQ=255)                                                                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg    >  1:319676/1‑71 (MQ=255)                                                                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg    >  1:1585647/1‑71 (MQ=255)                                                                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat     >  1:3293199/1‑70 (MQ=255)                                                                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg        <  1:1791464/67‑1 (MQ=255)                                                                   tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc   >  1:994847/1‑71 (MQ=255)                                                                    ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt  >  1:1545759/1‑71 (MQ=255)                                                                     |                                                                      AGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCT  >  minE/608026‑608163

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑