Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,132,787	A→C	100%	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%	76.4 / NA	23	I204S (ATC→AGC)	wcaA	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/11);  total (12/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGT  >  NC_000913/2132517‑2133039                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             gCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTAACCGTGATTGATATGCAGGctct                                                                                                                                                                                                                                                           >  1:352467/1‑274 (MQ=255)                   ggCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                                                                                                                                        >  2:236243/1‑275 (MQ=255)                                         tttAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGACATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTTATATGCAGGCTCTGCGTCGCTTCTTCTACTTTCCACGGTTCGCTGTATTccacc                                                                                                                                                                                                                  >  1:351171/1‑275 (MQ=255)                                           tAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                                                                                                                                                >  1:156532/1‑275 (MQ=255)                                                      aattacaccAGCCAGACGTTTGCCCGTAAGCACCGACAGCGGGGTCAACAAAGTGGGCCAGGTCATACGTTTATTTCGGATCTGGTAGGGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                                                                                                                                                                                                     <  2:351171/267‑1 (MQ=255)                                                                                    cACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                                                                                                                                       <  2:156532/275‑1 (MQ=255)                                                                                                cGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTc                                                                                                                                                                                                                               >  1:206994/1‑207 (MQ=255)                                                                                                 gTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaa                                                                                                                                                          >  2:105943/1‑275 (MQ=255)                                                                                                   cAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaac                                                                                                                                                        <  1:310229/275‑1 (MQ=255)                                                                                                    aGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                                                                                                                                           <  2:206994/207‑1 (MQ=255)                                                                                                                             tATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGACTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGTCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATGCGCAGGAAGATGTCGTACTCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGt                                                                                                                               >  2:169853/1‑274 (MQ=255)                                                                                                                                    gATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAAGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt                                                                                                                       <  1:236243/275‑1 (MQ=255)                                                                                                                                                    aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGatat                                                                                                        <  1:243225/274‑1 (MQ=255)                                                                                                                                                         cTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGt                                                                                                   >  1:174255/1‑274 (MQ=255)                                                                                                                                                            atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCAt                                                                                                                                       >  2:329272/1‑235 (MQ=255)                                                                                                                                                                        ggggcgGTCGAATTGTTCTTTGGGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAggcggc                                                                                   <  2:352467/272‑1 (MQ=255)                                                                                                                                                                                 gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACg                                                                          <  2:174255/275‑1 (MQ=255)                                                                                                                                                                                                                 aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGAATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCGGGTTGCAAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGt                                          <  1:105943/275‑1 (MQ=255)                                                                                                                                                                                                                                                      cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCCTGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGCAGAACAGGCGGCGTGAAAACGGCGATGTGGGATACAGCGGCAAGATTGCCGGGTAGGAAAAGACTTGCCCCTGGAAAACCTTGACGTTaa     >  2:27184/1‑274 (MQ=255)                                                                                                                                                                                                                                                         ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTGGGTATCGATCAAACACCCTTTGCAACGCCATGACCAGGTAAAGACCTGGTTGCCAACGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTCGCAAACGTCGTCGTTAGCGt  >  1:223438/1‑275 (MQ=255)                                                                                                                                                                                                                                                                 gattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACg               >  2:384298/1‑254 (MQ=255)                                                                                                                                                                                                                                                                     gatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGt           <  1:384298/254‑1 (MQ=255)                                                                                                                                                                                                                                                                        atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGc    >  2:188400/1‑258 (MQ=255)                                                                                                                                                                                                                                                                        atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGc    <  1:188400/258‑1 (MQ=255)                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             GCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGT  >  NC_000913/2132517‑2133039

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑