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%	97.6 / NA	37	I204S (ATC→AGC)	wcaA	putative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/22);  total (15/22)

GGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132659‑2132923                                                                                                                                 |                                                                                                                                         gggTGAACAGGTGGTGTTCGTTGCGGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACTGTGATTGATATGCAGGCTCTGCGTCGCctct                                                                                                                            <  2:615070/143‑1 (MQ=255)  ggTAAACCGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                           <  1:16178/143‑1 (MQ=255)  ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGCGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                           <  1:768266/143‑1 (MQ=255)       ataGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                      <  1:377854/141‑1 (MQ=255)       aCAGCTGATATGTTGTGCGGGCGCGGTCGAATTTGGCTGTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                      <  1:736582/143‑1 (MQ=255)          gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGACGCCTCTTCTCCTTTc                                                                                                                   >  2:339153/1‑143 (MQ=255)             gATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGTAGGCTCTGCGTCGCCTCTTCTACTTTCCAc                                                                                                                <  2:490735/143‑1 (MQ=255)                   tttgGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                          <  2:464728/143‑1 (MQ=255)                         tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGGATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTAATTGATAAAAAAGCACTACGTCAACTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                    >  1:452061/1‑143 (MQ=255)                                  cGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATATGCATTGCTCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                           <  1:717523/143‑1 (MQ=255)                                   gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacca                                                                                          >  1:510252/1‑143 (MQ=255)                                             ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCTCGGTTCGCCGTCCTCCCCCTCCGTCCTCAGGaa                                                                                >  1:147297/1‑143 (MQ=255)                                              cgggcgtGCGGTAAACTCGCAAGTCCCCGGAGAACTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg                                                                               <  1:490101/137‑1 (MQ=255)                                                   ttGCGGTAAAAGTGGAAGTTCCCGGATAATTTTTTCGTTGACGTGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTc                                                                          <  1:486146/143‑1 (MQ=255)                                                    tGCGCTAAAAGTGGATGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                         <  2:370646/143‑1 (MQ=255)                                                         tAAAAGTGGAAGTACCCGGATAATTTTTTCGGTGGCGAGGTGCTATGCCTTTCACCGTGGTTGATATGCAGGCCCTGCGTCGCCTCTTCTACTTTCCACGCCCCACCGTACTCCGCCACCATCCGCCTGCAGTCGTCGTAATc                                                                    >  2:434498/1‑143 (MQ=255)                                                         tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCGGGCAGAGGTCGTAAta                                                                    >  2:729401/1‑142 (MQ=255)                                                          caatgTGGACGTTCCCGGTGATTTTTTTCGGTGACGAGGTGATCTTCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACTCACGTTCTTTGGAAGATGTCGTAATcc                                                                   <  2:279372/139‑1 (MQ=255)                                                             ggTGGAAGCTCCCGGAGAATTGTTGCGGTGTCGAGGTGATCTGCATTTCACCGTGAGTGATATGCAGGCTGTGCGTCGCCTGTTCTACTGTCCACGGTGCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                <  2:147433/142‑1 (MQ=255)                                                             aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCAACAGGAAGACATCGTAAtacaat                                                                >  1:661853/1‑138 (MQ=255)                                                              gTGGTAGTTCCTTGAGAATTTTTTCGGTTACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCACCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                               <  2:599616/143‑1 (MQ=255)                                                                 gAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATGCGAACGCACATGTCGCAACCCAGTAccc                                                            >  2:30807/1‑142 (MQ=255)                                                                     tACCCGGAGAATTTTTTCGGTGACCAGGTGATCTGCATTTCACAGTAATCGCCACGCATGCTCTGCGGTCCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                        <  2:564304/143‑1 (MQ=255)                                                                        ccGGAGAATTTTTTCGGTGACGAGGTGATATGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTATTCTACTTTACACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGGTCGTAATCCTCTACAGCttcaag                                                     >  2:148220/1‑139 (MQ=255)                                                                          ggatgaTTTTTTCTGTGACGAGCTGATCTGCATTTCTCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                   <  2:307531/138‑1 (MQ=255)                                                                            cgaATTGTTTCGGTCACCAGGTGATTTGCATTTCACCGTGCTTTATATTCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGTTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                 <  2:270396/142‑1 (MQ=255)                                                                            agtaTTTTTTCGTTGACGATCTGAGCTGTATGTCTCCGTGATTGGTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGTCGTACTCCACCACCATCCGCACGATGCTGTAGGTGTGCGGTGCCGCTTTCAGTTCg                                                 <  2:279424/140‑1 (MQ=255)                                                                              aaTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATACGCAGGAAGATGGCGTAATCCTGTGCCGCCGTCAGTTcact                                               >  2:417766/1‑140 (MQ=255)                                                                                    ttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGCGCCTCTTCTACTTTCCACGGTTCACCGTACTCCACCACAATCCGCAGGAAGATATCGTAATCCTGTACCACTGGTAGGTCGGTATCGaa                                         >  2:152367/1‑143 (MQ=255)                                                                                         tGACGAGGTGATCTGCATTTCACCATGATTGATATGCAGGCTCGGCGTCGCCTCTTCTACTTTCCACGGGTCGCCGCACTCCCCCACCATCAGCAGGACGATGTCGAAAGCCTGTGCCGCGCACAGTCAGGTAGCgaagaacc                                    >  2:403532/1‑141 (MQ=255)                                                                                               ggTGATCTGCATGTGACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACCCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                              <  2:426865/143‑1 (MQ=255)                                                                                                   aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGAACAAACACTCTTTGaa                          >  2:443858/1‑143 (MQ=255)                                                                                                     cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCATCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGCTGTCGTAATCCTGTTCCGCTTTAAGTTCGGGATCGCCCATACCCTCTTTGAAAc                        >  1:234378/1‑143 (MQ=255)                                                                                                               ccGTGATTGCGATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCa              <  2:436824/143‑1 (MQ=255)                                                                                                                 gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTACGGAATAAAAAAAAAAATAAATGAAAAGCAATGCCAAgg            >  1:420610/1‑143 (MQ=255)                                                                                                                        atatCCAGGCTCTGCGTCTCCTGTTCTACTTTTCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGACGTCGAAATCCTGTGCCGCTTTCAGTTCGGTATCGAAAAAACACTCTTTGAAATGCCATGCCCAGGTAAAGAc     <  2:346695/143‑1 (MQ=255)                                                                                                                           ttcAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg  <  1:733141/141‑1 (MQ=255)                                                                                                                                 |                                                                                                                                         GGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132659‑2132923

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