Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,901	A→C	13.4%	V451V (GTT→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,901	0	A	C	13.4%	188.4 / 30.0	89	V451V (GTT→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (34/43);  new base C (6/6);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTT  >  minE/608835‑608971                                                                   |                                                                       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:1245158/71‑1 (MQ=255)        aGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGtt                                                             <  1:368580/71‑1 (MQ=255)         ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                            <  1:3280851/71‑1 (MQ=255)         ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                            <  1:355522/71‑1 (MQ=255)         ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                            <  1:1808459/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)                gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAc                                                                      >  1:2949624/1‑54 (MQ=255)                gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                      <  1:2592926/70‑1 (MQ=255)                gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                      >  1:3643924/1‑70 (MQ=255)                 ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                    <  1:3931434/71‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCATCCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                <  1:2740702/58‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                <  1:3593097/58‑1 (MQ=255)                  cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACgccgcc                                                                <  1:3753135/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)                   tttGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCaa                                                                       <  1:1231372/50‑1 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 >  1:271456/1‑70 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 >  1:447005/1‑70 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 >  1:344849/1‑70 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 >  1:3777530/1‑70 (MQ=255)                     tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 >  1:3899270/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)                                   ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAAc                                                                      <  1:2367990/35‑1 (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:1038274/71‑1 (MQ=255)                                        aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                             <  1:2491041/71‑1 (MQ=255)                                        aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                             <  1:1971130/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:1016534/71‑1 (MQ=255)                                           aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                          <  1:3040962/71‑1 (MQ=255)                                           aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                          <  1:2910320/71‑1 (MQ=255)                                           aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                          <  1:2681248/71‑1 (MQ=255)                                            aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                         >  1:1329070/1‑71 (MQ=255)                                             ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                        >  1:2043102/1‑71 (MQ=255)                                               tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 <  1:2522047/44‑1 (MQ=255)                                               tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                 <  1:1021816/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:3654931/1‑69 (MQ=255)                                                     ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                  >  1:2381497/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:3865116/68‑3 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                 <  1:1916967/68‑3 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                  <  1:3618961/67‑2 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag               <  1:1703505/70‑1 (MQ=255)                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag               <  1:1401060/70‑1 (MQ=255)                                                        aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa             >  1:595649/1‑71 (MQ=255)                                                        aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga              <  1:3160660/70‑1 (MQ=255)                                                        aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga              <  1:2728383/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:3681981/1‑45 (MQ=16)                                                              tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                 >  1:1831769/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)                                                                   |                                                                       ATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTT  >  minE/608835‑608971

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