Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,087	A→C	5.1%	A389A (GCT→GCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,087	0	A	C	5.1%	215.8 / 4.0	79	A389A (GCT→GCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base C (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGAC  >  minE/609018‑609156                                                                      |                                                                      gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                                                       <  1:845802/70‑1 (MQ=255) gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg                                                                      >  1:5548385/1‑71 (MQ=255)   tgttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                                                    >  1:456444/1‑71 (MQ=255)            ccGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAt                                                              <  1:1150458/68‑1 (MQ=255)             cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                                           <  1:5337889/70‑1 (MQ=255)             cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                                           <  1:2198931/70‑1 (MQ=255)               gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                        <  1:5231204/71‑1 (MQ=255)               gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                        <  1:3077275/71‑1 (MQ=255)               gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGAAACTTTAGCGGCGAATTTCtt                                                         >  1:249047/1‑70 (MQ=255)                aCGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                                                       <  1:3909367/71‑1 (MQ=255)                 cGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                        >  1:1036772/1‑69 (MQ=255)                 cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTg                                                        >  1:5175063/1‑69 (MQ=255)                 cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGa                                                       <  1:1731496/70‑1 (MQ=255)                   cATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                                                    >  1:1500114/1‑55 (MQ=255)                   cATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGAAgg                                                    <  1:1224327/71‑2 (MQ=255)                    aTCAGATCGTCGTTTTCATACTGCAAAGAGGATGTGACGATAGACACCTGCGCGGCATATTTTTTGaa                                                      >  1:3717520/1‑68 (MQ=255)                     tCAGTTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                                  <  1:3847743/71‑1 (MQ=255)                     tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                        <  1:2152486/65‑1 (MQ=255)                     tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                                  <  1:2073074/71‑1 (MQ=255)                     tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                                  <  1:4737692/71‑1 (MQ=255)                     tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                                  <  1:4503201/71‑1 (MQ=255)                      cAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGtt                                                   >  1:4818782/1‑69 (MQ=255)                       aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                                      >  1:1385118/1‑65 (MQ=255)                       aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                                      >  1:2111126/1‑65 (MQ=255)                       aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                                                       >  1:2873580/1‑64 (MQ=255)                         gTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                                                       <  1:5636241/46‑1 (MQ=255)                          tCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                             <  1:4179864/71‑1 (MQ=255)                           cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                             >  1:5694392/1‑70 (MQ=255)                            aTCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGc                                            <  1:5530304/70‑1 (MQ=255)                             tCGTTCTCATTCTGCAGAGAAGAGGTGTCGATGGACACTTTAGGGGCGAATttgttg                                                        <  1:3695301/57‑1 (MQ=255)                             tCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAg                                          <  1:1416829/71‑1 (MQ=255)                               gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGGGGCAGtt                                        >  1:145938/1‑71 (MQ=255)                               gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtt                                        >  1:3218500/1‑71 (MQ=255)                                ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                                       <  1:701365/71‑1 (MQ=255)                                ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                                       <  1:5655573/71‑1 (MQ=255)                                  ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                               >  1:3134255/1‑45 (MQ=255)                                  ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                               >  1:2721395/1‑45 (MQ=255)                                  ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                                      <  1:5499503/70‑1 (MQ=255)                                  ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                                      >  1:1640732/1‑70 (MQ=255)                                  ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                                      <  1:986105/70‑1 (MQ=255)                                   tcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                        <  1:2005458/51‑1 (MQ=255)                                   tcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                                  <  1:2449680/57‑1 (MQ=255)                                    cATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                                     >  1:3973005/1‑69 (MQ=255)                                    cATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAATTCTTGAAGTTCAGCGGCAGttttt                                     >  1:3973945/1‑69 (MQ=255)                                     aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                                     <  1:4208587/68‑1 (MQ=255)                                     aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                                     <  1:2829398/68‑1 (MQ=255)                                           cAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGt                                                    >  1:3659735/1‑47 (MQ=255)                                           cAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACca                              >  1:4997603/1‑69 (MQ=255)                                            agagAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAg                           <  1:2205120/71‑1 (MQ=255)                                             gagaAGAGGTGTCGATGGGCACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGa                          <  1:2248125/71‑1 (MQ=255)                                             gagaAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCa                                           >  1:3008616/1‑54 (MQ=255)                                             gagaAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCa                                           >  1:439162/1‑54 (MQ=255)                                              agaagaGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAg                           <  1:1388042/69‑1 (MQ=255)                                                 agagGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATgg                      <  1:143769/71‑1 (MQ=255)                                                     gtgtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGt                     <  1:944340/68‑1 (MQ=255)                                                      tgtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                             <  1:5172726/43‑1 (MQ=255)                                                       gtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGt                >  1:5575644/1‑71 (MQ=255)                                                        tCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGtt               >  1:3521916/1‑71 (MQ=255)                                                        tCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGt                <  1:2777022/70‑1 (MQ=255)                                                         cGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGaa                         >  1:5057120/1‑60 (MQ=255)                                                          gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAg                                                >  1:4742184/1‑36 (MQ=255)                                                          gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATg                 <  1:4558037/67‑1 (MQ=255)                                                          gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATg                 <  1:4903770/67‑1 (MQ=255)                                                          gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATg                 <  1:3906427/67‑1 (MQ=255)                                                          gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCg             >  1:1249222/1‑71 (MQ=255)                                                            tGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACcac                             <  1:4296934/53‑1 (MQ=255)                                                             ggACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGtt          >  1:908052/1‑71 (MQ=255)                                                              gACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCg             >  1:993016/1‑67 (MQ=255)                                                              gACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTc         >  1:1119123/1‑71 (MQ=255)                                                               acacTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                                       <  1:842993/40‑1 (MQ=255)                                                                cacTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCg       <  1:3235862/71‑1 (MQ=255)                                                                 acTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCgg      <  1:5056921/71‑1 (MQ=255)                                                                 acTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCgg      <  1:4839738/71‑1 (MQ=255)                                                                  ctatAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCgg            <  1:1968538/61‑1 (MQ=255)                                                                  cTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGtt          >  1:196748/1‑66 (MQ=255)                                                                   tttAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                                      <  1:65577/37‑1 (MQ=255)                                                                    ttAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGAc                                >  1:4845067/1‑42 (MQ=255)                                                                    ttAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGAc                                >  1:687612/1‑42 (MQ=255)                                                                    ttAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTc              >  1:3949843/1‑60 (MQ=255)                                                                    ttAGCGGCGAATTTCTTGAAGTACAGCGGCAGTTTTTCAGAc                                >  1:3069775/1‑42 (MQ=255)                                                                     tAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGAc  >  1:2382759/1‑71 (MQ=255)                                                                      |                                                                      GTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGAC  >  minE/609018‑609156

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