Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,526,050	A→G	T95A (ACC→GCC)	cysK →	cysteine synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,526,050	0	A	G	100.0%	68.7 / NA	20	T95A (ACC→GCC)	cysK	cysteine synthase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/10);  total (10/10)

TTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGG  >  NZ_CP009273/2525964‑2526132                                                                                       |                                                                                   ttGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCa                                                                                 <  2:524353/90‑1 (MQ=255)   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  1:310203/90‑1 (MQ=255)   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  2:4290/90‑1 (MQ=255)   gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                               <  1:661530/90‑1 (MQ=255)    aaCCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGa                                                                              >  1:589120/1‑90 (MQ=255)                            ttGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGa                                                      <  1:26258/90‑1 (MQ=255)                            ttGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGa                                                      <  2:763091/90‑1 (MQ=255)                            ttGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGa                                                      <  1:129036/90‑1 (MQ=255)                                        aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg                                          >  1:738878/1‑90 (MQ=255)                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  1:602572/1‑90 (MQ=255)                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  2:314832/1‑90 (MQ=255)                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  1:462885/1‑90 (MQ=255)                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  2:664028/1‑90 (MQ=255)                                             gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAAAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                     >  2:586537/1‑90 (MQ=255)                                                     tCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAGCGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCt                             >  1:404917/1‑90 (MQ=255)                                                            tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa                      <  2:175876/90‑1 (MQ=255)                                                             aCAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAg                     <  2:534132/90‑1 (MQ=255)                                                                          gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAt        >  1:509734/1‑90 (MQ=255)                                                                             cATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGaa     <  2:602572/90‑1 (MQ=255)                                                                                gCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAgg  >  2:369656/1‑90 (MQ=255)                                                                                       |                                                                                   TTGAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGG  >  NZ_CP009273/2525964‑2526132

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