Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,126,699	Δ1 bp	coding (335/879 nt)	pflC	pyruvate formate lyase II activase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	4,126,699	0	T→.	100.0%	48.5	16	coding (335/879 nt)	pflC	pyruvate formate lyase II activase
Reads supporting (aligned to +/- strand):  new base (6/10):  ref base (0/0):  total (6/10)

TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCACG  >  REL606/4126664‑4126730                                    |                                tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                 <  1:75775/36‑1 (MQ=255) tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                 <  1:2693199/36‑1 (MQ=255) tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                 <  1:4811343/36‑1 (MQ=255)   aCGCTCTGGAACGGGAAGTGATGAAAGATGACAttt                               <  1:755321/36‑1 (MQ=255)    cGCTCTGGAACGGGAAGGGATGAAAGATGACAtttt                              <  1:7104088/36‑1 (MQ=255)     gCTCTGGAACGGGAAGTGATGAAAGATGACAttttt                             <  1:313383/36‑1 (MQ=255)     gCTCTAGAACGGGAAGTGATGAAAGATGACAttttt                             <  1:3315030/36‑1 (MQ=255)      ctctGGAACGGGAAGTGATGAAAGATGACAtttttt                            <  1:6171431/36‑1 (MQ=255)       tctGGAACGGGAAGTGATGAAAGATGACAttttttt                           >  1:5898541/1‑36 (MQ=255)        ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                          <  1:4131841/36‑2 (MQ=255)         tGGAACGGGAAGTGATGAAAGATGACAtttttttcg                         <  1:4860010/36‑3 (MQ=255)          ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        >  1:4001711/1‑33 (MQ=255)          ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        <  1:1745955/36‑4 (MQ=255)          ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        <  1:6025660/36‑4 (MQ=255)          ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        <  1:5334793/36‑4 (MQ=255)            aaCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCAc                     <  1:5222821/36‑1 (MQ=25)              cGGGAAGTGATGAAAGATGACA.TTTTTTTCGCACGt                   <  1:3082321/36‑1 (MQ=25)               gggAAGTGATGATAGATGACA.TTTTTTTCGCACGTc                  <  1:3090296/36‑1 (MQ=21)                ggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTcc                 <  1:1567857/36‑1 (MQ=25)                  aaGTGATGAAAGATGACA.TTTTTTTCGCACGTCcgg               <  1:6817167/36‑1 (MQ=25)                  aaGTGATGAAAGATGACA.TTTTTTTCGCACGTCcgg               >  1:3827949/1‑36 (MQ=25)                  aaGTGATGAAAGATGACA.TTTTTATCGCACCTCcgg               >  1:2994927/1‑36 (MQ=18)                   aGTGATGAAAGATGACA.TTTTTTTCGCACGTCcggc              >  1:6164439/1‑36 (MQ=25)                    gTGATGAAAGATGACA.TTTTTTTCGCACGTCcggcg             <  1:3557922/36‑1 (MQ=25)                     tgatgaAAGATGACA.TTTTTTTCGCACGTCcggcgg            <  1:6555534/36‑1 (MQ=25)                      gatgaAAGATGACA.TTTTTTTCGCACGTCcggcggc           <  1:7491051/36‑1 (MQ=25)                       atgaAAGATGACA.TTTTTTTCGCACGTCcgccggcg          >  1:3437022/1‑36 (MQ=21)                            aGATGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTc     <  1:6263417/36‑1 (MQ=25)                             gATGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCa    >  1:3676875/1‑36 (MQ=25)                              aTGACA.TTTTTTTCGCACGTCCGGCGGAGGCGTCaa   >  1:550706/1‑35 (MQ=21)                               tGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCACg  <  1:5384542/36‑1 (MQ=25)                                    |                                TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCACG  >  REL606/4126664‑4126730 Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG