Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,124,626	G→T	A374E (GCG→GAG)	wcaI ←	putative colanic biosynthesis glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,124,626	0	G	T	100.0%	138.2 / NA	38	A374E (GCG→GAG)	wcaI	putative colanic biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (18/20);  total (18/20)

ACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCA  >  NC_000913/2124550‑2124702                                                                             |                                                                             aCGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCt                                                                           >  1:386016/1‑80 (MQ=255)   gTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGAc                                                                        <  1:385969/81‑1 (MQ=255)      acacGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCATGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACgg                                                                      <  2:385970/80‑1 (MQ=255)       cacGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGAt                                                                    <  1:385971/81‑1 (MQ=255)          gTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATccc                                                                 <  1:385972/81‑1 (MQ=255)             ctctTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGc                                                              >  1:386018/1‑81 (MQ=255)              tctTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGc                                                              >  1:386019/1‑80 (MQ=255)               ctTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGcca                                                            >  1:386015/1‑81 (MQ=255)                tttATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGccac                                                           >  1:386025/1‑81 (MQ=255)                tttATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGccac                                                           >  2:386030/1‑81 (MQ=255)                 ttATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGccac                                                           >  2:386031/1‑80 (MQ=255)                 ttATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGccac                                                           >  2:386028/1‑80 (MQ=255)                    tCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAgc                                                       <  2:385973/81‑1 (MQ=255)                       aGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGcc                                                    >  1:386029/1‑81 (MQ=255)                        gcgtgcgtTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCt                                                   <  1:385974/81‑1 (MQ=255)                        gcgtgcgtTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCt                                                   >  1:386034/1‑81 (MQ=255)                         cgtgcgtTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTc                                                  <  1:385975/81‑1 (MQ=255)                          gtgcgtTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTc                                                  >  1:386020/1‑80 (MQ=255)                          gtgcgtTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCg                                                 <  2:385976/81‑1 (MQ=255)                            gcgtTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGAc                                               >  1:386022/1‑81 (MQ=255)                             cgtTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGAcc                                              >  1:386024/1‑81 (MQ=255)                                   gCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTcc                                        >  2:386026/1‑81 (MQ=255)                                      tatTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGt                                     >  1:386013/1‑81 (MQ=255)                                           aCGTGCCACCGTGTTGTGCTTGGGCAGCAGGGGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAAc                                <  1:385977/81‑1 (MQ=255)                                               gCCACCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCa                             >  1:386033/1‑80 (MQ=255)                                                  aCCGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAAcc                         <  1:385978/81‑1 (MQ=255)                                                    cGTGTTGTGTTTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGc                       <  2:385979/81‑1 (MQ=255)                                                        ttgtgttTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATg                   <  2:385981/81‑1 (MQ=255)                                                        ttgtgttTGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATg                   <  1:385980/81‑1 (MQ=255)                                                             tttGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCgg              <  2:385982/81‑1 (MQ=255)                                                             tttGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCgg              <  1:385983/81‑1 (MQ=255)                                                               tGGGCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGa             <  1:385984/80‑1 (MQ=255)                                                                gggCAGCAGGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGa             >  1:386040/1‑79 (MQ=255)                                                                      cagGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTc     <  2:385985/81‑1 (MQ=255)                                                                      cagGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTc     <  1:385986/81‑1 (MQ=255)                                                                       agGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCg    <  2:385987/81‑1 (MQ=255)                                                                       agGAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCg    >  2:385988/1‑81 (MQ=255)                                                                         gAGCTCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCa  <  2:622209/81‑1 (MQ=255)                                                                             |                                                                             ACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCA  >  NC_000913/2124550‑2124702

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