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%	80.5 / NA	26	I204S (ATC→AGC)	wcaA	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (14/12);  total (14/12)

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAA  >  NC_000913/2132650‑2132917                                                                                                                                          |                                                                                                                                   gcTGGTATAGGTTAAGCAGCTGATATTTTTTGCCGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTCGACGCACCCGGAGAATTTTTTCGGTGACGATGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                               <  2:402665/142‑1 (MQ=255)    gggagggggTACACAGCTGATATTTTTTGCTGGCGCGGTCGAATGTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                            <  1:133058/137‑1 (MQ=255)           ggTATACAGCCGATATTTTTTGGCGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                     <  2:166488/143‑1 (MQ=255)              gaaCAGCTGATATTTGTTGCTGGCGCGGGCGAATTTGGCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                  <  2:177145/142‑1 (MQ=255)                aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                >  1:446871/1‑143 (MQ=255)                                    gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGGCGAGGTGATCTGCATTTCACGGTGATTGATATGCAGGCTCTGCGTCGCCTCTCGTACTTTCCGCGCTTCGCGGTACTc                                                                                            >  1:122588/1‑143 (MQ=255)                                    gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                            <  1:472726/143‑1 (MQ=255)                                    gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGACGCCTCTTCTACTTTCGAGGGTTCGCCGTACTc                                                                                            >  1:146896/1‑143 (MQ=255)                                     cgcgCTCGAATTTGTCTTTCTCCTTCCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTGCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                                                                           <  1:117580/143‑1 (MQ=255)                                          tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACCccaccac                                                                                      >  2:384173/1‑143 (MQ=255)                                               tttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACCCCCCCACCCTca                                                                                 >  2:374490/1‑142 (MQ=255)                                                         tgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGCAGAt                                                                       >  1:250541/1‑143 (MQ=255)                                                                ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                >  1:318997/1‑143 (MQ=255)                                                                 gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTGCACGGTTCGCCGTCCCACACCACCCTCCACAGGCCGATGTCCTAAt                                                               >  2:151714/1‑143 (MQ=255)                                                                      tgTGGACGTCCCCGGAGAATTGTTGCGGTGACGAGGTGCTCTGCCTTTCTCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                          <  2:311227/142‑1 (MQ=255)                                                                       gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCAGCCGCGGGAGGATGTCGTAAACCtgtg                                                         >  1:419510/1‑143 (MQ=255)                                                                             gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTAATCCTGTCCAGCtt                                                   >  2:505188/1‑143 (MQ=255)                                                                             gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtg                                                   >  1:487114/1‑142 (MQ=255)                                                                                    gagaGTTTTTTCGGTGACGGGGGGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                            <  1:72538/143‑1 (MQ=255)                                                                                                 gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                               <  1:381080/143‑1 (MQ=255)                                                                                                                 cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGACCAACCACTCTTTGAAACGcc               >  1:464292/1‑143 (MQ=255)                                                                                                                  aTTTCACCGCGATTGATATGCCGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTGATCCTGTGCCGCTTTCAGATCGGTATCGAATAAACACTCTTTGAAACGCCa              <  2:272916/143‑1 (MQ=255)                                                                                                                   tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCAt             >  1:394682/1‑143 (MQ=255)                                                                                                                       aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGCGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc         <  1:232494/143‑1 (MQ=255)                                                                                                                           tgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAGACGCCATGCCCAGGt     >  1:58466/1‑143 (MQ=255)                                                                                                                              tcgcgatGCGGGCTCTGCGTCGCCTCTTCTACGCTCCACGGTTCGCCGTACCCCACCACCATCCGGAGGAATATGTCGTAAGCCTGTGCCGCTTTCAGTCCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa  <  2:462877/138‑1 (MQ=255)                                                                                                                                          |                                                                                                                                   TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAA  >  NC_000913/2132650‑2132917

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