Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,904	G→C	15.0%	G450G (GGC→GGG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,904	0	G	C	15.0%	198.6 / 24.6	80	G450G (GGC→GGG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base C (6/6);  total (38/42)
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.82e-01

ATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAC  >  minE/608835‑608974                                                                      |                                                                       atAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgc                                                                       >  1:1091706/1‑71 (MQ=255)      tCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                   <  1:1676283/70‑1 (MQ=255)       cAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGt                                                                 >  1:3295514/1‑71 (MQ=255)       cAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGt                                                                 >  1:3908735/1‑71 (MQ=255)        aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                   <  1:2494603/68‑1 (MQ=255)        aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                                <  1:368580/71‑1 (MQ=255)        aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                                <  1:1245158/71‑1 (MQ=255)         ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                               <  1:3280851/71‑1 (MQ=255)         ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                               <  1:1808459/71‑1 (MQ=255)         ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                               <  1:355522/71‑1 (MQ=255)          gACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                   >  1:2014325/1‑66 (MQ=255)           aCATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                             >  1:3088825/1‑71 (MQ=255)            cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                            >  1:3686774/1‑71 (MQ=255)              tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                               <  1:1961296/66‑1 (MQ=255)                gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                         >  1:3643924/1‑70 (MQ=255)                gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                         <  1:2592926/70‑1 (MQ=255)                 ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                       <  1:3931434/71‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCATCCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                   <  1:2740702/58‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                   <  1:3753135/58‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                   <  1:3593097/58‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                   <  1:2102384/58‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                      <  1:1548144/71‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                      <  1:255345/71‑1 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                    >  1:271456/1‑70 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                    >  1:3899270/1‑70 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                    >  1:3777530/1‑70 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                    >  1:447005/1‑70 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                    >  1:344849/1‑70 (MQ=255)                       aTCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCAt                                                 <  1:3374726/71‑1 (MQ=255)                        tCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                >  1:3559816/1‑71 (MQ=255)                         cGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                               >  1:3071503/1‑71 (MQ=255)                          ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                                >  1:1746406/1‑53 (MQ=255)                          ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                              <  1:1785360/71‑1 (MQ=255)                           gTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAcgccgc                                                                    <  1:4080758/48‑1 (MQ=255)                            ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                               <  1:1368281/68‑1 (MQ=255)                             tCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGccgccg                                                                  >  1:3150346/1‑48 (MQ=255)                              cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                          >  1:962584/1‑71 (MQ=255)                              cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                          >  1:508069/1‑71 (MQ=255)                               aGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                         <  1:2526462/55‑1 (MQ=255)                                  cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                >  1:1254041/1‑61 (MQ=255)                                   ttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                     >  1:3006853/1‑71 (MQ=255)                                   ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                               <  1:2825671/45‑1 (MQ=255)                                    tCGGCCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                    >  1:1510627/1‑71 (MQ=255)                                    tCGGACCAACCTGCATTTTCAGTTTTTCGTCAACTCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                    >  1:2552676/1‑71 (MQ=255)                                        aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                <  1:2491041/71‑1 (MQ=255)                                        aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                <  1:1971130/71‑1 (MQ=255)                                        aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                <  1:1038274/71‑1 (MQ=255)                                        aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                <  1:2792016/71‑1 (MQ=255)                                         ccgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                    >  1:3826218/4‑66 (MQ=255)                                          cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtta                                 >  1:1234575/3‑67 (MQ=255)                                           aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                             <  1:3040962/71‑1 (MQ=255)                                           aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                             <  1:2910320/71‑1 (MQ=255)                                           aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                             <  1:2681248/71‑1 (MQ=255)                                           aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                             <  1:1016534/71‑1 (MQ=255)                                            aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                            >  1:1329070/1‑71 (MQ=255)                                             ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                           >  1:2043102/1‑71 (MQ=255)                                               tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                    <  1:1021816/44‑1 (MQ=255)                                               tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                    <  1:2522047/44‑1 (MQ=255)                                               tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                    <  1:1945260/44‑1 (MQ=255)                                               tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                          >  1:24167/1‑70 (MQ=255)                                                 cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACgcgc                         <  1:1895898/69‑1 (MQ=255)                                                     ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                     >  1:2381497/1‑69 (MQ=255)                                                     ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                     >  1:3654931/1‑69 (MQ=255)                                                     ttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                     >  1:1494808/1‑68 (MQ=255)                                                       cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                  <  1:102482/70‑1 (MQ=255)                                                       cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                  <  1:27440/70‑1 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                    <  1:1916967/68‑3 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                    <  1:3865116/68‑3 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                     <  1:3618961/67‑2 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                  <  1:1401060/70‑1 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                  <  1:1703505/70‑1 (MQ=255)                                                        aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                >  1:595649/1‑71 (MQ=255)                                                        aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                 <  1:2728383/70‑1 (MQ=255)                                                        aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                 <  1:3160660/70‑1 (MQ=255)                                                         gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc               <  1:2532648/71‑1 (MQ=255)                                                          tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGAACGCGCACCGAAGAACt              >  1:2008479/1‑71 (MQ=255)                                                           ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg             >  1:3852655/1‑71 (MQ=255)                                                           ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg             >  1:3491701/1‑71 (MQ=255)                                                            tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                            >  1:1668574/1‑55 (MQ=255)                                                            tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg             >  1:3339379/1‑70 (MQ=255)                                                             ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                  >  1:229808/1‑64 (MQ=255)                                                              tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                    >  1:1831769/1‑45 (MQ=16)                                                              tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                    >  1:3681981/1‑45 (MQ=16)                                                              tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa           <  1:2971132/70‑1 (MQ=255)                                                               cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGATCTGCAt          <  1:1270/70‑1 (MQ=255)                                                                gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt          <  1:52892/69‑1 (MQ=255)                                                                 tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg       <  1:2163221/71‑1 (MQ=255)                                                                   aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt     >  1:1856739/1‑71 (MQ=255)                                                                    aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt    <  1:3815239/71‑1 (MQ=255)                                                                      gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                              >  1:1648793/1‑43 (MQ=255)                                                                      gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                              >  1:2535825/1‑43 (MQ=255)                                                                      gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc  >  1:2708355/1‑71 (MQ=255)                                                                      |                                                                       ATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAC  >  minE/608835‑608974

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