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

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132650‑2132925                                                                                                                                          |                                                                                                                                           gcgggTAGAGGGTAAACAGCTGATATCTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTGGCGGTAAAAGTGGAGGGGCCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                       <  1:643436/140‑1 (MQ=255)      tAGAGGGTTAACAGCTGATATTTTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                  <  1:437136/143‑1 (MQ=255)           ggTAGACAGCTGATACTTTTTGTTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                             <  2:698306/143‑1 (MQ=255)           ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                             >  1:813416/1‑143 (MQ=255)            gTAAACGGCTTATAGTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAATGTGTAAGTACTCTGTGTATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            <  2:373083/143‑1 (MQ=255)                 cAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                       <  1:788753/143‑1 (MQ=255)                     tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTACTGGGTGACGAGGCGCTCTGCATTGCACCGTGATTGAGATGGAGGCTCTGCGTCGCCTCTTCTACTTACCa                                                                                                                   >  2:882135/1‑143 (MQ=255)                       ttatTTTTCGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGCGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                 <  1:681839/142‑1 (MQ=255)                                   ggCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTTATCTGCATTTCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTCGAGCTTTCCAGGGGTCGCCGGACt                                                                                                     >  2:512090/1‑143 (MQ=255)                                        ggTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                                <  1:491003/143‑1 (MQ=255)                                         gTCGAATTTGTCTTTGTGCTCGCGGTAAAAGTGGAAGTACCCGGAGCATTTTTTCTGGGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                               <  2:315033/143‑1 (MQ=255)                                           cGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATCTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCCACTTTCCAAGATTCTCCGTACTccaccacc                                                                                             >  2:694403/1‑143 (MQ=255)                                               tttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATcc                                                                                         >  2:175560/1‑143 (MQ=255)                                                   tCTTTGTGCTTGCTGTAAAAGTGGAAGGTCCCGGAGAGCGTCTTCTTTGTTGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                     <  2:638784/143‑1 (MQ=255)                                                         tgCTTGCGGTAAAAGTGGACGTACCGGGAGAATTTTGGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGGTCGCCGTACTCCACCACAATCCGCAGGAAggt                                                                               >  1:816151/1‑141 (MQ=255)                                                         tgCTTGCGGTAAAAGTGGAAGTACCCGGAAAATTTTTTCGGTGACGAGGTGATCAGCATTTGACCGTGGTTGATATGCAGACTCTGCGTCGCCTCTTCTACTTTCCACGGAACGCAGTACTCCACCACCATCCGCAGGAAgct                                                                               >  2:126225/1‑141 (MQ=255)                                                           cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaacacgt                                                                             >  2:28969/1‑138 (MQ=255)                                                           cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt                                                                             <  1:629848/143‑1 (MQ=255)                                                            ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            >  1:490100/1‑143 (MQ=255)                                                                ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACGATCCGCAGGAAGATGCCGTaa                                                                        >  2:729789/1‑143 (MQ=255)                                                                    aaggtcgAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTTTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt                                                                    <  1:404062/137‑1 (MQ=255)                                                                    aaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCTCCTTTCCCCGGTGCGCCGTCCTCCACCACCATCCGCAGCAAGATGTCGCAATCCt                                                                    >  2:679657/1‑143 (MQ=255)                                                                         tgAAATACCCGGAGAATCATTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGcc                                                               <  1:431873/142‑1 (MQ=255)                                                                         ggAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGTCTCTTCTACTGTCCACGGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGcc                                                               <  1:362439/143‑1 (MQ=255)                                                                           aaGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGc                                                             <  1:92755/143‑1 (MQ=255)                                                                             gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCACATGTCGTACTCCTGTGCCGCtt                                                           >  2:858661/1‑143 (MQ=255)                                                                               aCCCGGAGAATTTTTTTGGTGACGAGGTAATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTTGGTGTACTCCACCACCATCCGTATCAAGATTTTGTAATCCTTCGCCGCTTTc                                                         >  2:6089/1‑143 (MQ=255)                                                                                 ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTCCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAg                                                       <  1:779480/143‑1 (MQ=255)                                                                                    gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  2:748613/143‑1 (MQ=255)                                                                                      gaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                  >  1:106374/1‑143 (MQ=255)                                                                                      gaATTTTTGCGGTGAGGGGGTGATCTGCATTTCGCCGTGATTGGTATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                  <  2:178875/143‑1 (MQ=255)                                                                                         ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGGCGCTTTCAGTTCGGTaa                                               >  1:257266/1‑142 (MQ=255)                                                                                          tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                              <  1:673314/143‑1 (MQ=255)                                                                                           ttttCGGTGACGAGGTGATCTGCATTTCACCGTGACTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCACCACCATCCTCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTCCGGTATCg                                             >  1:320370/1‑143 (MQ=255)                                                                                             ttCGGGGACGAGGTGATCTGCATTTCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGGTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                           <  2:280467/143‑1 (MQ=255)                                                                                                 gTGGCGAGGTGATGTCCTTTCCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  2:825545/143‑1 (MQ=255)                                                                                                 gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGACGCTCTCACTTCGGTATCGaacaaa                                       >  2:24565/1‑143 (MQ=255)                                                                                                   gACGAGGTGCTCTGCACGTCACCGTGATTGATATGCAGGCTCTGTGATACCTCGTCTACTTACCGCGGCTTTCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                     <  1:169326/143‑1 (MQ=255)                                                                                                        ggtggtCTGCAGTGCTCCGTGATTGATATGCAGGCTCGGCGCCGCCTCTTCTACTTTGCACGGTTTGCCGTGCGCGACCACCATCCGCAGTAAGGTGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                <  2:668084/143‑1 (MQ=255)                                                                                                        ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCACGAAGATGTCGTAATCCTGTGCCCCTTTCAGTTCGGTATCGAACAAACACCCtt                                >  2:592727/1‑143 (MQ=255)                                                                                                                     tccccGTGATTGCTATGCTGGCTCTGCGTCGCCTCTTCTACTTTCCTCGGGTCGCTGTACTCCACCTCGATCTGCATTACGATGTCGTAATCCTTTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGc                   <  2:495325/140‑1 (MQ=255)                                                                                                                     tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGc                   <  2:329300/143‑1 (MQ=255)                                                                                                                         cGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg               <  2:813416/143‑1 (MQ=255)                                                                                                                         cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCGGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAGACGCCATGCCCAg               >  1:81521/1‑143 (MQ=255)                                                                                                                         cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCAACACCATCCGGGGGGAGATGTCGTAATCCTGTGCCGCTTTTCGTACGGTATCGAACAAACACTATTTCAACCGACATGtcccg               >  2:439331/1‑141 (MQ=255)                                                                                                                               tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCAACACCATCCACAGGACGCTGTCGTACTCCTCTGCCGCTTTCAGTTCGGTATCGACAAAACACTCTTTGACCCGCCATCCCCAGCGAAAg         >  2:710183/1‑143 (MQ=255)                                                                                                                                   atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGAACAGAccc     >  1:115928/1‑142 (MQ=255)                                                                                                                                      cAGGCTCTGCGTCGCCTCTGCGACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGACTGCTTACGCCTGTGCCGCTTTCAGCTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGt  <  1:209324/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                           TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132650‑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: ‑