Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,132,787	A→C	I204S (ATC→AGC)	wcaA ←	putative 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%	124.9 / NA	45	I204S (ATC→AGC)	wcaA	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/30);  total (15/30)

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGA  >  NC_000913/2132647‑2132919                                                                                                                                             |                                                                                                                                     ggAGCTGGTAGAGGGTAGACAGCTGATATTTTTTGCTGGCGCGGTCGTATTTGTCTTTGTGCTTGCGGTACAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                    <  2:788518/143‑1 (MQ=255)       ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATATGCAGGCTCTGCGTc                                                                                                                              >  1:673226/1‑143 (MQ=255)       ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGGCGAATTTGTCTCTGTGCTTGCGGTAACAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGACCCGCATTTCCCCCTGACCGATACGCAGGCGCTACGTc                                                                                                                              >  2:259602/1‑143 (MQ=255)         tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGCCTCTGCGTCGc                                                                                                                            >  2:229033/1‑143 (MQ=255)             gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCGGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCcgct                                                                                                                        >  1:453567/1‑143 (MQ=255)                 aagtaGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                    <  2:470000/139‑1 (MQ=255)                  aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAc                                                                                                                   >  1:675916/1‑143 (MQ=255)                      gCTGATATTTTTTGCTGGCGCGGTCGTATTTGTTTTTTTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                               <  1:334804/143‑1 (MQ=255)                      gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                               <  1:337412/143‑1 (MQ=255)                       cTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTcc                                                                                                              >  1:609502/1‑143 (MQ=255)                        tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                             <  1:99953/143‑1 (MQ=255)                          atatTTTTGGCTGGAGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                           <  1:580678/143‑1 (MQ=255)                           tatTGTTTGCTTGCGCGGTCGAAGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTGCCCGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACgg                                                                                                          <  2:564011/143‑1 (MQ=255)                                    cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                 >  2:52008/1‑143 (MQ=255)                                      ggCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACt                                                                                               >  1:3847/1‑143 (MQ=255)                                       gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                              <  1:491232/143‑1 (MQ=255)                                       gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                              >  1:755053/1‑143 (MQ=255)                                       gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                              <  1:275176/143‑1 (MQ=255)                                            gTCGATTTTGTCGTTGGGCTTGCTGTAAAAGTTGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                         <  2:107160/143‑1 (MQ=255)                                                 aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGGACTCCACCACCATc                                                                                    >  2:549348/1‑143 (MQ=255)                                                        tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                                                                             >  1:359994/1‑143 (MQ=255)                                                         ttggacTTGCGGTAAAAGACGACGTACCCCGAGCATGTTTTCGGCGACGAGGTGATCTGCATTTCACCGTGATTGCTACGCAGGCTCTGCGTCGCCTCTTCTACATTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                            <  1:766200/138‑1 (MQ=255)                                                          tgtgCTTGTGGTACAAGTGGAAGTACCCGGAGAATTTTGTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg                                                                           <  2:517088/143‑1 (MQ=255)                                                            tgCTTGCGGTAAAAGTGTAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                         <  2:209444/143‑1 (MQ=255)                                                              ccggCGGTCGAAGTGGAAGTACTCGGATAATGTGTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt                                                                       <  2:576984/140‑1 (MQ=255)                                                                tGCGCTAAAAGTCGACGTACCCTGAGAGTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                     <  1:340532/143‑1 (MQ=255)                                                                 gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGt                                                                    >  1:155864/1‑143 (MQ=255)                                                                     tAAAGGTGGAAGTGCCCGGAGAATTTTTTCGGGGGCGAGGTGATCTGCATTTCACCGTGATTGATAGGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                <  2:92525/143‑1 (MQ=255)                                                                           tgccaGTTCCCGGAGCAGTTTTTCGGTGGCGAGGTGTTCTGTATTTCCCCGTGATTGATATGCAGGCTGTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc                                                          <  2:234736/139‑1 (MQ=255)                                                                                 tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                    <  1:352329/143‑1 (MQ=255)                                                                                  aCCCGGAGACTGTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTc                                                   <  1:747122/143‑1 (MQ=255)                                                                                  aCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGGAGATGTCGTAATCCTGTGCCGCtcta                                                   >  1:685471/1‑140 (MQ=255)                                                                                   cccGGAGAAGTTTTTCGGTGACGCGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                  <  2:673361/143‑1 (MQ=255)                                                                                     cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAAGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACCCCCCCACCACCCGCAGGAACATAGCGTAATACAATGCCGCTCGCAGt                                                >  2:403449/1‑143 (MQ=255)                                                                                      ggAGAATTTTTTCGGTGACGAGGTGATCTGCATCTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                               <  1:156268/143‑1 (MQ=255)                                                                                         gaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                            <  1:562672/143‑1 (MQ=255)                                                                                           aTTTTTTCGGTCACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTa                                          <  1:811798/143‑1 (MQ=255)                                                                                                  cgttGACGGGGTGGTCTGCCTTTCTGCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGTAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaca                                   <  2:675916/140‑1 (MQ=255)                                                                                                      gtcGAGGTGCTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTCCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                               <  2:546621/141‑1 (MQ=255)                                                                                                            gtgtcttGCATTTCGCCGTGATTGTTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCTGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                         <  2:608988/137‑1 (MQ=255)                                                                                                               aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCGTTGaa                      >  1:812472/1‑143 (MQ=255)                                                                                                                 ctctaTTTCACCGTTGTTGTTATGCAGTCTCTGCGTCGCCTCTTGTACTTTCCACGTTTGTCCGTACTCCACCACCATTCGCAGGAATATTTCGTAATCCTGTGCCGCTTTCATTTCGGTATCGAACAAACACTCTTTGAAAc                    <  2:307217/139‑1 (MQ=255)                                                                                                                                 ttgatATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCCCTGTTCGCCGCAGTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa    <  2:115179/143‑1 (MQ=255)                                                                                                                                  tgctatGCAGGCTCTGCGTCGCCTCTGCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAg   <  2:748748/140‑1 (MQ=255)                                                                                                                                   gatACGCAGGCTCTGCGTCCCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGa  <  1:167877/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                     GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGA  >  NC_000913/2132647‑2132919

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