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

GACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGC  >  NZ_CP009273/2525971‑2526135                                                                                |                                                                                      gACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTAt                                                                             >  2:433954/1‑90 (MQ=255)  aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  >  1:401335/1‑84 (MQ=255)  aCCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  2:401335/84‑1 (MQ=255)             tACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCaa                                                                 <  2:26837/90‑1 (MQ=255)                 gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  >  1:229665/1‑69 (MQ=255)                 gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  2:229665/69‑1 (MQ=255)                 gggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctg                                                             >  1:159974/1‑90 (MQ=255)                  ggATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAGCGCCGCAAgctgc                                                            >  1:509672/1‑90 (MQ=255)                   gATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAgctgct                                                           >  1:570256/1‑90 (MQ=255)                               cTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  <  2:531624/55‑1 (MQ=255)                               cTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATg                                                                                  >  1:531624/1‑55 (MQ=255)                                      gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTGTTGAACGCCGCGAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  2:237191/1‑90 (MQ=255)                                      gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  1:306078/1‑90 (MQ=255)                                      gCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                        >  1:119124/1‑90 (MQ=255)                                        tGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGaa                                                                         >  1:503322/1‑55 (MQ=255)                                        tGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGaa                                                                         <  2:503322/55‑1 (MQ=255)                                         gCCGCTCGCGGTTACAAATTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTg                                     >  1:252187/1‑90 (MQ=255)                                           cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt                                   >  1:429402/1‑90 (MQ=255)                                           cGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGt                                   <  2:306078/90‑1 (MQ=255)                                                     tACAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGaa                         <  2:252187/90‑1 (MQ=255)                                                      aCAAACTCACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAg                        <  2:570256/90‑1 (MQ=255)                                                             cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 >  1:450580/1‑90 (MQ=255)                                                             cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 <  1:376945/90‑1 (MQ=255)                                                             cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 <  2:54450/90‑1 (MQ=255)                                                             cACCCTGACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTaa                 >  1:314768/1‑90 (MQ=255)                                                                   gACCATGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCAt           >  1:547497/1‑90 (MQ=255)                                                                        tGCCAGAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAg      >  1:531951/1‑90 (MQ=255)                                                                             gAAGCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGgcgc  <  2:520228/89‑1 (MQ=255)                                                                                |                                                                                      GACCAGCGGTAATACCGGGATTGCACTGGCCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGCGC  >  NZ_CP009273/2525971‑2526135

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