Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,913	G→A	18.5%	I447I (ATC→ATT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,913	0	G	A	18.5%	194.4 / 33.9	81	I447I (ATC→ATT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base A (7/8);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCA  >  minE/608843‑608983                                                                       |                                                                       ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                                        <  1:1808459/71‑1 (MQ=255) ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                                        <  1:355522/71‑1 (MQ=255) ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                                        <  1:3280851/71‑1 (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:2592926/70‑1 (MQ=255)        gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                                  >  1:3643924/1‑70 (MQ=255)         ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                                <  1:3931434/71‑1 (MQ=255)          cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                               <  1:255345/71‑1 (MQ=255)          cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                               <  1:1548144/71‑1 (MQ=255)             tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                             >  1:344849/1‑70 (MQ=255)             tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                             >  1:271456/1‑70 (MQ=255)             tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                             >  1:3777530/1‑70 (MQ=255)             tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                             >  1:3899270/1‑70 (MQ=255)             tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                             >  1:447005/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)                  ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                                       <  1:1785360/71‑1 (MQ=255)                    ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                                        <  1:1368281/68‑1 (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:2792016/71‑1 (MQ=255)                                aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                         <  1:1971130/71‑1 (MQ=255)                                aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                         <  1:1038274/71‑1 (MQ=255)                                aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                         <  1:2491041/71‑1 (MQ=255)                                 ccgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                             >  1:3826218/4‑66 (MQ=255)                                  cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtta                                          >  1:1234575/3‑67 (MQ=255)                                   aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                      <  1:2681248/71‑1 (MQ=255)                                   aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                      <  1:2910320/71‑1 (MQ=255)                                   aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                      <  1:3040962/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:1945260/44‑1 (MQ=255)                                       tGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                             <  1:2522047/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:1916967/68‑3 (MQ=255)                                               cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAccaa                             <  1:3865116/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:3491701/1‑71 (MQ=255)                                                   ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                      >  1:3852655/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)                                                             ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa            <  1:725150/69‑1 (MQ=255)                                                             ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa            <  1:2023686/69‑1 (MQ=255)                                                              gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                       >  1:2535825/1‑43 (MQ=255)                                                              gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                       >  1:1648793/1‑43 (MQ=255)                                                              gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc           >  1:2708355/1‑71 (MQ=255)                                                               ccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaa                            >  1:2585727/1‑53 (MQ=255)                                                                 gccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCaa        >  1:128065/1‑71 (MQ=255)                                                                   cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc       >  1:3075564/1‑70 (MQ=255)                                                                   cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc       >  1:953683/1‑70 (MQ=255)                                                                   cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg      <  1:677044/71‑1 (MQ=255)                                                                      tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                       >  1:2935016/1‑35 (MQ=255)                                                                      tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat    >  1:657558/1‑70 (MQ=255)                                                                       gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                      >  1:2819634/1‑35 (MQ=255)                                                                       gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca  >  1:587614/1‑71 (MQ=255)                                                                       gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca  >  1:1852293/1‑71 (MQ=255)                                                                       gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc   <  1:4005748/70‑1 (MQ=255)                                                                       gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTCGTTTACCAACGatca  >  1:685199/1‑71 (MQ=255)                                                                       |                                                                       GGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCA  >  minE/608843‑608983

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