Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,286,033	G→A	G186S (GGC→AGC)	agaD →	galactosamine‑specific PTS enzyme IID component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,286,033	0	G	A	100.0%	89.7 / NA	26	G186S (GGC→AGC)	agaD	galactosamine‑specific PTS enzyme IID component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/14);  total (12/14)

GTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCGCT  >  NC_000913/3285955‑3286109                                                                               |                                                                             gTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATcagc                                                                            <  1:553312/81‑1 (MQ=255) gTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATcagc                                                                            <  2:553312/81‑1 (MQ=255)  tGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATcagcg                                                                           <  2:553313/81‑1 (MQ=255)    aaGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATcagcgg                                                                          <  2:553314/80‑1 (MQ=255)     aGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGc                                                                        <  1:553316/81‑1 (MQ=255)         gatcgatAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGAt                                                                    <  1:553318/81‑1 (MQ=255)         gatcgatAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGAt                                                                    >  1:553358/1‑81 (MQ=255)         gatcgatAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGAt                                                                    <  2:553316/81‑1 (MQ=255)           acgatAAAGTGCGGGTGAACAGGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCg                                                                  <  2:553318/80‑1 (MQ=255)              atAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCtt                                                               <  2:553320/81‑1 (MQ=255)              atAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCtt                                                               <  2:553319/81‑1 (MQ=255)                          agaACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCata                                                   <  1:553322/81‑1 (MQ=255)                          agaACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCata                                                   <  1:553323/81‑1 (MQ=255)                          agaACTCGCAAATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCata                                                   >  2:553363/1‑81 (MQ=255)                                 gCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACgt                                            >  2:553358/1‑81 (MQ=255)                                    gatgatTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACgtggt                                         <  2:553323/81‑1 (MQ=255)                                           gCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATc                                   <  2:553324/80‑1 (MQ=255)                                           gCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATc                                   <  2:553325/80‑1 (MQ=255)                                             ccGTTCGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGtt                                >  1:553375/1‑81 (MQ=255)                                                  cGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGcc                            <  1:553328/80‑1 (MQ=255)                                                  cGGCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGcc                            <  1:553327/80‑1 (MQ=255)                                                   ggCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                          >  1:553365/1‑81 (MQ=255)                                                   ggCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                          >  1:553361/1‑81 (MQ=255)                                                   ggCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                          >  1:553360/1‑81 (MQ=255)                                                   ggCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                          >  2:553362/1‑81 (MQ=255)                                                    gCAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATc                         <  2:553328/81‑1 (MQ=255)                                                     cAACCATCCTCGGGATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATc                         >  1:553368/1‑80 (MQ=255)                                                                  gATCACGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAg           >  1:553378/1‑81 (MQ=255)                                                                       cGGTAATCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTg      >  1:553343/1‑81 (MQ=255)                                                                           aaTCAGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCGCt  >  1:553369/1‑81 (MQ=255)                                                                               |                                                                             GTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCGCT  >  NC_000913/3285955‑3286109

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑