Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,220,044	A→G	S236P (TCG→CCG)	yohF ←	SDR family oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,220,044	0	A	G	100.0%	61.9 / NA	19	S236P (TCG→CCG)	yohF	SDR family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/4);  total (15/4)

ACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGC  >  NZ_CP009273/2219957‑2220126                                                                                        |                                                                                   aCCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCt                                                                                  >  1:788992/1‑90 (MQ=255)     aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                              <  2:368212/90‑1 (MQ=255)      aTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGGGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCg                                                                             >  1:322285/1‑90 (MQ=255)       tGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                            <  1:815256/90‑1 (MQ=255)         aaaTATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCggtg                                                                          >  1:396616/1‑90 (MQ=255)                cACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAAt                                                                   >  1:387379/1‑90 (MQ=255)                               gcTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                                    <  2:397634/90‑1 (MQ=255)                                            gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcg                                                              >  1:153658/1‑67 (MQ=255)                                            gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcg                                                              <  2:153658/67‑1 (MQ=255)                                            gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                       >  1:571410/1‑90 (MQ=255)                                            gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                       >  1:77881/1‑90 (MQ=255)                                            gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                       >  1:122229/1‑90 (MQ=255)                                            gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                       >  1:438878/1‑90 (MQ=255)                                            gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                       >  1:288860/1‑90 (MQ=255)                                            gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                       >  1:13767/1‑90 (MQ=255)                                                 gtgGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                                  >  2:281494/1‑90 (MQ=255)                                                                       gccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGcc            >  2:78341/1‑90 (MQ=255)                                                                       gccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGcc            >  1:575056/1‑90 (MQ=255)                                                                                 aTCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAAcgccgc  >  1:880173/1‑90 (MQ=255)                                                                                        |                                                                                   ACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGC  >  NZ_CP009273/2219957‑2220126

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