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%	148.7 / NA	48	R81C (CGT→TGT)	wcaK	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/25);  total (23/25)

TGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCC  >  NC_000913/2118044‑2118297                                                                                                                         |                                                                                                                                      tGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAAcgcg                                                                                                                     >  2:1009512/1‑139 (MQ=255) tGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAAcgcg                                                                                                                     >  2:208689/1‑139 (MQ=255)                aCAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagag                                                                                                                                    >  2:119663/1‑109 (MQ=255)                aCAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagag                                                                                                                                    <  1:119663/109‑1 (MQ=255)                   gcagGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCcgccgc                                                                                                   >  1:695187/1‑139 (MQ=255)                    cagGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCcgcg                                                                                                  <  1:6322/139‑1 (MQ=255)                     agGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCgcgc                                                                                                 <  1:317262/139‑1 (MQ=255)                                         tCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGc                                                                             <  1:484538/139‑1 (MQ=255)                                                  gatggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGc                                                                    <  2:13974/139‑1 (MQ=255)                                                   atggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCg                                                                   <  1:1009512/139‑1 (MQ=255)                                                    tggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCgg                                                                  >  1:1492898/1‑139 (MQ=255)                                                    tggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCgg                                                                  >  2:58077/1‑139 (MQ=255)                                                    tggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCgg                                                                  >  2:512268/1‑139 (MQ=255)                                                     ggcgatATTGTGCAGCTTGACAGAGTCGGGGACGCGTGAGAGCAATACCTGATCCTGGTAGCGGCGACAGAGTACTTTTTTACCGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                                 <  1:851364/139‑1 (MQ=255)                                                     ggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                                 <  1:632779/139‑1 (MQ=255)                                                     ggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                                 <  2:1492898/139‑1 (MQ=255)                                                       cgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATc                                                               <  2:695187/139‑1 (MQ=255)                                                         atatTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGc                                                             <  1:208689/139‑1 (MQ=255)                                                                 cAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTAcc                                                     >  1:1279845/1‑139 (MQ=255)                                                                    cTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAg                                                                           >  2:1349284/1‑114 (MQ=255)                                                                    cTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAg                                                                           <  1:1349284/114‑1 (MQ=255)                                                                    cTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGTATCGCCCATTACCggg                                                  >  2:1093921/1‑139 (MQ=255)                                                                    cTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCggg                                                  >  1:460917/1‑139 (MQ=255)                                                                      tGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGTGCg                                                <  1:874371/139‑1 (MQ=255)                                                                      tGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCg                                                >  1:254812/1‑139 (MQ=255)                                                                              tCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                                    <  1:241045/79‑1 (MQ=255)                                                                              tCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                                    >  2:241045/1‑79 (MQ=255)                                                                                 gTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAAcgc                                                                                                                      >  1:733609/1‑58 (MQ=255)                                                                                 gTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAAcgc                                                                                                                      <  2:733609/58‑1 (MQ=255)                                                                                 gTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGagcagc                                      <  1:541827/138‑1 (MQ=255)                                                                                 gTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGagcagc                                      >  2:541827/1‑138 (MQ=255)                                                                                      gcgTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                          >  1:245584/1‑129 (MQ=255)                                                                                      gcgTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                          <  2:245584/129‑1 (MQ=255)                                                                                           agagCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGagc                                         <  1:206326/125‑1 (MQ=255)                                                                                           agagCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGagc                                         >  2:206326/1‑125 (MQ=255)                                                                                                       gATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                          >  2:1087918/1‑112 (MQ=255)                                                                                                       gATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                          <  1:1087918/112‑1 (MQ=255)                                                                                                         tGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCTAGGAAGAACTGACCGGATAGCGGc             <  1:1093921/139‑1 (MQ=255)                                                                                                          gCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgttgt                                                                                           >  2:1101719/1‑60 (MQ=255)                                                                                                          gCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgttgt                                                                                           <  1:1101719/60‑1 (MQ=255)                                                                                                            tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtca          <  2:460917/139‑1 (MQ=255)                                                                                                              gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAg                                                                           <  2:74/72‑1 (MQ=255)                                                                                                              gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAg                                                                           >  1:74/1‑72 (MQ=255)                                                                                                              gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCgg              >  1:1011491/1‑133 (MQ=255)                                                                                                              gTAGCGGCGACAGAGGACTTTTTTAACGCGACCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCgg              <  2:1011491/133‑1 (MQ=255)                                                                                                                  cggcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGc                                                                             >  2:1514052/1‑66 (MQ=255)                                                                                                                  cggcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGc                                                                             <  1:1514052/66‑1 (MQ=255)                                                                                                                    gcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                          >  2:538243/1‑99 (MQ=255)                                                                                                                    gcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                          <  1:538243/99‑1 (MQ=255)                                                                                                                    gcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATcc  >  1:104887/1‑139 (MQ=255)                                                                                                                         |                                                                                                                                      TGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACATCC  >  NC_000913/2118044‑2118297

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