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

TGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132652‑2132925                                                                                                                                         |                                                                                                                                           tggtggAGGGTGAACAGCTGATATTTTGTGCTGGCGCGGTCGAACATACAGCTCAATTTTCTGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                     <  2:512362/143‑1 (MQ=255) gggTAGGGGGGAAAGAGCTGCTATTTTTTGCTGGCGCGGTCGAAGTTGTCGTGGTGCGTGCGGTAAAAGTGGACGTGCACGGAGAATTTTTTCGGTG‑ACGTGGTGCTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                     <  2:529049/142‑1 (MQ=255)  ggTAGAGGGTAAACAGCTAATATATTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGCATTTTTTCGGTG‑ACGAGGTGATCTGCATTTAACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                    >  2:522581/1‑143 (MQ=255)    tAGAGGGTAAGCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                  <  1:712715/143‑1 (MQ=255)    tAGAGGGTAAACAGTGGAGATTTATTGCTAGCGCGGTCTAATTTGTCTTTGCGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                  <  2:559165/143‑1 (MQ=255)            aagcAGCTGATATTGTTGGCTGGCGCGGTCGGAGTTGGCTTTCTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                          <  2:124973/140‑1 (MQ=255)            aaaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTCTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                          <  1:278674/143‑1 (MQ=255)                   tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                   >  2:584176/1‑143 (MQ=255)                   ggATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGGGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                   <  2:167926/142‑1 (MQ=255)                            cggCTGGCGCGGGCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTCACCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                          <  1:769450/141‑1 (MQ=255)                                          gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacca                                                                                            >  2:715074/1‑143 (MQ=255)                                          gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATATCACCGCGAATGATATGCAGGCTCCGCCTCGCCCCTTCTACTTTCCACGCTTAGACGGACTccaccccca                                                                                            >  2:666773/1‑143 (MQ=255)                                                       tgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGCAGCAATCGTTGCAAAG‑TGATGGCCGTGACTACGCCGCTGGCTCTGCGTCGCCCCTTCTGCTGTCCACGGTGCGCCGTACTCCACCACCATCCGCAGGAAggt                                                                               >  1:396938/1‑141 (MQ=255)                                                          ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            <  1:722305/143‑1 (MQ=255)                                                              ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTTGCCGTACTCCACCACCACACGCAGGAAGTAGTCGTaa                                                                        >  2:780962/1‑143 (MQ=255)                                                              ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                        <  1:155519/143‑1 (MQ=255)                                                                    aGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  >  1:179241/1‑143 (MQ=255)                                                                     gTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATATGCATTTCACAGGGATTGATATGCAAGCTCTGCGTCGCCTCTTCTACTTTCCACAGTTCTCCGTACTCCACCACCATCCGCAGCAAGATGTCGTAAACCggtg                                                                 >  2:237054/1‑143 (MQ=255)                                                                        gAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCg                                                              <  1:416532/143‑1 (MQ=255)                                                                        gAAGTACCCGGAGAATTTTTTCGGTG‑ACAAGATGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGACATACTCCACCACCATCAGCAGGAAGATGTCGTAATccactccag                                                              >  2:509784/1‑136 (MQ=255)                                                                             aCCCGGAGAATTTTTTCTGTG‑TCTAGTTGTTCTGCATTTCACCGTGATTGTTATTCAGGCTCTGCGTCGGCGCATCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTc                                                         <  1:715074/143‑1 (MQ=255)                                                                                cGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAGGATGTCGTAAGTCTGTGCCGGTTTCGGt                                                      >  1:634849/1‑143 (MQ=255)                                                                                        tttttCGGTG‑ACGAGGTGCTCTCCATCTCACCGTGTTTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                              <  1:154717/143‑1 (MQ=255)                                                                                           ttCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                           <  1:340683/143‑1 (MQ=255)                                                                                            tCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGaac                                          >  1:89259/1‑143 (MQ=255)                                                                                                       ggTGATCGCCATTCCGCCGTGGTTGAGATGCAGGCTCTGCGTCGCTCGGTCGACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTGCGGTATCGAACAAACACGCtt                                <  2:546058/143‑1 (MQ=255)                                                                                                              tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGTGCCGCCTTCAGTCCGGTACCGGACAAACAGTCTCTGAAACg                         >  2:775791/1‑143 (MQ=255)                                                                                                                cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCCCTTGCAGCTCAGCATAGAACAAACACTCTCTGAACCGcc                       >  2:785450/1‑143 (MQ=255)                                                                                                                 aTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCa                      >  1:774876/1‑143 (MQ=255)                                                                                                                  tttCACGGGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAGTTTCGACGGTTCGCCGTACTCCACCAGCATCCGCAGGAAGATGTGGTAATCCTGTGCGGCGTTCAGTTAGGTAGCGAACTAACCATCTTTGAACCGCCag                     >  2:620807/1‑142 (MQ=255)                                                                                                                      aCCGTGCTTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc                 <  2:231064/143‑1 (MQ=255)                                                                                                                          tcattgatATGCAGGCCCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCGTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGt             <  2:467361/141‑1 (MQ=255)                                                                                                                                     cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGCCCTGGt  >  2:190497/1‑143 (MQ=255)                                                                                                                                         |                                                                                                                                           TGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTG‑ACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132652‑2132925

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