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

CACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGT  >  NZ_CP009273/3010931‑3011095                                                                                |                                                                                      cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                             >  2:271664/1‑90 (MQ=255) cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                             >  2:198793/1‑90 (MQ=255)   ccAGTGCAAATACAGCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                           <  1:83119/90‑1 (MQ=255)               cAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa                                                               <  1:38795/90‑1 (MQ=255)               cAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa                                                               <  1:461394/90‑1 (MQ=255)                aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                          <  1:459630/78‑1 (MQ=255)                aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                          >  2:459630/1‑78 (MQ=255)                   tGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGc                                                           <  2:285661/90‑1 (MQ=255)                    gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCt                                                          >  2:297087/1‑90 (MQ=255)                         cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                     >  2:521222/1‑90 (MQ=255)                                  tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                            <  1:397633/90‑1 (MQ=255)                                   aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                           >  2:627298/1‑90 (MQ=255)                                   aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                           >  1:285974/1‑90 (MQ=255)                                        aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                                      >  1:616952/1‑90 (MQ=255)                                        aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                                      >  1:438885/1‑90 (MQ=255)                                         gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGtt                                                                               >  1:616542/1‑48 (MQ=255)                                         gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGtt                                                                               <  2:616542/48‑1 (MQ=255)                                         gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                     >  2:182679/1‑90 (MQ=255)                                         gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                     >  1:564983/1‑90 (MQ=255)                                                 aaTATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGc                             <  2:564983/90‑1 (MQ=255)                                                     tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc                         >  2:598208/1‑90 (MQ=255)                                                        gcgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCgg                      >  1:488910/1‑90 (MQ=255)                                                         cgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGt                     <  1:297087/90‑1 (MQ=255)                                                              tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg                <  2:285974/90‑1 (MQ=255)                                                               gAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGt               >  1:493928/1‑90 (MQ=255)                                                                   aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt           <  1:478501/90‑1 (MQ=255)                                                                        cgcgCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGc      <  1:474185/90‑1 (MQ=255)                                                                            cTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGt  <  2:28330/90‑1 (MQ=255)                                                                                |                                                                                      CACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGCAGGT  >  NZ_CP009273/3010931‑3011095

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