Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,122,677	A→G	F152S (TTC→TCC)	wcaE ←	colanic acid biosynthesis glycosyltransferase WcaE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,122,677	0	A	G	100.0%	37.3 / NA	12	F152S (TTC→TCC)	wcaE	colanic acid biosynthesis glycosyltransferase WcaE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (4/8);  total (4/8)

GCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTA  >  NZ_CP009273/2122590‑2122739                                                                                        |                                                               gCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGgaa                                                              <  2:335654/90‑1 (MQ=255)   tttATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTAAAGCCGGATGCAGGGgaaaa                                                            >  1:256163/1‑90 (MQ=255)           ttttGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCt                                                    >  2:263107/1‑90 (MQ=255)            tttGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCtg                                                   <  2:92724/90‑1 (MQ=255)                             aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGc                                  <  1:263107/90‑1 (MQ=255)                             aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGc                                  <  1:275417/90‑1 (MQ=255)                             aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGc                                  <  2:394664/90‑1 (MQ=255)                               tCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCtg                                <  1:391634/90‑1 (MQ=255)                                            tataTTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATAcc                   >  2:59366/1‑90 (MQ=255)                                                 tCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGccc              <  2:14134/90‑1 (MQ=255)                                                     ggTCATAACGCCATTTTTTCAAGCCGGATACGGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGGTAAATATACCAGCCCgggt          >  1:308140/1‑88 (MQ=255)                                                          tAACGCCATTTTTTCAAACCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTTgcg     >  1:344285/1‑90 (MQ=255)                                                             cGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTa  <  1:243105/90‑1 (MQ=255)                                                             cGCCATTTTTTCAAGCCGGATACAGGGGAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTa  <  1:59366/90‑1 (MQ=255)                                                                                        |                                                               GCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTA  >  NZ_CP009273/2122590‑2122739

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