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%	111.0 / NA	34	E531G (GAA→GGA)	ygfK	putative selenate reductase subunit YgfK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/7);  total (27/7)

ATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATT  >  NZ_CP009273/3010929‑3011086                                                                                  |                                                                             aTCACAAGTCCAAATAAAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGtt                                                                      <  2:495734/90‑1 (MQ=255)        gTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGa                                                               >  2:514354/1‑90 (MQ=255)             aaTACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAAtata                                                          >  2:271186/1‑90 (MQ=255)                  aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCaa                                                     >  2:646221/1‑90 (MQ=255)                      gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCt                                                 >  1:206103/1‑90 (MQ=255)                        aCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTgg                                               >  1:324066/1‑90 (MQ=255)                           cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                            >  2:509774/1‑90 (MQ=255)                              cTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACaaa                                         >  2:123350/1‑90 (MQ=255)                                ggATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAAcc                                       >  1:709228/1‑90 (MQ=255)                                    tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                   <  2:488958/90‑1 (MQ=255)                                     aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                  >  1:548291/1‑90 (MQ=255)                                     aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                  >  1:714260/1‑90 (MQ=255)                                          aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                             >  2:713201/1‑90 (MQ=255)                                           gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:625585/1‑90 (MQ=255)                                           gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:477677/1‑90 (MQ=255)                                           gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:428292/1‑90 (MQ=255)                                           gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:106527/1‑90 (MQ=255)                                           gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:375624/1‑90 (MQ=255)                                           gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  2:297779/1‑90 (MQ=255)                                           gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:245967/1‑90 (MQ=255)                                           gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGATGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                            >  1:76537/1‑90 (MQ=255)                                                 tGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCAc                      <  1:486709/90‑1 (MQ=255)                                                 tGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGATGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCAc                      <  2:76537/90‑1 (MQ=255)                                                  gAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                   <  2:135787/76‑1 (MQ=255)                                                  gAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                   >  1:135787/1‑76 (MQ=255)                                                      atCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCAt                 >  2:124022/1‑90 (MQ=255)                                                       tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc                >  2:546265/1‑90 (MQ=255)                                                       tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc                >  1:189811/1‑90 (MQ=255)                                                           cgAACTGAAAAAAGTCGCGCTGGGGAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGt            <  2:95224/90‑1 (MQ=255)                                                            gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt           >  2:393824/1‑90 (MQ=255)                                                            gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt           >  1:369148/1‑90 (MQ=255)                                                            gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt           >  2:708820/1‑90 (MQ=255)                                                            gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt           >  1:104239/1‑90 (MQ=255)                                                                     aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGCTTCTGGTTCACGCCATCCGGTTGCCGTGAtt  <  2:548291/90‑1 (MQ=255)                                                                                  |                                                                             ATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATT  >  NZ_CP009273/3010929‑3011086

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