Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,807	C→T	9.2%	V149V (GTG→GTA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,807	0	C	T	9.2%	233.2 / 10.7	76	V149V (GTG→GTA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (39/30);  new base T (4/3);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

TGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTG  >  minE/609738‑609876                                                                      |                                                                      tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                      >  1:2104948/1‑71 (MQ=255) tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                      >  1:1967169/1‑71 (MQ=255) tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                      >  1:2272082/1‑71 (MQ=255) tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                      >  1:1758809/1‑71 (MQ=255) tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                                                      >  1:2561492/1‑71 (MQ=255)     tCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGcc                                                                  >  1:2871227/1‑71 (MQ=255)     tCTGGCAGACC‑GTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                                                 >  1:3835619/1‑71 (MQ=255)       tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                                                >  1:2667979/1‑71 (MQ=255)       tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                                                >  1:1463254/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:2484232/1‑71 (MQ=255)           aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                                            >  1:583371/1‑71 (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:1610704/1‑69 (MQ=255)               cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                                                          >  1:1938058/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)                          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:1061055/67‑1 (MQ=255)                                   cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg                                        <  1:2792215/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: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:3825851/53‑1 (MQ=255)                                                               cGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                          <  1:156395/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)                                                                     acacGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTg  <  1:1031598/71‑1 (MQ=255)                                                                     acacGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTg  <  1:2557017/71‑1 (MQ=255)                                                                     acacGCCCTAGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTg                          <  1:3930945/47‑1 (MQ=255)                                                                      cacGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttt                               <  1:988322/41‑1 (MQ=255)                                                                      cacGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAttt                               <  1:3071578/41‑1 (MQ=255)                                                                      |                                                                      TGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCATCGGATCCAGTTCGCGGTTG  >  minE/609738‑609876

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