Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,542,969	A→G	F187S (TTC→TCC)	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,542,969	0	A	G	100.0%	154.4 / NA	42	F187S (TTC→TCC)	cysP	thiosulfate/sulfate ABC transporter periplasmic binding protein CysP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/14);  total (28/14)

TGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAG  >  NC_000913/2542837‑2543103                                                                                                                                     |                                                                                                                                       tGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGt                                                                                                                                 >  2:858495/1‑140 (MQ=255)   cgcTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCa                                                                                                                               >  2:616975/1‑140 (MQ=255)        cATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAAc                                                                                                                          >  2:790293/1‑140 (MQ=255)                  aCGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGttt                                                                                                                >  2:325067/1‑140 (MQ=255)                  aCGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCACGGCAAAAGTGGTGGACGCGCCACGACCGCCAGTATCGAACACATCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGttt                                                                                                                >  2:463106/1‑140 (MQ=255)                    ggATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTg                                                                                                              >  1:26169/1‑140 (MQ=255)                         ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCttt                                                                                                         >  2:247093/1‑140 (MQ=255)                         ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCttt                                                                                                         >  1:129281/1‑140 (MQ=255)                         ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCttt                                                                                                         >  1:1035049/1‑140 (MQ=255)                         ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCttt                                                                                                         >  1:385830/1‑140 (MQ=255)                         ttgttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCttt                                                                                                         >  2:834430/1‑140 (MQ=255)                            ttCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTc                                                                                                      >  1:749609/1‑140 (MQ=255)                                 ttCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTcaccac                                                                                                 >  1:121713/1‑140 (MQ=255)                                      ttccGAACCTAACCACCCCACCCCCCGGCCCGCCCCCGGAAAATGTGGTGTCCCCCCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCa                                                                                            <  2:305892/137‑1 (MQ=255)                                                   cAGCCCTTCGCCCAGGCCGCGCTCGGAAAAAGTGGTGGTCGCGCCACGGCCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTg                                                                               <  1:387152/140‑1 (MQ=255)                                                           cGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCAt                                                                       >  1:631239/1‑140 (MQ=255)                                                           cGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCAt                                                                       >  2:86544/1‑140 (MQ=255)                                                           cGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCAt                                                                       >  1:221922/1‑140 (MQ=255)                                                           cGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCAt                                                                       >  2:157152/1‑140 (MQ=255)                                                           cGCCCAGGCCGCGCTCGGCAAAAGTGGTGGCCGCCCCACGACCGCCAGTATCGAACACTTCACCGTTTTTCAGGGACTGGGCCATAAACGGTTCGGTTTTGCTTTTTCCACCCCCGTCAGCTTTATCCGCTGCCCCCCag                                                                       >  2:846173/1‑139 (MQ=255)                                                              ccAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATgcc                                                                    <  1:834430/140‑1 (MQ=255)                                                                  gCCGCGCTCGGCAAAGGTGGTGGCCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCa                                                                <  1:270904/140‑1 (MQ=255)                                                                  gCCGCGCTCGGAAAAGGGGGGGGTCGCCCCCCGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCa                                                                <  2:340174/140‑1 (MQ=255)                                                                    cgcgCTCGGCAAATGTGGGGGCCGCCCCGCTCCCGGCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGa                                                              <  2:93061/140‑1 (MQ=255)                                                                      cgcTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATa                                                            >  1:94140/1‑140 (MQ=255)                                                                            gCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGtata                                                      >  1:801054/1‑140 (MQ=255)                                                                                   gggtggtCGCCCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAAcgcgcg                                               <  1:616975/139‑1 (MQ=255)                                                                                     gtggtCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATTACGCGCGtt                                             >  2:994044/1‑140 (MQ=255)                                                                                         ccgccCCAGGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTAccc                                         <  2:121713/139‑1 (MQ=255)                                                                                                 cGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGtttt                                 >  1:305463/1‑140 (MQ=255)                                                                                                      gCCAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgggt                            <  2:631239/140‑1 (MQ=255)                                                                                                       ccAGTATCGAACACTTCAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgggtt                           >  1:763325/1‑140 (MQ=255)                                                                                                       ccAGTATCGAACACTTCAACGATTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTttcgggtt                           >  2:825812/1‑140 (MQ=255)                                                                                                                       cAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTc           >  1:833454/1‑140 (MQ=255)                                                                                                                       cAACGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTc           >  1:350240/1‑140 (MQ=255)                                                                                                                          cGTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACg        <  1:185600/140‑1 (MQ=255)                                                                                                                           gTTTTTCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGt       <  1:247093/140‑1 (MQ=255)                                                                                                                             ttttCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCGCGTCg     >  2:298285/1‑140 (MQ=255)                                                                                                                               ttCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGa   >  2:683100/1‑140 (MQ=255)                                                                                                                               ttCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGa   >  1:698684/1‑140 (MQ=255)                                                                                                                               ttCAGGGACTGGGTCATAAACTGTTCGGTTTTGCC‑TTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAg  <  1:858495/140‑1 (MQ=255)                                                                                                                                tCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAg  <  2:104657/140‑1 (MQ=255)                                                                                                                                tCAGGGACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAg  <  2:125423/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       TGCGCTTCATACTGTTTACGGATGTTGTTCACTTCCGATTCGAAGCTAATCAGCACATCGCCCAGGCCGCGCTCGGCAAAAGTGGTGGTCGCGCCACGACCGCCAGTATCGAACACTTCAACGTTTTTCAGGAACTGGGTCATAAACTGTTCGGTTTTGCCTTTGTCACCACCGTCAGCTTTATCCGCTGCGCCCCATGCCGCCAGATAGGTATAACGCGCGTTACCCGACGTTTTCGGGTTCGGGAAAATCAGCTTCACGTCGGAG  >  NC_000913/2542837‑2543103

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