Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	2,501,962	0	T→G	23.1%	10.7	39	V500G (GTG→GGG)	tktB	transketolase 2, thiamin‑binding
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (3/27):  total (12/27)

CCTGGCGACCGTGCGATCAGGTGGAAGCGGCGGTGGGCTGGAAGCTGGCGGTTGAGCGCCACAACGGAC  >  REL606/2501929‑2501997                                  |                                    cgtGGCGACCGTGCGATCAGGTGGAAGCGGCGGTgg                                   <  1:18096/34‑1 (MQ=255) ccTGGCGACCGTGCGATCAGGTGGAAGCGGCGGTgg                                   <  1:1807346/36‑1 (MQ=255) ccTGGCGACCGTGCGATCAGGTGGAAGCGGCGGGgg                                   >  1:5745819/1‑36 (MQ=255)  cTGGCGACCGTGCGATCAGGTGGAAGCGGCGGGggg                                  >  1:5983135/1‑36 (MQ=255)  cTGGCGACCGTGCGATCAGGTGGAAGCGGCGGGggg                                  >  1:6955437/1‑36 (MQ=255)  cTGGCGACCGTGCCATCAGGTGGAAGCGGCGGGggg                                  >  1:1805536/1‑36 (MQ=25)   tGGCGACCGTGCGATCAGGTGGAAGCGGCGGGgggg                                 >  1:4049194/1‑35 (MQ=255)    ggCGACCGTGCGATCAGGTGGAAGCGGCGGTGGGCt                                <  1:3701527/36‑1 (MQ=255)    ggCGACCGTGCGATCAGGTGGAAGCGGCGGGggggg                                >  1:4379753/1‑34 (MQ=37)      cGACCGTGCGATCAGGTGGAAGCGGCGGTGGGCTgg                              <  1:2332210/36‑1 (MQ=255)       gACCGTGCGATCAGGTGGAAGCGGCGGTGGGCTGGa                             <  1:5612027/36‑1 (MQ=255)       gACCGTGCGATCAGGTGGAAGCGGCGGTGGGCTGGa                             <  1:3808627/36‑1 (MQ=255)       gACCGTGCGATCAGGTGGAAGCGGCGGTGGGCTGGa                             <  1:1691946/36‑1 (MQ=255)        aCCGTGGGATCAGGTGGAAGCGGCGGTGGGCTGGaa                            <  1:6116563/36‑1 (MQ=255)        aCCGTGCGATCAGGTGGAAGCGGCGGGGGGGTgggg                            >  1:1646191/1‑34 (MQ=25)         ccGTGCGATCAGGTGGAAGCGGCGGTGGGCTGGAAg                           <  1:1730669/36‑1 (MQ=255)         ccGTGCGATCAGGTGGAAGCGGCGGGGGGCGggggg                           >  1:4779344/1‑33 (MQ=21)          cGTGCGATGAGGTGGAAGCGGCGGTGGGCTGGAAGc                          <  1:2626837/36‑1 (MQ=255)          cGTGCGATCAGGTGGAAGCGGCGGGGGGCGGGGAgg                          >  1:5008782/1‑35 (MQ=21)           gTGCGATCAGGTGGAAGCGGCGGTGGGCTGGAAGCt                         <  1:6747257/36‑1 (MQ=255)           gTGCGATCAGGTGGAAGCGGCGGTGGGCTGGAAGCt                         <  1:7296342/36‑1 (MQ=255)            tGCGATCAGGTGGAAGCGGCGGTGGGCTGGAAGCTg                        <  1:889813/36‑1 (MQ=255)             gCGATCAGGTGGAAGCGGCGGTGGGCTGGAAGCTgg                       <  1:6297215/36‑1 (MQ=255)              cGATCAGGTGGAAGCGGCGGTGGGCGGGAAGCTGGc                      <  1:5839349/36‑1 (MQ=255)               gATCAGGTGGAAGCGGCGGTGGGCTGGAAGCTGGCg                     <  1:1276088/36‑1 (MQ=255)                aTCAGGAGGAAGCGGCGGTGGGCTGGAAGCTGGCgg                    <  1:6354625/36‑1 (MQ=255)                 tCAGGTGGAAGCGGCGGTGGGCGGGAAGCTGGCGGt                   <  1:6766297/36‑1 (MQ=255)                  ctggTGGAAGCGGCGGTGGGCTGGAAGCTGGCGGtt                  <  1:1712720/34‑1 (MQ=255)                     gTGGGAGCGGCGGTGGGCTGGAAGCTGGCGGTTGAg               <  1:7197159/36‑1 (MQ=255)                       gaaaGCGGCGGTGGGCTGGAAGCTGGCGGTTGAgcg             <  1:1105947/34‑1 (MQ=255)                        gAAGCGGCGGTGGGCTGGAAGCTGGCGGTTGAgcgc            <  1:734679/36‑1 (MQ=255)                           gcggcgGTGGGCTGGAAGCTGGCGGTTGAGCGCcac         <  1:5363746/36‑1 (MQ=255)                           gcggcgGTGGGCTGGAAGCTGGCGGTTGAGCGCcac         <  1:257700/36‑1 (MQ=255)                           gcgccgCTGTGCTGGAAGCTGGCGGTTGAGCGCcac         <  1:6550039/36‑1 (MQ=18)                             ggcggTGGGCTGGAAGCTGGCGGTTGAGCGCCACaa       <  1:7028983/36‑1 (MQ=255)                               cggTGGGCTGGAAGCTGGCGGGTCAACGCCCCAAcc     >  1:677422/1‑35 (MQ=17)                                 gTGGGCTGGAAGCTGGCGGTTGAGCGCCACAACGGa   <  1:425116/36‑1 (MQ=255)                                 gTGGGCTGGAAGCTGGCGGGTGAGCGCCCCCACgcc   >  1:1304315/1‑34 (MQ=18)                                 gTGGGCTGGAAGCTGGCGGGGGGGCGGCGCAACggg   >  1:5798857/1‑35 (MQ=12)                                  tGGGCTGGAAGCGGGCGGTTGAGCGCCACAACggtc  <  1:6457285/36‑3 (MQ=25)                                  |                                    CCTGGCGACCGTGCGATCAGGTGGAAGCGGCGGTGGGCTGGAAGCTGGCGGTTGAGCGCCACAACGGAC  >  REL606/2501929‑2501997 Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG