Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,117,185	T→C	E407G (GAG→GGG)	wcaK ←	putative colanic acid biosynthesis pyruvyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,117,185	0	T	C	100.0%	146.3 / NA	41	E407G (GAG→GGG)	wcaK	putative colanic acid biosynthesis pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/20);  total (21/20)

TAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTGCCATCTCCGGTAGTCCCAGCTGCTGCAT  >  NC_000913/2117052‑2117321                                                                                                                                      |                                                                                                                                         tAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGAc                                                                                                                                    >  1:886853/1‑140 (MQ=255) tAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGAc                                                                                                                                    >  1:222127/1‑140 (MQ=255) tAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGAc                                                                                                                                    >  1:711111/1‑140 (MQ=255)   aaCGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTg                                                                                                                                  <  2:886853/140‑1 (MQ=255)     cgcgGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAGAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGAc                                                                                                                                >  1:772716/1‑140 (MQ=255)       cgGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACgg                                                                                                                              <  2:11174/140‑1 (MQ=255)        ggTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGc                                                                                                                             >  2:135626/1‑140 (MQ=255)                           aaGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTc                                                                                                          >  1:545742/1‑140 (MQ=255)                            aGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCa                                                                                                         >  1:175585/1‑140 (MQ=255)                                 tCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGcc                                                                                                    >  1:488106/1‑140 (MQ=255)                                      cgacAGCGGGAATTTCAGTAAAGAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAg                                                                                               <  2:772716/140‑1 (MQ=255)                                      cgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAg                                                                                               <  1:135626/140‑1 (MQ=255)                                            cGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGcc                                                                                         >  1:137390/1‑140 (MQ=255)                                             gggAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCt                                                                                        <  2:711111/140‑1 (MQ=255)                                                  tttCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt                                                                                   <  2:175585/140‑1 (MQ=255)                                                           aaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAAc                                                                          >  2:224314/1‑140 (MQ=255)                                                              gaagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCAt                                                                       <  2:488106/140‑1 (MQ=255)                                                                     aCCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGc                                                                >  1:740094/1‑140 (MQ=255)                                                                      ccTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                                                               <  2:713796/140‑1 (MQ=255)                                                                      ccTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCa                                                               <  2:740094/140‑1 (MQ=255)                                                                             ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGcc                                                        <  2:137390/140‑1 (MQ=255)                                                                                    ccccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCtaat                                                 >  2:589053/1‑140 (MQ=255)                                                                                     cccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCt                                                           >  1:806387/1‑129 (MQ=255)                                                                                     cccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCt                                                           >  1:357044/1‑129 (MQ=255)                                                                                     cccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCt                                                           <  2:806387/129‑1 (MQ=255)                                                                                     cccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCt                                                           <  2:357044/129‑1 (MQ=255)                                                                                     cccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCg                                                       >  2:608178/1‑133 (MQ=255)                                                                                     cccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCg                                                       <  1:608178/133‑1 (MQ=255)                                                                                       cGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATaaa                                              <  1:917435/140‑1 (MQ=255)                                                                                         aTGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATg                                            <  2:115823/140‑1 (MQ=255)                                                                                          tGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGa                                           <  1:197345/140‑1 (MQ=255)                                                                                              cTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGa                                       >  1:362570/1‑140 (MQ=255)                                                                                               tCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGat                                      >  1:848320/1‑140 (MQ=255)                                                                                                   gCACAGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCa                                  <  2:545742/140‑1 (MQ=255)                                                                                                       aGACTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTg                              <  1:589053/140‑1 (MQ=255)                                                                                                         aCTGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTGcc                            >  1:355396/1‑140 (MQ=255)                                                                                                           tGCACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTGCCAt                          <  2:222127/140‑1 (MQ=255)                                                                                                             cACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTGCCAtct                        >  1:891143/1‑140 (MQ=255)                                                                                                             cACCATCTGCATTCCTGTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTGCCAtct                        <  2:891143/140‑1 (MQ=255)                                                                                                                             gTCTGACGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTGCCATCTCCGGTAGTCCCAgctg        >  2:607273/1‑140 (MQ=255)                                                                                                                                   cGCCCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTGCCATCTCCGGTAGTCCCAGCTGCTGCAt  >  1:659525/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                         TAAACGCGGTAATTTGGTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCAACCATCGCTTGCAGGCTGCCGTCTAATAAATGACGGATATCAATTGCCATCTCCGGTAGTCCCAGCTGCTGCAT  >  NC_000913/2117052‑2117321

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