Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,011,010	A→G	E531G (GAA→GGA)	ygfK →	putative selenate reductase subunit YgfK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,011,010	0	A	G	100.0%	102.9 / NA	32	E531G (GAA→GGA)	ygfK	putative selenate reductase subunit YgfK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/11);  total (21/11)

TTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATT  >  NZ_CP009273/3010923‑3011086                                                                                        |                                                                             tttGCAACCCCCAGTCCAAAAACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGgaa                                                                            <  2:424839/90‑1 (MQ=255)         cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                    >  2:418493/1‑90 (MQ=255)           ccAGTCCAAATACACCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                  <  2:581036/90‑1 (MQ=255)            cAGTGCAAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                 <  2:53272/90‑1 (MQ=255)              gTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGa                                                               >  2:144264/1‑90 (MQ=255)                        aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCaa                                                     >  1:173516/1‑90 (MQ=255)                            gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCt                                                 >  2:189203/1‑90 (MQ=255)                              aCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTgg                                               >  2:539374/1‑90 (MQ=255)                              aCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTgg                                               >  2:399950/1‑90 (MQ=255)                              aCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTgg                                               >  2:290870/1‑90 (MQ=255)                                 cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                            >  1:517937/1‑90 (MQ=255)                                           aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                  <  2:79125/90‑1 (MQ=255)                                           aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                  >  2:286747/1‑90 (MQ=255)                                               cAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTc                              <  2:390258/90‑1 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATg                                                                >  1:260802/1‑54 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATg                                                                <  2:260802/54‑1 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:14933/1‑90 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:460474/1‑90 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:307414/1‑90 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:314966/1‑90 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:293453/1‑90 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:294115/1‑90 (MQ=255)                                                 gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGATGGCACAAACCAGCCGGTTCTg                            >  1:287663/1‑90 (MQ=255)                                                           tatCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                  <  1:16170/42‑1 (MQ=255)                                                           tatCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                  >  2:16170/1‑42 (MQ=255)                                                            atCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCAt                 <  2:315928/90‑1 (MQ=255)                                                                gcgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCgg             >  2:241423/1‑90 (MQ=255)                                                                 cgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGt            <  1:387139/90‑1 (MQ=255)                                                                 cgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGt            <  1:241423/90‑1 (MQ=255)                                                                  gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt           >  1:561836/1‑90 (MQ=255)                                                                  gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt           >  1:23940/1‑90 (MQ=255)                                                                          aaaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAt   <  2:561836/90‑1 (MQ=255)                                                                           aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt  <  2:319793/90‑1 (MQ=255)                                                                                        |                                                                             TTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATT  >  NZ_CP009273/3010923‑3011086

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑