Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,541,004	C→T	G226S (GGC→AGC)	cysW ←	sulfate/thiosulfate ABC transporter inner membrane subunit CysW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,541,004	0	C	T	100.0%	198.0 / NA	57	G226S (GGC→AGC)	cysW	sulfate/thiosulfate ABC transporter inner membrane subunit CysW
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/25);  total (32/25)

CATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGC  >  NC_000913/2540867‑2541134                                                                                                                                          |                                                                                                                                   cATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTgg                                                                                                                                  <  2:1237659/140‑1 (MQ=255)    aCTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGaaa                                                                                                                               <  1:204437/140‑1 (MQ=255)      tttttAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAacc                                                                                                                             >  1:839411/1‑140 (MQ=255)                  gggTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACcgcg                                                                                                                 <  1:116135/140‑1 (MQ=255)                  gggTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACcgcg                                                                                                                 >  1:68505/1‑140 (MQ=255)                  gggTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACcgcg                                                                                                                 >  1:194748/1‑140 (MQ=255)                  gggTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACcgcg                                                                                                                 >  1:399586/1‑140 (MQ=255)                  gggTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCCCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAACCGAGCTGGAAACCGCCGACACcgcg                                                                                                                 <  2:495537/140‑1 (MQ=255)                             gCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGcc                                                                                                      >  2:365319/1‑140 (MQ=255)                                 tCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAAtt                                                                                                  >  2:708628/1‑140 (MQ=255)                                        tAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCggg                                                                                           >  1:658652/1‑140 (MQ=255)                                          aCAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         <  1:1160528/140‑1 (MQ=255)                                          aCAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         <  2:695617/140‑1 (MQ=255)                                          aCAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         <  2:658652/140‑1 (MQ=255)                                          aCAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         <  2:109678/140‑1 (MQ=255)                                             gcgcCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         >  1:773524/1‑137 (MQ=255)                                             gcgcCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         <  2:1259430/137‑1 (MQ=255)                                             gcgcCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         >  1:1259430/1‑137 (MQ=255)                                             gcgcCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         <  2:773524/137‑1 (MQ=255)                                             gcgcCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTg                                                                                      >  2:551110/1‑140 (MQ=255)                                             gcgcCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTg                                                                                      >  2:553651/1‑140 (MQ=255)                                              cgccgcAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTgg                                                                                     >  2:995285/1‑140 (MQ=255)                                               gccgcAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         <  1:1030141/135‑1 (MQ=255)                                               gccgcAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCg                                                                                         >  2:1030141/1‑135 (MQ=255)                                               gccgcAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGt                                                                                    >  1:192681/1‑140 (MQ=255)                                                       ggTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAAcacca                                                                            <  1:319502/140‑1 (MQ=255)                                                       ggTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAAcacca                                                                            <  2:11333/140‑1 (MQ=255)                                                         tAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACg                                                                          <  2:1045480/140‑1 (MQ=255)                                                              gAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATa                                                                     <  2:426152/140‑1 (MQ=255)                                                                gCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATaaa                                                                   >  2:1158566/1‑140 (MQ=255)                                                                 ccGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAag                                                                  >  1:635536/1‑140 (MQ=255)                                                                        tgttgtAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGcc                                                           >  2:491293/1‑140 (MQ=255)                                                                              aGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCgg                                                     <  1:971777/140‑1 (MQ=255)                                                                                  cTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGc                                                       <  2:826147/134‑1 (MQ=255)                                                                                  cTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGc                                                       >  1:826147/1‑134 (MQ=255)                                                                                  cTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGt                                                 >  1:871111/1‑140 (MQ=255)                                                                                  cTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGt                                                 >  2:971681/1‑140 (MQ=255)                                                                                  cTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGt                                                 >  2:1198707/1‑140 (MQ=255)                                                                                  cTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGt                                                 >  1:348017/1‑140 (MQ=255)                                                                                   tGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGtt                                                <  2:635536/140‑1 (MQ=255)                                                                                    gCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTc                                               >  1:936816/1‑140 (MQ=255)                                                                                           cAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAAtg                                        <  2:348017/140‑1 (MQ=255)                                                                                                cAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGAcg                                   <  1:1200409/140‑1 (MQ=255)                                                                                                  aTCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGAcgcg                                 >  1:767989/1‑140 (MQ=255)                                                                                                    cTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGAc                               <  2:1239675/140‑1 (MQ=255)                                                                                                     tGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACg                              <  1:551110/140‑1 (MQ=255)                                                                                                     tGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACg                              <  2:767989/140‑1 (MQ=255)                                                                                                         aCGGCAGCGACAGGGTTTCGCCTAGAATCGAGCTGGAAGCCACCGACACCGCGCCAAACCCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAAc                          <  2:839411/140‑1 (MQ=255)                                                                                                             cAGCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCt                      >  1:305266/1‑140 (MQ=255)                                                                                                               gCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGc                    >  1:675557/1‑140 (MQ=255)                                                                                                               gCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGc                    >  1:1269477/1‑140 (MQ=255)                                                                                                               gCGACAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGc                    >  1:776599/1‑140 (MQ=255)                                                                                                                  aCAGGGTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAg                 <  2:871111/140‑1 (MQ=255)                                                                                                                      ggTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCgg             >  1:8689/1‑140 (MQ=255)                                                                                                                       gTTTCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGa            <  1:775502/140‑1 (MQ=255)                                                                                                                        tttCGCCGCGAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGAc           >  1:985494/1‑140 (MQ=255)                                                                                                                                cgAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAg   >  1:24521/1‑140 (MQ=255)                                                                                                                                 gAATCGAGCTGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGc  >  1:365251/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                   CATACTTTTTAAAAACAGGGTGATAATCGCCATCAGCGTTAACAGCGCCGCAGCGGTAAAGGAGCCGACGGTGTTGTAGTCCTGCTCCAGCAATTCAATCTGTAACGGCAGCGACAGGGTTTCGCCGCGAATCGAGCCGGAAACCACCGACACCGCGCCAAACTCGCCAATTGCGCGGGCGTTGGTCAACACCACGCCATAAAGCAGCGCCCAGCGGATGTTCGGTAATGTGACGCGACGGAACATCTGCCAGCCGGACGCGCCAAGC  >  NC_000913/2540867‑2541134

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: ‑