Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,801	G→A	10.3%	D151D (GAC→GAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,801	0	G	A	10.3%	245.1 / 15.1	87	D151D (GAC→GAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (36/42);  new base A (6/3);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGC  >  minE/609731‑609871                                                                       |                                                                       tgccaGATGCATATGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                                        >  1:2894211/2‑71 (MQ=255) ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                                        >  1:96175/1‑71 (MQ=255) ggccaTATGCATCTGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                                        <  1:1154638/71‑1 (MQ=255)        tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                 <  1:1557305/71‑1 (MQ=255)        tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                 <  1:3863764/71‑1 (MQ=255)        tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                 <  1:3713503/71‑1 (MQ=255)          cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACg                                                               <  1:3294142/71‑1 (MQ=255)           aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATcaa                                                                    <  1:2573065/65‑3 (MQ=255)           aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGc                                                              >  1:934271/1‑71 (MQ=255)           aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                                    >  1:622401/1‑65 (MQ=255)           aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                                    >  1:3472631/1‑65 (MQ=255)             cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                            >  1:1200076/1‑71 (MQ=255)             cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                            >  1:683712/1‑71 (MQ=255)             cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                            >  1:3202735/1‑71 (MQ=255)             cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                            >  1:3124160/1‑71 (MQ=255)               ggCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                                    <  1:3232434/61‑1 (MQ=255)                gCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                                    <  1:2694772/60‑1 (MQ=255)                   aaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGtta                                                      >  1:3413941/2‑70 (MQ=255)                    aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                   <  1:540759/57‑1 (MQ=255)                    aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                     <  1:794064/71‑1 (MQ=255)                    aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                     >  1:1962760/1‑71 (MQ=255)                     ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg                                                    <  1:1080095/71‑1 (MQ=255)                      cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                                   <  1:472274/71‑1 (MQ=255)                       ggTCAGAACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCACTGGTtatggg                                                  >  1:2884392/1‑69 (MQ=255)                        gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                                                   >  1:2894459/1‑69 (MQ=255)                         tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGGCGAACACGCCCTGGTTGTGAGtt                                                <  1:1754112/71‑1 (MQ=255)                          cAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                               >  1:10591/1‑71 (MQ=255)                           aGAACTCCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                                              <  1:2585645/71‑1 (MQ=255)                           aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTTTAAACGTCGAACACGCCCTGGGTGTGAGtttt                                              <  1:3339867/71‑1 (MQ=255)                           aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                                     <  1:3287562/48‑1 (MQ=255)                           aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                               <  1:186373/70‑1 (MQ=255)                           aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                               <  1:2490538/70‑1 (MQ=255)                           aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                                               <  1:858471/70‑1 (MQ=255)                             aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                                             <  1:816445/70‑1 (MQ=255)                              acgccAGATTTACGGCATCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                           >  1:1815777/4‑71 (MQ=255)                               cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                                          >  1:1129451/1‑71 (MQ=255)                               cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                                          >  1:437384/1‑71 (MQ=255)                                gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGttt                                               <  1:1631026/64‑1 (MQ=255)                                gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                           >  1:2061190/2‑69 (MQ=255)                                  cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                                       >  1:1531336/1‑71 (MQ=255)                                   aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                                       >  1:1452631/1‑54 (MQ=255)                                    gATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                           >  1:1326924/1‑65 (MQ=255)                                    gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACg                                                               >  1:1028313/1‑45 (MQ=255)                                    gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGt                                         >  1:1099139/1‑67 (MQ=255)                                    gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                                      >  1:593872/1‑70 (MQ=255)                                    gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc                                     <  1:1668607/71‑1 (MQ=255)                                       ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                                          >  1:2056657/1‑47 (MQ=255)                                       ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                  >  1:1708068/1‑71 (MQ=255)                                       ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                  >  1:3031529/1‑71 (MQ=255)                                        tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                                                     <  1:627472/35‑1 (MQ=255)                                        tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                  <  1:2080340/70‑1 (MQ=255)                                         aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                                 <  1:2314118/70‑1 (MQ=255)                                         aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg                                 <  1:3194530/70‑1 (MQ=255)                                          cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGa                                <  1:1607697/70‑1 (MQ=255)                                              aaCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtt                           >  1:1050678/1‑71 (MQ=255)                                               aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                            >  1:1502410/1‑37 (MQ=255)                                               aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                            >  1:649819/1‑37 (MQ=255)                                               aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGa                             >  1:3117557/1‑68 (MQ=255)                                                 gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAg                                                  <  1:82200/45‑1 (MQ=255)                                                 gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                        <  1:1930890/71‑1 (MQ=255)                                                  cAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                       >  1:1831887/1‑71 (MQ=255)                                                   aGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttttt                      >  1:1797344/1‑71 (MQ=255)                                                    ggATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                     <  1:1372368/71‑1 (MQ=255)                                                     gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                    <  1:2008984/71‑1 (MQ=255)                                                     gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                    <  1:2074140/71‑1 (MQ=255)                                                     gATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                    <  1:2937464/71‑1 (MQ=255)                                                      aTGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat                   >  1:3690858/1‑71 (MQ=255)                                                       tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatc                  <  1:1569361/71‑1 (MQ=255)                                                       tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat                   >  1:2696638/1‑70 (MQ=255)                                                          ccGGAGTGTAAACGTCTAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                  <  1:1123984/52‑1 (MQ=255)                                                           cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcatc               <  1:3007321/70‑1 (MQ=255)                                                           cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcatc               <  1:392250/70‑1 (MQ=255)                                                           cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg              <  1:418031/71‑1 (MQ=255)                                                           cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg              >  1:622611/1‑71 (MQ=255)                                                            ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg             <  1:351944/71‑1 (MQ=255)                                                            ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg             <  1:2932642/71‑1 (MQ=255)                                                            ggAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCgg             <  1:3602155/71‑1 (MQ=255)                                                             gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg                                   <  1:1069834/48‑1 (MQ=255)                                                             gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGa            >  1:985011/1‑71 (MQ=255)                                                                tgtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATcc         <  1:2847058/71‑1 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttt                         <  1:3734233/53‑1 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAg       >  1:3697829/1‑71 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAg       >  1:2603706/1‑71 (MQ=255)                                                                    aaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGtt     >  1:3641481/1‑71 (MQ=255)                                                                     aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGa            >  1:3677054/1‑63 (MQ=255)                                                                     aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc    >  1:1798891/1‑71 (MQ=255)                                                                     aCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTc    >  1:3148315/1‑71 (MQ=255)                                                                      cGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                        <  1:1078831/50‑1 (MQ=255)                                                                      cGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcg   <  1:3790389/71‑1 (MQ=255)                                                                      cGTCGAACACGCCCTGGTTGTGAGATTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcg   >  1:1820147/1‑71 (MQ=255)                                                                       gTCGAACACGCCCTGGTTGTGAGTTTTAGGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcgc  >  1:1263605/1‑71 (MQ=255)                                                                       gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcgc  >  1:253230/1‑71 (MQ=255)                                                                       gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcg   >  1:3811134/1‑70 (MQ=255)                                                                       |                                                                       GGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGC  >  minE/609731‑609871

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