Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,280	A→G	9.1%	A658A (GCT→GCC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,280	0	A	G	9.1%	178.8 / 7.5	66	A658A (GCT→GCC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base G (2/4);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGT  >  minE/608212‑608350                                                                     |                                                                       aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg                                                                       >  1:3709296/1‑70 (MQ=255)  cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:1712302/71‑1 (MQ=255)  cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:2255876/71‑1 (MQ=255)     cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                                  >  1:1025265/1‑71 (MQ=255)     cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                                  >  1:3407785/1‑71 (MQ=255)      gtcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                                  <  1:3206845/70‑1 (MQ=255)         gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                                                >  1:2550577/1‑69 (MQ=255)          tcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACAATCTTTAGCGTAAGCAAACGGCAg                                                             <  1:1804639/71‑1 (MQ=255)           cTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGt                                                            >  1:1149312/1‑71 (MQ=255)              tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGttt                                                          <  1:3813684/70‑1 (MQ=255)              tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGttt                                                          <  1:821148/70‑1 (MQ=255)                cccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                         <  1:1176120/69‑1 (MQ=255)                  cAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                         >  1:1048478/1‑67 (MQ=255)                   aGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                    >  1:7101/1‑71 (MQ=255)                    gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                   <  1:3339879/71‑1 (MQ=255)                     tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  <  1:269442/71‑1 (MQ=255)                     tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  <  1:1722830/71‑1 (MQ=255)                     tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTATAGc                                                       >  1:2105193/1‑66 (MQ=255)                       cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  >  1:3867587/1‑69 (MQ=255)                       cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  >  1:3834838/1‑69 (MQ=255)                        gTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAg                                               >  1:369184/1‑71 (MQ=255)                          tCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                              >  1:361587/1‑70 (MQ=255)                           cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACg                                                                 >  1:101027/1‑50 (MQ=255)                           cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                              >  1:3140270/1‑69 (MQ=255)                            ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                     >  1:1670850/1‑61 (MQ=255)                            ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           <  1:2044591/71‑1 (MQ=255)                            ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           <  1:694873/71‑1 (MQ=255)                            ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           <  1:3333443/71‑1 (MQ=255)                             gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           <  1:3827985/70‑1 (MQ=255)                              aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                            >  1:3352267/1‑68 (MQ=255)                              aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                            >  1:1412528/1‑68 (MQ=255)                              aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                            >  1:1982606/1‑68 (MQ=255)                               aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                          >  1:2992609/1‑67 (MQ=255)                               aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                        <  1:324297/71‑1 (MQ=255)                               aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          >  1:3302705/1‑69 (MQ=255)                               aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          >  1:3283272/1‑69 (MQ=255)                                  attgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                     <  1:267180/68‑1 (MQ=255)                                  attgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                     <  1:2742794/68‑1 (MQ=255)                                  aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     <  1:2445083/71‑1 (MQ=255)                                   tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  >  1:1607299/1‑57 (MQ=255)                                   tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     >  1:428323/1‑70 (MQ=255)                                    agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagag                                         >  1:2221826/1‑65 (MQ=255)                                      agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     <  1:2212037/67‑1 (MQ=255)                                      agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACt                                 >  1:2131970/1‑71 (MQ=255)                                        aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                >  1:523529/1‑54 (MQ=255)                                        aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                >  1:1123783/1‑70 (MQ=255)                                            tgtACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGac                                                <  1:2837786/50‑2 (MQ=255)                                             gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                          >  1:953217/1‑71 (MQ=255)                                              tacgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCttt                         <  1:1753908/68‑1 (MQ=255)                                                ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                        <  1:2401660/70‑1 (MQ=255)                                                    aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                   <  1:1654394/71‑1 (MQ=255)                                                      aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                              <  1:2029132/42‑1 (MQ=255)                                                      aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                   >  1:1331851/1‑69 (MQ=255)                                                        cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGa               >  1:720036/1‑71 (MQ=255)                                                         aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     <  1:3942664/48‑1 (MQ=255)                                                         aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACATTTCTGGTCACGAc              >  1:3260619/1‑71 (MQ=255)                                                          tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAc              >  1:1051499/1‑70 (MQ=255)                                                            tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt           >  1:1980820/1‑71 (MQ=255)                                                              tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGTACCTTTCTGGTCACGACCGTg          >  1:87529/1‑70 (MQ=255)                                                              tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc             <  1:2037412/67‑1 (MQ=255)                                                                 cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc      <  1:3250035/71‑1 (MQ=255)                                                                 cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTAACGACCGTGc         <  1:913851/68‑1 (MQ=255)                                                                 cGTAAGAAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc      <  1:3164314/71‑1 (MQ=255)                                                                   tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                        <  1:3559758/51‑1 (MQ=255)                                                                   tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg    >  1:1019022/1‑71 (MQ=255)                                                                   tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTAATg                                <  1:355408/43‑1 (MQ=255)                                                                     aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                               >  1:2251034/1‑42 (MQ=255)                                                                     aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCATCTTTCTGGTCACGACCGTGCATCGGGt  <  1:3696793/71‑1 (MQ=255)                                                                     aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt  >  1:850934/1‑71 (MQ=255)                                                                     aGCAAACGGCAGTTTAGCAACGGAAGGCAGCGCGGCTACTGCACCTTTCTGGTCACGACCGTGCATCggg   <  1:3344432/70‑1 (MQ=255)                                                                     |                                                                       ACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGT  >  minE/608212‑608350

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