Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,705	A→G	5.5%	R183R (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,705	0	A	G	5.5%	205.5 / 2.6	73	R183R (CGT→CGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (37/32);  new base G (2/2);  total (39/34)
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.75e-01

CAGCCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGG  >  minE/609635‑609774                                                                       |                                                                      cAGCCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGa                                                                       >  1:3330366/1‑71 (MQ=255) cAGCCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGa                                                                       >  1:3288259/1‑71 (MQ=255)    ccTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACgg                                                                    <  1:62663/71‑1 (MQ=255)     cTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACgg                                                                    <  1:3103040/70‑1 (MQ=255)     cTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGt                                                                   <  1:2508675/71‑1 (MQ=255)     cTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGt                                                                   <  1:2508047/71‑1 (MQ=255)        cAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGt                                                                >  1:2879555/1‑71 (MQ=255)        cAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGt                                                                >  1:2624355/1‑71 (MQ=255)             aaGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCAcc                                                            >  1:1064537/1‑70 (MQ=255)                tGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                        <  1:2263731/71‑1 (MQ=255)                 gAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTc                                                               <  1:3950375/63‑1 (MQ=255)                 gAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                       >  1:2835699/1‑71 (MQ=255)                   aCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                       <  1:1435398/69‑1 (MQ=255)                    cTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCa                                                              >  1:2789970/1‑61 (MQ=255)                     tgtgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACgg                                                                    <  1:3229497/54‑1 (MQ=255)                       tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCTCAACGCGACGGTAGTCACCGATGATacgg                                                  <  1:3307852/70‑1 (MQ=255)                       tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacgg                                                  <  1:2629354/70‑1 (MQ=255)                       tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                   <  1:404804/69‑1 (MQ=255)                          cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAg                                                                 >  1:2959841/1‑52 (MQ=255)                          cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGt                                                                >  1:3442974/1‑53 (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)                             tttGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGa                                                                       <  1:3752511/43‑1 (MQ=255)                             tttGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGa                                                                       <  1:839794/43‑1 (MQ=255)                                gTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacggcc                                          <  1:1633531/69‑1 (MQ=255)                                  cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCa                                                              <  1:3920702/47‑1 (MQ=255)                                  cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACAGCCACGGCCatat                                      >  1:3121457/1‑71 (MQ=255)                                   tttCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                        >  1:1992020/1‑52 (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:618009/50‑1 (MQ=255)                                          aGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                   <  1:1682036/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:3429532/1‑63 (MQ=255)                                               gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                 >  1:2225540/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:1437964/1‑71 (MQ=255)                                                                   cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa     <  1:3746778/71‑1 (MQ=255)                                                                       aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACgg  >  1:3381941/1‑70 (MQ=255)                                                                       |                                                                      CAGCCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGG  >  minE/609635‑609774

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