Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,490,112	C→T	A105T (GCG→ACG)	glpE	thiosulfate sulfurtransferase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	3,490,112	0	C→T	100.0%	164.0	54	A105T (GCG→ACG)	glpE	thiosulfate sulfurtransferase
Reads supporting (aligned to +/- strand):  new base (28/24):  ref base (0/1):  total (29/25)

AAGGGGACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGTCGCTGCCAGGCTT  >  REL606/3490079‑3490144                                  |                                 aaGGGGACAGTATAAAGCGTTACGCGCCGTACGTCa                                >  1:1650378/1‑36 (MQ=255) aaGCGGACAGTATAAAGCGTTACGCGCCGTACGTCa                                >  1:4227837/1‑36 (MQ=25)  aGGGGACAGTATAAAGCGTTACGCGCCGTACGTCAc                               >  1:1867902/1‑36 (MQ=255)    gggACAGTATAAAGCGTTACGCGCCGTACGTCACct                             <  1:745871/36‑1 (MQ=255)     ggACAGTATAAAGCGTTACGCGCCGTACGTCACctc                            >  1:6717639/1‑36 (MQ=255)       aCAGTATAAAGCGTTACGCGCCGTACGTCACCTCTg                          <  1:3495867/36‑1 (MQ=255)       aCAGTATAAAGCGTTACGCGCCGTACGTCACCTCTg                          <  1:5217823/36‑1 (MQ=255)        cAGTATAAAGCGTTACGCGCCGTACGTCACCTCTGc                         >  1:1326779/1‑36 (MQ=255)         aGTATAAAGCGTTACGCGCCGTACGTCACCTCTGCg                        <  1:6414380/36‑1 (MQ=255)         aGTATAAAGCGTTACGCGCCGTACGTCACCTCTGCg                        >  1:1969532/1‑36 (MQ=255)           tataAAGCGTTACGCGCCGTACGTCACCTCTGCggg                      <  1:830706/36‑1 (MQ=255)             taAAGCGTTACGCGCCGTACGTCACCTCTGCGGGaa                    <  1:5341611/36‑1 (MQ=255)             taAAGCGTTACGCGCCGTACGTCACCTCTGCGGGaa                    <  1:2057461/36‑1 (MQ=255)             taAAGCGTTACGCGCCGTACGTCACCTCTGCGGGaa                    <  1:5072387/36‑1 (MQ=255)              aaaGCGTTATGCGCCGTACGTCACCTCTGCGGGaaa                   >  1:3040800/1‑36 (MQ=25)              aaaGCGTTACGCGCCGTACGTCACCTCTGCGGGaaa                   >  1:5597295/1‑36 (MQ=255)              aaaCCGCTACGCGCCGTACATCACCTCTGCGGGaca                   <  1:5397547/36‑3 (MQ=17)               aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                  >  1:5548933/1‑36 (MQ=255)               aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                  <  1:6599176/36‑1 (MQ=255)               aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                  >  1:1143469/1‑36 (MQ=255)               aaGCGTTACGCGCCGTACGTCACCTCTGCGGGAAAc                  >  1:1473750/1‑36 (MQ=255)                aGCGTGACGCGCCGTACGTCACCTCTGCGGGAAACt                 <  1:3315203/36‑1 (MQ=25)                 gCGTTACGCGCCGTACGTCACCTCTGCGGGGAACTg                >  1:6108374/1‑36 (MQ=25)                 gCGTTACGCGCCGTACGTCACCTCTGCGGGAAACTg                <  1:7465959/36‑1 (MQ=255)                 gCGTTACGCGCCGTACGTCACCTCTGCGGGAAACTg                >  1:281499/1‑36 (MQ=255)                  cGTTTCGCGCCGTACGTCAGCTCTGCGGGAAACTGt               <  1:2678739/36‑1 (MQ=21)                   gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              >  1:418318/1‑36 (MQ=255)                   gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              >  1:2943615/1‑36 (MQ=255)                   gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              <  1:4419835/36‑1 (MQ=255)                   gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              >  1:3663757/1‑36 (MQ=255)                   gTTACGCGCCGTACGTCACCTCTGCGGGAAACTGTc              >  1:5493009/1‑36 (MQ=255)                    ttACGCGCCGTACGTCTCCTCTGCGGGAAACTGTCg             <  1:6711266/36‑1 (MQ=25)                    ttACGCGCCGTACGTCACCTCTGCGGGAAACTGTCg             >  1:673623/1‑36 (MQ=255)                    ttACGCGCCGTACGTCACCTCTGCGGGAAACTGTCg             >  1:5185007/1‑36 (MQ=255)                     tACGCGCCGTACGTCACCTCTGCGGGAAACTGTCGc            >  1:4389826/1‑36 (MQ=255)                      aCGCGCCGTACGTCACCTCTGCGGGAAACTGTCGCt           <  1:1881160/36‑1 (MQ=255)                      aCGCGCCGTACGTCACCTCTGCGGGAAACTGTCGCt           >  1:3748736/1‑36 (MQ=255)                      aCGCGCCGTACGGCACCTCTGCGGGAAACTGTCGCt           >  1:7323983/1‑36 (MQ=255)                      aCGCCCCGTACGTCACCTCTGCGGGAAACTGTCGCt           <  1:2439653/36‑1 (MQ=25)                       cgcgCCGTACGTCACCTCTGCGGGAAACTGTCGCTg          <  1:2704761/36‑1 (MQ=255)                         cgcCGTACGTCACCTCTGCGGGAAACTGTCGCTGcc        <  1:1710121/36‑1 (MQ=255)                          gcCGTACGTCGCTTCTGCGGGAAACTGTCGCTGCCa       <  1:529802/36‑1 (MQ=21)                          gcCGTACGTCACCTCTGCGGGAAACTGTTGCTGCCa       >  1:5986895/1‑36 (MQ=25)                          gcCGTACGTCACCTCTGCGGGAAACTGTCGCTGctt       >  1:660059/1‑34 (MQ=25)                          gcCGTACGTCACCTCTGCGGGAAACTGTCGCTGCCa       >  1:6354661/1‑36 (MQ=255)                           ccGTACGTCACCTCTGCGGGAAACTGTCGCTGCCAg      >  1:1360071/1‑36 (MQ=255)                           ccGTACGTCACCTCTGCGGGAAACTGTCGCTACCAg      <  1:4821406/36‑1 (MQ=25)                            cGTACGTCACCTCTGCGGGAAACTGTCGCTGCCAgg     <  1:6862634/36‑1 (MQ=255)                             gTACGTCACCTCTGCGGGAAACTGTCGCTGCCAGGc    <  1:606757/36‑1 (MQ=255)                             gTACGTCACCTCTGCGGGAAACTGTCGCTGCCAGGc    >  1:492415/1‑36 (MQ=255)                             gTACGTCACCTCTGCGGGAAACTGTCGCTGCCAGGc    >  1:4581802/1‑36 (MQ=255)                              tACGTCACCTCTGCGGGAAACTGTCGCTGCCAAGCt   >  1:4667512/1‑36 (MQ=25)                              tACGCCACCCCTGCGGGAAACTGTCGCTGCCAGGCt   <  1:2247892/36‑1 (MQ=255)                               aCGTCACCTCTGCGGGAAACTGTCGCTGCCAGGCtt  <  1:2285899/36‑1 (MQ=255)                                  |                                 AAGGGGACAGTATAAAGCGTTACGCGCCGTACGCCACCTCTGCGGGAAACTGTCGCTGCCAGGCTT  >  REL606/3490079‑3490144 Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG