Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,499,566	A→G	S258P (TCC→CCC)	ydcO ←	BenE family transporter YdcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,499,566	0	A	G	100.0%	88.5 / NA	28	S258P (TCC→CCC)	ydcO	BenE family transporter YdcO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/12);  total (16/12)

TGATCTTTATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGAAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATTGCTGCGA  >  NZ_CP009273/1499481‑1499646                                                                                      |                                                                                 tGATCTTTATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGgaaaa                                                                              >  1:1193086/1‑90 (MQ=255)       ttATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGt                                                                        >  2:612821/1‑90 (MQ=255)       ttATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGt                                                                        >  2:105252/1‑90 (MQ=255)         aTCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGc                                                                      >  2:514974/1‑90 (MQ=255)             ggATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGc                                                                  <  1:870605/90‑1 (MQ=255)              gATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGGGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCa                                                                 >  1:1404469/1‑90 (MQ=255)               aTGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCaa                                                                <  2:70217/90‑1 (MQ=255)                 gcgcTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATc                                                              <  1:1471768/90‑1 (MQ=255)                 gcgcTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATc                                                              <  1:786106/90‑1 (MQ=255)                 gcgcTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATc                                                              <  1:484329/90‑1 (MQ=255)                   gcTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCa                                                            >  2:865837/1‑90 (MQ=255)                      tCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTa                                                         <  1:373289/90‑1 (MQ=255)                          ggCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATa                                                     <  1:491394/90‑1 (MQ=255)                              ttGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAAt                                                 >  1:469913/1‑90 (MQ=255)                              ttGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAAt                                                 >  1:1309402/1‑90 (MQ=255)                                        gCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGaa                                       >  2:1410201/1‑90 (MQ=255)                                          cgcgGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAAc                                     >  2:1254005/1‑90 (MQ=255)                                          cgcgGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAAc                                     >  2:1343462/1‑90 (MQ=255)                                          cgcgGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAAc                                     >  2:712107/1‑90 (MQ=255)                                               tGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAg                                <  1:102696/90‑1 (MQ=255)                                                gATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGc                               <  2:357908/90‑1 (MQ=255)                                                             ccGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAgct                  <  2:1309402/90‑1 (MQ=255)                                                               gACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAgctgc                >  1:1400225/1‑90 (MQ=255)                                                                    aaTAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCa           >  1:456674/1‑90 (MQ=255)                                                                    aaTAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCa           >  1:1178694/1‑90 (MQ=255)                                                                      tAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGCAACAGGAGCCGAATATCCAGCTGCTTTCAtt         <  2:1193086/90‑1 (MQ=255)                                                                        aaCGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATtgc       <  1:466083/90‑1 (MQ=255)                                                                             cGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATTGCTGCGa  >  1:300009/1‑90 (MQ=255)                                                                             cGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATTGCTGCGa  >  1:207134/1‑90 (MQ=255)                                                                                      |                                                                                 TGATCTTTATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGAAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATTGCTGCGA  >  NZ_CP009273/1499481‑1499646

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