Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,804	G→A	8.3%	F150F (TTC→TTT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,804	0	G	A	8.3%	262.7 / 10.4	84	F150F (TTC→TTT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (43/34);  new base A (4/3);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

ATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGT  >  minE/609735‑609874                                                                      |                                                                       atatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                                                       <  1:1666750/71‑1 (MQ=255)  tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                                                      <  1:1281979/71‑1 (MQ=255)    tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                    >  1:1967169/1‑71 (MQ=255)    tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                    >  1:2104948/1‑71 (MQ=255)    tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                    >  1:1758809/1‑71 (MQ=255)    tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                    >  1:2561492/1‑71 (MQ=255)    tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                    >  1:2272082/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)        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:588441/3‑70 (MQ=255)              aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt                                                           >  1:2203540/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: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: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:113296/1‑65 (MQ=255)                                                    gTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtttttt                          >  1:638912/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:2828502/64‑1 (MQ=255)                                                              tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCg                 <  1:3410771/64‑1 (MQ=255)                                                              tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAg          >  1:3723050/1‑71 (MQ=255)                                                              tAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAg          >  1:3723209/1‑71 (MQ=255)                                                                aaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                        <  1:1112008/55‑1 (MQ=255)                                                                aaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGtt        >  1:3585666/1‑71 (MQ=255)                                                                aaCGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGtt        >  1:218940/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)                                                                    tCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTcgcg    <  1:3860315/71‑1 (MQ=255)                                                                      gAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGatcat                   >  1:3603849/1‑54 (MQ=255)                                                                      gAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGt  >  1:390382/1‑71 (MQ=255)                                                                      gAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGt  <  1:123865/71‑1 (MQ=255)                                                                      |                                                                       ATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGT  >  minE/609735‑609874

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