Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,928	G→C	19.7%	T442T (ACC→ACG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,928	0	G	C	19.7%	174.7 / 37.1	81	T442T (ACC→ACG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base C (8/8);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACG  >  minE/608858‑608994                                                                       |                                                                   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)           cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                                    >  1:1254041/1‑61 (MQ=255)            ttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                                         >  1:3006853/1‑71 (MQ=255)             tCGGCCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                                        >  1:1510627/1‑71 (MQ=255)             tCGGACCAACCTGCATTTTCAGTTTTTCGTCAACTCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                                        >  1:2552676/1‑71 (MQ=255)                 aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                    <  1:1971130/71‑1 (MQ=255)                 aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                    <  1:2491041/71‑1 (MQ=255)                 aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                    <  1:2792016/71‑1 (MQ=255)                 aCCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                    <  1:1038274/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)                        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:27440/70‑1 (MQ=255)                                cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                      <  1:102482/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)                                              ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                       <  1:2023686/69‑1 (MQ=255)                                              ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                       <  1:725150/69‑1 (MQ=255)                                               gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                  >  1:1648793/1‑43 (MQ=255)                                               gccgccGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                  >  1:2535825/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:953683/1‑70 (MQ=255)                                                    cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                  >  1:3075564/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:1852293/1‑71 (MQ=255)                                                        gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca             >  1:587614/1‑71 (MQ=255)                                                        gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc              <  1:4005748/70‑1 (MQ=255)                                                        gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTCGTTTACCAACGatca             >  1:685199/1‑71 (MQ=255)                                                          ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc           <  1:910646/71‑1 (MQ=255)                                                          ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc           <  1:1823578/71‑1 (MQ=255)                                                            gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaagaa                                    >  1:3238564/1‑44 (MQ=12)                                                            gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                            >  1:2117140/1‑36 (MQ=255)                                                                aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat            >  1:111947/1‑64 (MQ=255)                                                                aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAAGatcat            >  1:1450022/1‑64 (MQ=255)                                                                   gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt      >  1:1288974/1‑67 (MQ=255)                                                                   gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg  <  1:43335/71‑1 (MQ=255)                                                                   gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg  <  1:2648335/71‑1 (MQ=255)                                                                   gCATGGTTTTCGCCAGGTTTGCACCCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg  >  1:3004311/1‑71 (MQ=255)                                                                   gCATGGGTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAttt                            >  1:538428/1‑45 (MQ=255)                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                    <  1:1360880/50‑1 (MQ=255)                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg         <  1:2636176/61‑1 (MQ=255)                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg         <  1:4123172/61‑1 (MQ=255)                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg         <  1:1040688/61‑1 (MQ=255)                                                                      tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCGTACg  >  1:3760010/1‑68 (MQ=255)                                                                       ggTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                               >  1:119737/1‑38 (MQ=255)                                                                       |                                                                   TCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACG  >  minE/608858‑608994

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