Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,118,164	G→A	R81C (CGT→TGT)	wcaK ←	colanic acid biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,118,164	0	G	A	100.0%	105.1 / NA	33	R81C (CGT→TGT)	wcaK	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

AATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTC  >  NC_000913/2118042‑2118288                                                                                                                           |                                                                                                                             aaTGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAAcg                                                                                                              <  1:1042055/139‑1 (MQ=255)    ggCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACgcgc                                                                                                           >  1:566359/1‑139 (MQ=255)     gCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCACTACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGcc                                                                                                          >  2:936962/1‑139 (MQ=255)      cgtcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGccc                                                                                                         >  1:274658/1‑139 (MQ=255)        tcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagag                                                                                                                           <  2:322689/119‑1 (MQ=255)        tcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagag                                                                                                                           >  1:322689/1‑119 (MQ=255)        tcgtAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGa                                                                                                       >  1:69028/1‑139 (MQ=255)                                 aGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgttgtgttg                                                                              >  1:540953/1‑139 (MQ=255)                                         aaTCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCAtct                                                                      <  2:566359/139‑3 (MQ=255)                                                 ggcgatggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCcgcg                                                                                         >  1:318950/1‑112 (MQ=255)                                                 ggcgatggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCcgcg                                                                                         <  2:318950/112‑1 (MQ=255)                                                     atggcggtATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                           <  2:911111/106‑1 (MQ=255)                                                     atggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                           >  1:911111/1‑106 (MQ=255)                                                      tggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCcgcg                                                                                         >  1:324991/1‑107 (MQ=255)                                                      tggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCcgcg                                                                                         <  2:324991/107‑1 (MQ=255)                                                                gcgcAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCAtt                                               >  1:310351/1‑139 (MQ=255)                                                                                tCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                        <  1:811761/114‑1 (MQ=255)                                                                                tCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                        >  2:811761/1‑114 (MQ=255)                                                                                tCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGagca                               >  1:1048421/1‑139 (MQ=255)                                                                                   gTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGcc                            <  2:274658/139‑1 (MQ=255)                                                                                    tGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCa                           >  1:233639/1‑139 (MQ=255)                                                                                         cgTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCcgcg                                                                                         <  2:480826/72‑1 (MQ=255)                                                                                         cgTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCcgcg                                                                                         >  1:480826/1‑72 (MQ=255)                                                                                              gagCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTA‑CGGGCGGTTGAGCAGCCAGGAAGAACTGa                <  2:540953/139‑1 (MQ=255)                                                                                                gCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGTAGAACTGAc               <  2:469529/139‑1 (MQ=255)                                                                                                gCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGAc               <  2:310351/139‑1 (MQ=255)                                                                                                           tGCTGGTAGCGGCGACAGAGGACTTTTTTAACgcgc                                                                                                           <  2:985765/36‑1 (MQ=255)                                                                                                           tGCTGGTAGCGGCGACAGAGGACTTTTTTAACgcgc                                                                                                           >  1:985765/1‑36 (MQ=255)                                                                                                             cTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtc  <  2:233639/139‑1 (MQ=255)                                                                                                                gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgtt                                                                                    >  2:1102931/1‑54 (MQ=255)                                                                                                                gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgtt                                                                                    <  1:1102931/54‑1 (MQ=255)                                                                                                                gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                 >  2:861055/1‑105 (MQ=255)                                                                                                                gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                 <  1:861055/105‑1 (MQ=255)                                                                                                                     ggcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCgg     >  2:244293/1‑128 (MQ=255)                                                                                                                     ggcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCgg     <  1:244293/128‑1 (MQ=255)                                                                                                                           |                                                                                                                             AATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTC  >  NC_000913/2118042‑2118288

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