Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,004	G→A	8.9%	D750D (GAC→GAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,004	0	G	A	8.9%	252.8 / 10.5	90	D750D (GAC→GAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (54/28);  new base A (3/5);  total (57/33)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCG  >  minE/607938‑608074                                                                   |                                                                       cGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGt                                                                       <  1:258612/68‑1 (MQ=255) cGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCct                                                                    >  1:1132147/1‑71 (MQ=255)   cTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgc                                                                  >  1:5226863/1‑71 (MQ=255)    tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGt                                                                       >  1:4540037/1‑65 (MQ=255)     ttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctg                                                                >  1:2236666/1‑71 (MQ=255)      tttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgc                                                               >  1:4449022/1‑71 (MQ=255)        tACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgctg                                                             >  1:3002028/1‑71 (MQ=255)           tGTACTATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctg                                                                <  1:970631/65‑1 (MQ=255)             tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGtt                                                           >  1:3642488/1‑68 (MQ=255)             tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCtt                                                        >  1:2897703/1‑71 (MQ=255)                gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAg                                                     <  1:4346096/71‑1 (MQ=255)                 aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGt                                                            >  1:1604128/1‑63 (MQ=255)                 aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGt                                                            >  1:1914988/1‑63 (MQ=255)                 aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGt                                                            >  1:1146529/1‑63 (MQ=255)                  tttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTc                                                   <  1:355191/71‑1 (MQ=255)                   ttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgct                                                              >  1:3999315/1‑59 (MQ=255)                   ttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgct                                                              >  1:2868653/1‑59 (MQ=255)                      aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGc                                                <  1:2413541/70‑1 (MQ=255)                      aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGc                                                <  1:2762980/70‑1 (MQ=255)                          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCct                                                                    >  1:1269137/1‑46 (MQ=255)                          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                           <  1:2151974/71‑1 (MQ=255)                          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                           >  1:229397/1‑71 (MQ=255)                          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                           <  1:645509/71‑1 (MQ=255)                          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                           >  1:4352840/1‑71 (MQ=255)                          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                           >  1:973397/1‑71 (MQ=255)                          aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                           >  1:2591745/1‑71 (MQ=255)                           aaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCtt                                                        <  1:894144/57‑1 (MQ=255)                            aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTc                                                                      >  1:919959/1‑42 (MQ=255)                            aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTc                                                                      >  1:380619/1‑42 (MQ=255)                            aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTc                                                                      >  1:5653643/1‑42 (MQ=255)                            aTCTAATTACATAGATTGACTGAAGGTACGAGTAATAACGt                                                                       >  1:1233973/1‑41 (MQ=255)                             tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGa                                              <  1:1309261/65‑1 (MQ=255)                              cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                       >  1:4477486/1‑71 (MQ=255)                              cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                       >  1:2170650/1‑71 (MQ=255)                              cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                       >  1:982748/1‑71 (MQ=255)                                 aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                       <  1:4419698/68‑1 (MQ=255)                                 aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                     >  1:2298041/1‑70 (MQ=255)                                 aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                     <  1:1962975/70‑1 (MQ=255)                                 aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                    <  1:3026259/71‑1 (MQ=255)                                  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                    >  1:4396549/1‑70 (MQ=255)                                  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                   >  1:2300988/1‑71 (MQ=255)                                  ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGATGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                    >  1:3494523/1‑70 (MQ=255)                                     cATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACt                                 <  1:4815956/70‑1 (MQ=255)                                      aTAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTa                                                      <  1:756690/48‑1 (MQ=255)                                      aTAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTa                                                      <  1:3725039/48‑1 (MQ=255)                                      aTAGATTGAGTGAAGGTACGAGTAATAAAGTCCTGCTGCTGTTCTTTa                                                      <  1:4201361/48‑1 (MQ=255)                                       tAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGGCAGCGAGTTGAAACGTACTGc                               <  1:3749529/70‑1 (MQ=255)                                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                            >  1:3197352/1‑56 (MQ=255)                                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                            >  1:71531/1‑56 (MQ=255)                                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                            >  1:2561674/1‑56 (MQ=255)                                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                            >  1:2338224/1‑56 (MQ=255)                                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                            >  1:2244981/1‑56 (MQ=255)                                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                            >  1:210154/1‑56 (MQ=255)                                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                            >  1:176854/1‑56 (MQ=255)                                        aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                              >  1:3536577/1‑70 (MQ=255)                                           ttGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                          <  1:2617571/71‑1 (MQ=255)                                            tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                     >  1:1953076/1‑59 (MQ=255)                                            tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                     >  1:1801633/1‑59 (MQ=255)                                             gAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGAGAGTTGAAACGTACTGCGTAGcc                         <  1:83408/70‑1 (MQ=255)                                               gtaaaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGa                      <  1:350587/68‑1 (MQ=255)                                                tGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                    <  1:1294776/56‑1 (MQ=255)                                                  aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                    <  1:2969174/70‑3 (MQ=255)                                                  aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga                    <  1:404675/70‑3 (MQ=255)                                                  aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                    >  1:227011/1‑70 (MQ=255)                                                  aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                    >  1:5315911/1‑70 (MQ=255)                                                   aGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                     >  1:805418/1‑52 (MQ=255)                                                   aGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                     >  1:5010636/1‑52 (MQ=255)                                                      tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                 >  1:3279435/1‑69 (MQ=255)                                                      tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                >  1:398507/1‑70 (MQ=255)                                                          aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg           <  1:87811/71‑1 (MQ=255)                                                          aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg           <  1:3814349/71‑1 (MQ=255)                                                            taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg               <  1:3326983/65‑1 (MQ=255)                                                            taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg               <  1:5158368/65‑1 (MQ=255)                                                            taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                <  1:366473/64‑1 (MQ=255)                                                             aataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGCATTGTCa        <  1:569736/71‑1 (MQ=255)                                                              ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt             >  1:2715948/1‑65 (MQ=255)                                                              ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt             >  1:531529/1‑65 (MQ=255)                                                              ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt             >  1:5538430/1‑65 (MQ=255)                                                              ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt             >  1:3593189/1‑65 (MQ=255)                                                              ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg       <  1:472346/71‑1 (MQ=255)                                                              ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg       >  1:1064665/1‑71 (MQ=255)                                                              ataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAg       >  1:4040379/1‑71 (MQ=255)                                                               taaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa        <  1:3450180/69‑1 (MQ=255)                                                                aaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                   >  1:3235411/1‑41 (MQ=255)                                                                aaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                   >  1:739180/1‑41 (MQ=255)                                                                aaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                   >  1:2849456/1‑41 (MQ=255)                                                                aaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                           <  1:1961296/49‑1 (MQ=255)                                                                 aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                             >  1:1115134/1‑46 (MQ=255)                                                                 aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg               >  1:2975706/1‑60 (MQ=255)                                                                 aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa        >  1:2340301/1‑67 (MQ=255)                                                                   gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg  >  1:5527442/1‑71 (MQ=255)                                                                   gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg  >  1:3563115/1‑71 (MQ=255)                                                                   |                                                                       CGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCG  >  minE/607938‑608074

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