Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,132,787	A→C	I204S (ATC→AGC)	wcaA ←	putative colanic acid biosynthesis glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,132,787	0	A	C	100.0%	68.1 / NA	26	I204S (ATC→AGC)	wcaA	putative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (10/16);  total (10/16)

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132647‑2132923                                                                                                                                             |                                                                                                                                         ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCTTTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                        <  1:274698/143‑1 (MQ=255)  gATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTCAAAGTGGAAGTACCCGGCGAATTTTTTCGGTGACGAAGTGATCTACATTTCACCCTGATTGAGATGCCGGctct                                                                                                                                       >  2:440714/1‑143 (MQ=255)      tGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAGGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCGCTGCGt                                                                                                                                   >  1:515273/1‑143 (MQ=255)        gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGGATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCAGTTCCCAGTGATTGATATGCAGGCGCAGCGTCg                                                                                                                                 >  2:54573/1‑143 (MQ=255)        atatagtgTATACAGCTGATATTTGTTGTTGGTGTGTTCGTATTTGTCTTTTCGCTTTCGGTAAAAGTTGATGTACTCGGAGAATTTTTTCGGTGACGAGGTTATCTGCTTTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                 <  2:450155/136‑1 (MQ=255)              gttAAACAGCTTATCTTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                           <  2:93103/141‑1 (MQ=255)              ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                           >  1:223638/1‑143 (MQ=255)                        tgttagttttgGCGGGCGCGTTCGCATTTGTCTTTGTGCTTGCGGTAACAGTGGAAGTACTCGGAGAATTTTTGCGGTGACTAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                 <  2:148954/137‑1 (MQ=255)                        tGATATCTTTTGCTGGCGCGGTCTAATTTGTCTTTGTGTTTGCGGTAAAAGTGGAAGTACCCGCAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                 <  2:638399/143‑1 (MQ=255)                           tatTTTTTGCTGGCGCGGTCGAATTTGTGGTATTTCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACgg                                                                                                              <  1:446662/143‑1 (MQ=255)                                     tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                    >  1:500349/1‑143 (MQ=255)                                             tCGAGTTTGTTTTTGTGCTTGCCGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                            <  1:368597/143‑1 (MQ=255)                                             tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTCTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGGTCGCAGTCGTccacgac                                                                                            >  2:533960/1‑143 (MQ=255)                                               gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccaccc                                                                                          >  1:445410/1‑142 (MQ=255)                                                     ggcgttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                    <  1:659893/139‑1 (MQ=255)                                                           gtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTACACGGTTCGTCGTAATCCACCACCATCCGTAGTAAGa                                                                              >  2:210583/1‑143 (MQ=255)                                                                    gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCATCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt                                                                     <  1:310840/143‑1 (MQ=255)                                                                                 tccccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCGCGGTTCGCCGGACTCCCCCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                        <  2:614800/141‑1 (MQ=255)                                                                                    ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTATGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGCTGTCGTAGTCCTGTGCAGCATTCAg                                                     >  2:506105/1‑143 (MQ=255)                                                                                     cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTACACGGTTCACCGTACTCCACCACCATCCACAGTAAGATGTCGTACTCCTGTGCCGCTTTCAGt                                                    >  2:514954/1‑143 (MQ=255)                                                                                         gaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                <  1:85144/143‑1 (MQ=255)                                                                                           atgtcttCCGTGACGAGGTGATCTGCATTCCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTa                                              <  1:238364/138‑1 (MQ=255)                                                                                             tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTGCAGTTCGGTATc                                            >  1:554932/1‑143 (MQ=255)                                                                                                                                attgctATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTACGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCGGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaa         <  1:440714/143‑1 (MQ=255)                                                                                                                                attgctATCCAGGCTCTGCGTCGCCTCTTCTACTTCCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaa         <  1:364937/143‑1 (MQ=255)                                                                                                                                       tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg  <  1:618129/143‑1 (MQ=255)                                                                                                                                       tGCAGGCTCTGCGTCGCCTCTCCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg  <  1:225835/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                         GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132647‑2132923

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