Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,786	G→A	17.3%	D156D (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,786	0	G	A	17.3%	200.7 / 31.2	81	D156D (GAC→GAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (39/28);  new base A (9/5);  total (48/33)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCA  >  minE/609716‑609856                                                                       |                                                                       cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                                                                        <  1:906833/71‑1 (MQ=255)  gatgatACGGCCACGGCCATATGCATCTGGCAGCCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                       >  1:1728849/1‑71 (MQ=255)  gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                       >  1:170203/1‑71 (MQ=255)  gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                                       >  1:426803/1‑71 (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:1758809/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:2104948/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)                        gCATCTGGCAGCCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                                                      >  1:958567/1‑50 (MQ=255)                         cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATca                                                      >  1:1686080/1‑64 (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)                             tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                           >  1:4108083/1‑56 (MQ=255)                             tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                           >  1:4131320/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:583371/1‑71 (MQ=255)                                 aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                                        >  1:2484232/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: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: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: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: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)                                                                       |                                                                       CGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATTTTTTTGATCA  >  minE/609716‑609856

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