Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,657,262	C→T	E273E (GAG→GAA)	wzyE ←	predicted Wzy protein involved in ECA polysaccharide chain elongation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,657,262	0	C	T	100.0%	222.7 / NA	68	E273E (GAG→GAA)	wzyE	predicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/40);  total (28/40)

GGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCAT  >  W3110S.gb/3657206‑3657317                                                         |                                                        ggAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:3450039/62‑1 (MQ=255) ggAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:3365309/62‑1 (MQ=255) ggAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:3897219/62‑1 (MQ=255) ggAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:1452357/62‑1 (MQ=255) ggAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:3831097/62‑1 (MQ=255)           ccTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAg                                          <  1:2841264/62‑1 (MQ=255)                 aGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcg                                    <  1:3529301/62‑1 (MQ=255)                 aGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcg                                    <  1:4673551/62‑1 (MQ=255)                 aGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcg                                    <  1:4126270/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:2667166/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:4076459/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:2683760/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:4078619/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:215198/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:2144019/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:3206875/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:4468886/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:3329206/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:1424155/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:527719/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:705577/62‑1 (MQ=255)                  gTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:1059469/62‑1 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTcc                                                      >  1:3011469/1‑40 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:30497/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:2856212/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:2671468/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:2275897/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:3968705/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:2545909/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:2212333/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:1879140/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:1866253/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:1714106/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:1451553/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGa                                 >  1:1444572/1‑61 (MQ=255)                     gATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTTTCCAACGGCGAGAAGGTGTcgc                                   >  1:2208420/1‑59 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:103577/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:3885131/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:4284600/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:3764026/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:903742/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:963899/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:3442539/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:3335897/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:1616050/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:2434524/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCa                                                    <  1:2572802/36‑1 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga                             >  1:1025130/1‑59 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga                             >  1:110754/1‑59 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga                             >  1:1259731/1‑59 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga                             >  1:1952420/1‑59 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga                             >  1:2908780/1‑59 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga                             >  1:3842076/1‑59 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga                             >  1:3620716/1‑59 (MQ=255)                           gtCGTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTga                             >  1:378799/1‑59 (MQ=255)                              gTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATag                        <  1:361088/61‑1 (MQ=255)                              gTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATag                        <  1:3159935/61‑1 (MQ=255)                              gTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATag                        <  1:669578/61‑1 (MQ=255)                              gTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATag                        <  1:1667604/61‑1 (MQ=255)                              gTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATag                        <  1:560662/61‑1 (MQ=255)                              gTAGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATag                        <  1:1372753/61‑1 (MQ=255)                                aGTTCTGCAACAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTcgc                                   <  1:3904364/48‑1 (MQ=255)                                        aaCAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGAt                          >  1:1533107/1‑49 (MQ=255)                                        aaCAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGAt                          >  1:4643124/1‑49 (MQ=255)                                        aaCAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGAt                          >  1:4143566/1‑49 (MQ=255)                                        aaCAGCAACGCCAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGAt                          >  1:3279476/1‑49 (MQ=255)                                                   cAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCAt  >  1:2589245/1‑62 (MQ=255)                                                   cAGATTTTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCAt  >  1:3568069/1‑62 (MQ=255)                                                         |                                                        GGAGCCAGGCCCTGGAAGTCGATGTTGTCGTAGTTCTGCAACAGCAACGCCAGATTCTCCCACGGCGAGAAGGTGTCGCGAGTGAGATAGAGAAACGTATAGAACGCTTCAT  >  W3110S.gb/3657206‑3657317

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