Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,084	T→C	5.4%	K390K (AAA→AAG) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,084	0	T	C	5.4%	212.6 / 2.0	74	K390K (AAA→AAG) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (30/40);  new base C (2/2);  total (32/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 = 8.57e-01

ATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGG  >  minE/609015‑609152                                                                      |                                                                     aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                     >  1:705335/1‑71 (MQ=255) aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                     >  1:3379225/1‑71 (MQ=255) aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                     >  1:445539/1‑71 (MQ=255) aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                     >  1:3179307/1‑71 (MQ=255)  tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                     >  1:3116851/1‑70 (MQ=255)  tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                     >  1:3110437/1‑70 (MQ=255)   cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                     >  1:2781174/1‑69 (MQ=255)    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:2198931/70‑1 (MQ=255)                cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                                       <  1:5337889/70‑1 (MQ=255)                  gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                    <  1:3077275/71‑1 (MQ=255)                  gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                                    <  1:5231204/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:2111126/1‑65 (MQ=255)                          aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                                  >  1:1385118/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)                              cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                                     >  1:3452629/1‑42 (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:5655573/71‑1 (MQ=255)                                   ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                                   <  1:701365/71‑1 (MQ=255)                                     ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                           >  1:2721395/1‑45 (MQ=255)                                     ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                           >  1:3134255/1‑45 (MQ=255)                                     ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                                  <  1:986105/70‑1 (MQ=255)                                     ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                                  <  1:5499503/70‑1 (MQ=255)                                     ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                                  >  1:1640732/1‑70 (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:4903770/67‑1 (MQ=255)                                                             gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATg             <  1:3906427/67‑1 (MQ=255)                                                             gATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATg             <  1:4558037/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)                                                                     cTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGtt      >  1:196748/1‑66 (MQ=255)                                                                      tttAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                                  <  1:65577/37‑1 (MQ=255)                                                                      |                                                                     ATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGG  >  minE/609015‑609152

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑