Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,723	A→G	14.9%	R177R (CGT→CGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,723	0	A	G	14.9%	160.7 / 19.5	67	R177R (CGT→CGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (37/20);  new base G (3/7);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 7.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGG  >  minE/609657‑609791                                                                   |                                                                     tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCTCAACGCGACGGTAGTCACCGATGATacgg                                                                   <  1:3307852/70‑1 (MQ=255) tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacgg                                                                   <  1:2629354/70‑1 (MQ=255) tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                                    <  1:404804/69‑1 (MQ=255)    cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                                               >  1:2075395/1‑71 (MQ=255)    cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                                               >  1:3202839/1‑71 (MQ=255)    cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                                               >  1:3306986/1‑71 (MQ=255)     aGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacg                                                              >  1:1162447/1‑71 (MQ=255)     aGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacg                                                              >  1:1433351/1‑71 (MQ=255)      gTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacgg                                                             >  1:3550676/1‑71 (MQ=255)      gTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacgg                                                             >  1:3906498/1‑71 (MQ=255)          gTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacggcc                                                           <  1:1633531/69‑1 (MQ=255)            cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACAGCCACGGCCatat                                                       >  1:3121457/1‑71 (MQ=255)               tcatcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCa                                                    >  1:2098277/1‑71 (MQ=255)                catcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                                   >  1:3428461/1‑71 (MQ=255)                  tcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt                                                 <  1:2305776/71‑1 (MQ=255)                    aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                                    <  1:1682036/50‑1 (MQ=255)                    aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                                    <  1:618009/50‑1 (MQ=255)                    aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                               >  1:3538684/1‑71 (MQ=255)                    aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                               >  1:1033102/1‑71 (MQ=255)                     ggTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggcc                                                                 >  1:2380777/1‑52 (MQ=255)                     ggTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTgg                                               >  1:1976427/1‑70 (MQ=255)                      gTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCa                                             <  1:2268395/71‑1 (MQ=255)                        aGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                                           >  1:878512/1‑71 (MQ=255)                         gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCata                                                        >  1:3657247/1‑57 (MQ=255)                         gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                                  >  1:2225540/1‑63 (MQ=255)                         gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                                  >  1:3429532/1‑63 (MQ=255)                         gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGa                                           <  1:1835672/70‑1 (MQ=255)                         gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                          >  1:2999190/1‑71 (MQ=255)                         gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                          >  1:1429495/1‑71 (MQ=255)                          tCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATac                                                                     <  1:1141364/43‑1 (MQ=255)                          tCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATAcggccacggcca                                                          <  1:914329/54‑1 (MQ=255)                          tCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCata                                                        >  1:961492/1‑56 (MQ=255)                          ctgATGCCATACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTg                                                >  1:941580/3‑64 (MQ=255)                          ctgATGCCATACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTg                                                >  1:208337/3‑64 (MQ=255)                            gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                       >  1:1975028/1‑71 (MQ=255)                            gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                       >  1:772263/1‑71 (MQ=255)                            gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                       <  1:651234/71‑1 (MQ=255)                            gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCg                                        <  1:3947610/70‑1 (MQ=255)                              tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacggcc                                                           <  1:2008917/49‑1 (MQ=255)                              tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                                          >  1:1150230/1‑66 (MQ=255)                              aaCCGTACAGCGCACCGCGACGGTAGTCACCGATGATacggccacggcc                                                           <  1:3550818/48‑1 (MQ=255)                               aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                                       >  1:3559219/1‑68 (MQ=255)                               aCCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGt                                      >  1:2029271/1‑69 (MQ=255)                                ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGAATCTGGCAGACCGGTCa                                    <  1:2138738/70‑1 (MQ=255)                                 cGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGa                                  <  1:614699/71‑1 (MQ=255)                                  gTACAGCGCAACGCGACGGTAGTCACCGATGATacggcca                                                                >  1:3048815/1‑40 (MQ=255)                                  gTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                 <  1:2199452/71‑1 (MQ=255)                                     cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                 <  1:1353960/68‑1 (MQ=255)                                      aGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGcc                             <  1:3794751/71‑1 (MQ=255)                                         gctaCGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGa                          <  1:2344681/68‑1 (MQ=255)                                         gcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt                                                 <  1:1224188/48‑1 (MQ=255)                                            aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACAcc                             <  1:3513508/65‑1 (MQ=255)                                            aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAttt                       <  1:3367157/71‑1 (MQ=255)                                             cgcgACGGTAGTCACCGATGATACGGCCACGGCCata                                                        >  1:1296742/1‑37 (MQ=255)                                             cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                      <  1:3746778/71‑1 (MQ=255)                                             cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa                      >  1:1437964/1‑71 (MQ=255)                                                 aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACgg                   >  1:3381941/1‑70 (MQ=255)                                                  cGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc                  >  1:2347202/1‑70 (MQ=255)                                                   ggTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCaa                >  1:277157/1‑71 (MQ=255)                                                     tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg              <  1:3363619/71‑1 (MQ=255)                                                       gTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCg              >  1:1254477/1‑69 (MQ=255)                                                        tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                        <  1:2711657/58‑1 (MQ=255)                                                        tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                     <  1:3535560/61‑1 (MQ=255)                                                        tCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg           <  1:2118094/71‑1 (MQ=255)                                                          aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCa                                    >  1:2631320/1‑44 (MQ=255)                                                          aCCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAgc                                  <  1:3941157/46‑2 (MQ=255)                                                           ccGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt        <  1:1036060/71‑1 (MQ=255)                                                           ccGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGt                                      <  1:2146723/41‑1 (MQ=255)                                                            cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg           >  1:2603756/1‑67 (MQ=255)                                                            cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg           >  1:762822/1‑67 (MQ=255)                                                            cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg       >  1:1591625/1‑71 (MQ=255)                                                            cGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGAtt                        <  1:1316290/54‑1 (MQ=255)                                                             gatgatACGGCCACGGCCATATGGATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt      >  1:3618208/1‑71 (MQ=255)                                                             gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt      >  1:2490979/1‑71 (MQ=255)                                                             gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat      <  1:3795941/71‑3 (MQ=255)                                                              atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTc                                     >  1:1744646/1‑39 (MQ=255)                                                              atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTc                                     >  1:2591318/1‑39 (MQ=255)                                                              atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc             >  1:563614/1‑63 (MQ=255)                                                              atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc             >  1:1713932/1‑63 (MQ=255)                                                              atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc     >  1:8738/1‑71 (MQ=255)                                                              atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAAGAtag      >  1:1777735/1‑68 (MQ=255)                                                                 atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg  <  1:2110042/71‑1 (MQ=255)                                                                 atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg   <  1:264934/70‑1 (MQ=255)                                                                  tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg   >  1:1819363/1‑69 (MQ=255)                                                                  tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg   >  1:403082/1‑69 (MQ=255)                                                                   |                                                                     TGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGG  >  minE/609657‑609791

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