Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,543,052	A→G	R159R (CGT→CGC)	cysP ←	thiosulfate/sulfate ABC transporter periplasmic binding protein CysP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,543,052	0	A	G	100.0%	112.8 / NA	33	R159R (CGT→CGC)	cysP	thiosulfate/sulfate ABC transporter periplasmic binding protein CysP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/17);  total (16/17)

TGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGC  >  NC_000913/2542919‑2543185                                                                                                                                      |                                                                                                                                      tggtggTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGcgcgcg                                                                                                                                 <  2:184673/140‑1 (MQ=255)     ggtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTAc                                                                                                                             <  1:127799/140‑1 (MQ=255)                   cGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTttcggg                                                                                                               <  1:97310/140‑1 (MQ=255)                               acacTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATc                                                                                                   <  2:677114/140‑1 (MQ=255)                                     cAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTc                                                                                             >  1:145271/1‑140 (MQ=255)                                                             tAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTc                                                                     <  2:668393/140‑1 (MQ=255)                                                             tAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTc                                                                     >  1:54241/1‑140 (MQ=255)                                                             tAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTc                                                                     >  2:669422/1‑140 (MQ=255)                                                                aCTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCaa                                                                  >  1:142270/1‑140 (MQ=255)                                                                    ttCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGt                                                              <  2:54241/140‑1 (MQ=255)                                                                    ttCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGt                                                              >  1:568436/1‑140 (MQ=255)                                                                    ttCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATACGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGt                                                              >  2:595041/1‑140 (MQ=255)                                                                      cGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAgg                                                                           <  1:11937/125‑1 (MQ=255)                                                                      cGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAgg                                                                           >  2:11937/1‑125 (MQ=255)                                                                      cGGTTTTGCC‑TTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGa                                                           <  1:405119/140‑1 (MQ=255)                                                                       ggTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGa                                                           <  2:568436/140‑1 (MQ=255)                                                                         ttttGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTGCACGTCGTAGCGCAGCAGGTCGTTCCAATCGTGGata                                                         >  1:225962/1‑140 (MQ=255)                                                                           ttGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAtt                                                       >  2:175280/1‑140 (MQ=255)                                                                             gCCTTTGTCACCACCATCAGCTTTAGCCGGTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttct                                                     <  1:308344/140‑1 (MQ=255)                                                                                  tGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGGTTCACGTCGGAGCGCACCAGGTCGTACCAAACGTGGATATTCTTCggg                                                >  2:308344/1‑140 (MQ=255)                                                                                   gTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGTCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttctt                                                    >  1:665427/1‑135 (MQ=255)                                                                                   gTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGTCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATAttctt                                                    <  2:665427/135‑1 (MQ=255)                                                                                             tCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACg                                     >  2:698961/1‑140 (MQ=255)                                                                                             tCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACg                                     >  1:502164/1‑140 (MQ=255)                                                                                              cAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGc                                    >  2:86661/1‑140 (MQ=255)                                                                                                    tATCCGCTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAgg                              <  1:576691/140‑1 (MQ=255)                                                                                                          cTGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGccc                        <  2:729670/140‑1 (MQ=255)                                                                                                           tGCGCCCCATGCCGCCAGATATGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGGTTCACGTCGTAGCGCAGCAGGTCGTTCCAATCGTGGATATGCTTCGGGTTACCCTTACGCACCAGAAAGCTCa                       >  2:457421/1‑140 (MQ=255)                                                                                                           tGCGCCCCATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCa                       <  2:502164/140‑1 (MQ=255)                                                                                                             cgcCCCATGCCGCCAGATAGGTATAGCGCGCGTTCCCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCAtg                     <  2:520584/140‑1 (MQ=255)                                                                                                                cccATGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCAtggtg                  <  1:698961/140‑1 (MQ=255)                                                                                                                   aTGCCGCCAGATAGGTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAg               >  2:165634/1‑140 (MQ=255)                                                                                                                                ggTATAGCGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGc  <  2:145271/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      TGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAGCGCACCAGGTCGTTCCAATCGTGGATATTCTTCGGGTTACCCTTACGCACCAGGAAGCCCATGGTGGAGTAGAACGGCGAGC  >  NC_000913/2542919‑2543185

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑