Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,801	G→A	12.8%	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	12.8%	239.3 / 21.7	86	D151D (GAC→GAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (42/33);  new base A (8/3);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCG  >  minE/609731‑609870                                                                       |                                                                      ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                                                       <  1:2958091/71‑1 (MQ=255)     atatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                                   <  1:1666750/71‑1 (MQ=255)      tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                  <  1:1281979/71‑1 (MQ=255)        tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                >  1:2272082/1‑71 (MQ=255)        tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                >  1:2104948/1‑71 (MQ=255)        tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                >  1:2561492/1‑71 (MQ=255)        tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                >  1:1967169/1‑71 (MQ=255)        tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                >  1:1758809/1‑71 (MQ=255)        gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                                  >  1:1757432/3‑69 (MQ=255)        gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                                  >  1:721075/3‑69 (MQ=255)        gccATCTGGTAAACAGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                                                  >  1:1214781/3‑69 (MQ=255)          cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATca                                                                    >  1:1686080/1‑64 (MQ=255)            tCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGcc                                                            >  1:2871227/1‑71 (MQ=255)            tCTGGCAGACC‑GTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                           >  1:3835619/1‑71 (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:583371/1‑71 (MQ=255)                  aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                                      >  1:2484232/1‑71 (MQ=255)                   gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                  <  1:3222016/58‑1 (MQ=255)                   gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                  <  1:482496/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)                      cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                    >  1:1938058/1‑69 (MQ=255)                      cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                    >  1:1610704/1‑69 (MQ=255)                            gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                >  1:2404495/1‑51 (MQ=255)                            gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                                             <  1:32959/70‑1 (MQ=255)                             aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAc                                           >  1:677873/1‑71 (MQ=255)                              acgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                          <  1:241043/68‑1 (MQ=255)                              acacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                          <  1:3980117/71‑1 (MQ=255)                              acacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg                                          <  1:1709024/71‑1 (MQ=255)                               cgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                          >  1:952727/3‑70 (MQ=255)                               cgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                                          <  1:1354525/68‑1 (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:3127472/1‑71 (MQ=255)                                  cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt                                      >  1:1292000/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:838025/1‑71 (MQ=255)                                       ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt                                 >  1:3992948/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:2791061/71‑1 (MQ=255)                                                 gCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttt                       <  1:1141544/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)                                                      aTGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                   <  1:2790006/70‑1 (MQ=255)                                                       tGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat                  >  1:3114934/1‑70 (MQ=255)                                                        gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                      >  1:638912/1‑65 (MQ=255)                                                        gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                      >  1:113296/1‑65 (MQ=255)                                                        gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtatttt                      >  1:903910/1‑65 (MQ=255)                                                        gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatca                >  1:477762/1‑71 (MQ=255)                                                           cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                      >  1:2178093/1‑62 (MQ=255)                                                           cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcatc              <  1:107495/70‑1 (MQ=255)                                                           cGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg             >  1:1080820/1‑71 (MQ=255)                                                             gAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGa                   >  1:1803670/1‑63 (MQ=255)                                                               gtgtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATc         <  1:280475/71‑1 (MQ=255)                                                                tgtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                                         <  1:2589183/38‑1 (MQ=255)                                                                 gtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGt                                        >  1:3401071/1‑38 (MQ=255)                                                                 gtAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTa                                       >  1:33971/1‑39 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt                                     <  1:2543922/40‑1 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttttttt                     <  1:1892635/56‑1 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                    >  1:2979279/1‑57 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat                  <  1:2900973/59‑1 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg             <  1:3410771/64‑1 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg             <  1:2828502/64‑1 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAg      >  1:3723209/1‑71 (MQ=255)                                                                  tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAg      >  1:3723050/1‑71 (MQ=255)                                                                    aaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                    <  1:1112008/55‑1 (MQ=255)                                                                    aaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGtt    >  1:218940/1‑71 (MQ=255)                                                                    aaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGtt    >  1:3585666/1‑71 (MQ=255)                                                                      cGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                    <  1:156395/53‑1 (MQ=255)                                                                      cGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                    <  1:3825851/53‑1 (MQ=255)                                                                      cGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcg  <  1:1436682/71‑1 (MQ=255)                                                                       gTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGat                  >  1:2427187/1‑54 (MQ=255)                                                                       |                                                                      GGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCG  >  minE/609731‑609870

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