Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,783	G→C	16.6%	I157M (ATC→ATG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,783	0	G	C	16.6%	199.0 / 29.7	84	I157M (ATC→ATG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (44/26);  new base C (9/5);  total (53/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTG  >  minE/609714‑609852                                                                      |                                                                      aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:2124549/1‑71 (MQ=255) aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:1913836/1‑71 (MQ=255)   cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                      >  1:589932/1‑69 (MQ=255)   cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                                                                    <  1:906833/71‑1 (MQ=255)    gatgatACGGCCACGGCCATATGCATCTGGCAGCCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                   >  1:1728849/1‑71 (MQ=255)    gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                                                     >  1:1716252/1‑69 (MQ=255)    gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                   >  1:426803/1‑71 (MQ=255)    gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                   >  1:170203/1‑71 (MQ=255)         tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                                                       >  1:312090/1‑62 (MQ=255)             gccacggccaTATGCATCTGGCAGACCGGTCAGAAGACCAGATTTACGGCAACGCAGGATGTCCGGAgtg                                                           >  1:530807/1‑70 (MQ=255)              ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                         <  1:3791377/71‑1 (MQ=255)               cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCggtga                                                            >  1:1923519/1‑64 (MQ=255)               cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                         <  1:293023/70‑1 (MQ=255)                  ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                     <  1:2958091/71‑1 (MQ=255)                   gccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                       >  1:2523725/1‑68 (MQ=255)                      ttatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                         <  1:4005554/62‑1 (MQ=255)                      atatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                 <  1:1666750/71‑1 (MQ=255)                       tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                <  1:1281979/71‑1 (MQ=255)                        atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg                                                               <  1:3523259/55‑1 (MQ=255)                         tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                                                             >  1:2575142/1‑56 (MQ=255)                         tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:2104948/1‑71 (MQ=255)                         tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:2561492/1‑71 (MQ=255)                         tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:2272082/1‑71 (MQ=255)                         tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:1967169/1‑71 (MQ=255)                         tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:1758809/1‑71 (MQ=255)                         gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                >  1:721075/3‑69 (MQ=255)                         gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                >  1:1757432/3‑69 (MQ=255)                         gccATCTGGTAAACAGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                >  1:1214781/3‑69 (MQ=255)                          gCATCTGGCAGCCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                                                  >  1:958567/1‑50 (MQ=255)                           cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATca                                                  >  1:1686080/1‑64 (MQ=255)                           cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                                                       >  1:541739/1‑44 (MQ=255)                             tCTGGCAGAGCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg                                                               <  1:2699775/50‑1 (MQ=255)                             tCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGcc                                          >  1:2871227/1‑71 (MQ=255)                             tCTGGCAGACC‑GTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                         >  1:3835619/1‑71 (MQ=255)                              cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                       >  1:3630091/1‑57 (MQ=255)                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                                                     >  1:674260/1‑42 (MQ=255)                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                                                     >  1:679048/1‑42 (MQ=255)                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                       >  1:4131320/1‑56 (MQ=255)                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                       >  1:4108083/1‑56 (MQ=255)                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                        >  1:1463254/1‑71 (MQ=255)                               tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                        >  1:2667979/1‑71 (MQ=255)                                  cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGt                                     >  1:3626158/1‑71 (MQ=255)                                   aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt                                     >  1:2203540/3‑70 (MQ=255)                                   aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt                                     >  1:588441/3‑70 (MQ=255)                                   aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                    >  1:2484232/1‑71 (MQ=255)                                   aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                    >  1:583371/1‑71 (MQ=255)                                    gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                <  1:482496/58‑1 (MQ=255)                                    gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                <  1:3222016/58‑1 (MQ=255)                                    gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtg                                   >  1:2602522/1‑71 (MQ=255)                                      ccGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTtatg                                 >  1:916180/1‑71 (MQ=255)                                      ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                  >  1:635880/1‑70 (MQ=255)                                       cGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                   <  1:1073497/36‑1 (MQ=255)                                       cGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                   <  1:3565293/36‑1 (MQ=255)                                       cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                  >  1:1938058/1‑69 (MQ=255)                                       cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                  >  1:1610704/1‑69 (MQ=255)                                          tCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                                        <  1:1840715/44‑1 (MQ=12)                                             gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                              >  1:2404495/1‑51 (MQ=255)                                             gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                           <  1:32959/70‑1 (MQ=255)                                              aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                       >  1:3022398/1‑41 (MQ=255)                                              aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAc                         >  1:677873/1‑71 (MQ=255)                                               acgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                        <  1:241043/68‑1 (MQ=255)                                               acacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                        <  1:1709024/71‑1 (MQ=255)                                               acacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                        <  1:3980117/71‑1 (MQ=255)                                                cgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                        <  1:1354525/68‑1 (MQ=255)                                                cgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                        >  1:952727/3‑70 (MQ=255)                                                  ccAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                <  1:2034711/44‑1 (MQ=255)                                                  ccAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTa                     >  1:2956282/1‑71 (MQ=255)                                                   cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAGACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                    >  1:3123276/1‑71 (MQ=255)                                                   cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                    >  1:1292000/1‑71 (MQ=255)                                                   cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                    >  1:3127472/1‑71 (MQ=255)                                                    aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACg                                            >  1:3724182/1‑46 (MQ=255)                                                    aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                   >  1:2109699/1‑71 (MQ=255)                                                     gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTgg                                      <  1:746411/51‑1 (MQ=255)                                                      aTTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGttt                            <  1:3061216/60‑1 (MQ=12)                                                      aTTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCa                 <  1:2775210/71‑1 (MQ=255)                                                        ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt               >  1:3992948/1‑71 (MQ=255)                                                        ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt               >  1:838025/1‑71 (MQ=255)                                                           cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg                <  1:2792215/67‑1 (MQ=255)                                                           cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg                <  1:1061055/67‑1 (MQ=255)                                                             gCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                <  1:577511/49‑1 (MQ=255)                                                              cAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAt         <  1:787331/71‑1 (MQ=255)                                                                 cGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt      <  1:180482/71‑1 (MQ=255)                                                                  gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt     <  1:1141544/71‑1 (MQ=255)                                                                  gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt     <  1:2791061/71‑1 (MQ=255)                                                                  gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt      >  1:3853816/1‑70 (MQ=255)                                                                    aGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg              <  1:3829167/60‑1 (MQ=255)                                                                     ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg  <  1:2070941/71‑1 (MQ=255)                                                                      gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg              <  1:1436360/58‑1 (MQ=255)                                                                      |                                                                      ACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTG  >  minE/609714‑609852

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