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%	98.6 / NA	28	T95A (ACC→GCC)	cysK	cysteine synthase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/12);  total (16/12)

GAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGC  >  NZ_CP009273/2525966‑2526135                                                                                     |                                                                                      gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  2:724617/90‑1 (MQ=255) gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  2:573448/90‑1 (MQ=255) gAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  1:38648/90‑1 (MQ=255)       aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTAtt                                                                            >  2:113473/1‑90 (MQ=255)          aGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGaa                                                                         <  1:868859/90‑1 (MQ=255)                aaTACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAAcgccgc                                                                   >  1:369439/1‑90 (MQ=255)                 aTACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCa                                                                  <  1:208631/90‑1 (MQ=255)                     cGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct                                                              <  2:776799/90‑1 (MQ=255)                     cGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgct                                                              >  1:271822/1‑90 (MQ=255)                        gATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctgct                                                           >  1:475243/1‑90 (MQ=255)                        gATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctgct                                                           >  2:249357/1‑90 (MQ=255)                           tGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGaa                                                        <  1:214565/90‑1 (MQ=255)                                 ggCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGGATTGAACGCCGCAAGCTGCTGAAAGCGtt                                                  >  1:734260/1‑90 (MQ=255)                                  gCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTa                                                 <  1:427782/90‑1 (MQ=255)                                      aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg                                             >  1:822238/1‑90 (MQ=255)                                      aTGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg                                             >  1:164188/1‑90 (MQ=255)                                           gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  2:870655/1‑90 (MQ=255)                                           gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  2:497021/1‑90 (MQ=255)                                           gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  2:412384/1‑90 (MQ=255)                                           gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  2:305396/1‑90 (MQ=255)                                                cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt                                   >  2:490109/1‑90 (MQ=255)                                                cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt                                   >  2:537525/1‑90 (MQ=255)                                                          tccAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa                         <  2:16475/88‑1 (MQ=255)                                                                  cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 >  1:498160/1‑90 (MQ=255)                                                                             tGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAg      >  2:779660/1‑90 (MQ=255)                                                                                 aGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGgcgc  <  1:296808/90‑1 (MQ=255)                                                                                 aGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGgcgc  <  1:599477/90‑1 (MQ=255)                                                                                 aGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGgcgc  <  2:822238/90‑1 (MQ=255)                                                                                     |                                                                                      GAACCGACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGC  >  NZ_CP009273/2525966‑2526135

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: ‑