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%	100.2 / NA	29	I204S (ATC→AGC)	wcaA	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (16/13);  total (16/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTG  >  NC_000913/2132507‑2133064                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                      tGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                                                                                                                                                               >  2:229640/1‑289 (MQ=255) tGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                                                                                                                                                               <  1:326576/289‑1 (MQ=255)                    ccAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                                                                                                                                                            >  2:18239/1‑289 (MQ=255)                            tGGCGCGCAGAGCTTCCAGCATCTTTAGCGGCCCCGAATCCCATCAGCCAGCCGTTTGCCGTTACGCACCGGCAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                                                                                                                                                                    <  1:229640/289‑1 (MQ=255)                                                aTCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                                                                                                                                                                >  2:372401/1‑289 (MQ=255)                                                     tAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                                                                                                                                                           <  2:293541/289‑1 (MQ=255)                                                     tAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                                                                                                                                                           <  1:241746/289‑1 (MQ=255)                                                                               aCGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                                                                                                                                                           <  1:172138/263‑1 (MQ=255)                                                                               aCGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                                                                                                                                                           >  2:172138/1‑263 (MQ=255)                                                                               aCGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCACCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCGGt                                                                                                                                                                                                 >  1:72925/1‑289 (MQ=255)                                                                                               aCCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                                                                                                                                                                  <  1:28099/288‑1 (MQ=255)                                                                                                  gccggcagcGTCAGAAAAGGGGGCCAGGTCTTACGTTGATTGCGGGTCTGTTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACAct                                                                                                                                                                              <  1:372401/285‑1 (MQ=255)                                                                                                                                       tATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCCCCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGACCAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACATGGTTGc                                                                                                                                          >  2:54266/1‑288 (MQ=255)                                                                                                                                       tATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGc                                                                                                                                          >  2:2979/1‑288 (MQ=255)                                                                                                                                        aTTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCa                                                                                                                                        <  1:18239/289‑1 (MQ=255)                                                                                                                                                         gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTa                                                                                                                       >  1:365056/1‑289 (MQ=255)                                                                                                                                                           gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg                                                                                                                                              <  2:236078/264‑1 (MQ=255)                                                                                                                                                           gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGAGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg                                                                                                                                              >  1:236078/1‑264 (MQ=255)                                                                                                                                                            tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACTTGGTTGCCAATGATATTCCGTTTGTCGa                                                                                                                     >  1:129157/1‑288 (MQ=255)                                                                                                                                                              aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                                                                                                                                                                       >  1:292739/1‑156 (MQ=255)                                                                                                                                                              aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                                                                                                                                                                       <  2:292739/156‑1 (MQ=255)                                                                                                                                                                             ttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGAGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCAGTACTCCACCACCATCCGTCGCAAGATGTCTTAATCCTTTGCCGCTTTCAGTTCGGTCTCGAACAATCGCTATTTGACACGTCATGCACAGTGAATGACGTGTTTGACAATGACAATGCGTTCGTATACCAggcgtctggaag                                                                                                      >  1:231124/1‑278 (MQ=255)                                                                                                                                                                                                            tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCgg                                                                    >  1:86810/1‑289 (MQ=255)                                                                                                                                                                                                            tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCgg                                                                    >  2:85839/1‑289 (MQ=255)                                                                                                                                                                                                                                                gTGACGAGGTGGTCTGCATTTCACCGGAATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGGTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGt                                 <  1:2979/288‑1 (MQ=255)                                                                                                                                                                                                                                                              ggCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGAAGGAAGAGGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCTAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAgg                  <  1:280948/288‑1 (MQ=255)                                                                                                                                                                                                                                                                          tgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCa      >  1:199695/1‑289 (MQ=255)                                                                                                                                                                                                                                                                              tgatAAGCAGGCTCTGCGGCGCCGCGTCTACTTTCCACGGTTCGCCCTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTGTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAGGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTg  <  1:421269/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAACAAGGCGTGTGTGACCAGtt   >  2:62619/1‑286 (MQ=255)                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                      TGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTTGCCAATGATATTGCGTTTGTAGAACAGGCGGCGTGAATACGGCGATTTGGGATACAGCGGCAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTG  >  NC_000913/2132507‑2133064

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