Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,540,599	A→G	R72R (CGT→CGC)	cysA ←	sulfate/thiosulfate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,540,599	0	A	G	100.0%	178.8 / NA	50	R72R (CGT→CGC)	cysA	sulfate/thiosulfate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/25);  total (25/25)

ACTTTCGCTTTGATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCACGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTG  >  NC_000913/2540465‑2540735                                                                                                                                       |                                                                                                                                         aCTTTCGCTTTGATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGcgtgc                                                                                                                                     >  1:593347/1‑140 (MQ=255)        cTTTGATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAgg                                                                                                                              <  1:107146/140‑1 (MQ=255)        cTTTGATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAgg                                                                                                                              <  1:403709/140‑1 (MQ=255)                 cTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACg                                                                                                                     <  2:593347/140‑1 (MQ=255)                   gCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTc                                                                                                                   >  1:675090/1‑140 (MQ=255)                                 cTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCTAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGc                                                                                                     >  2:866573/1‑140 (MQ=255)                                 cTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGc                                                                                                     >  2:766809/1‑140 (MQ=255)                                  tCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCg                                                                                                    <  1:698107/140‑1 (MQ=255)                                          cgGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGcc                                                                                            <  2:675090/140‑1 (MQ=255)                                           ggCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGccc                                                                                           <  1:305515/140‑1 (MQ=255)                                              agcaCCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCt                                                                                        >  1:317926/1‑140 (MQ=255)                                                     tCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGa                                                                                 >  1:726510/1‑140 (MQ=255)                                                          ccAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTc                                                                            <  2:922171/140‑1 (MQ=255)                                                             aaaGCGATATTGTCTAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAg                                                                         <  1:866573/140‑1 (MQ=255)                                                             aaaGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAg                                                                         <  1:766809/140‑1 (MQ=255)                                                                   atatTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATTGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCgg                                                                    <  2:709752/139‑1 (MQ=255)                                                                   atatTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCgg                                                                    >  1:709752/1‑139 (MQ=255)                                                                     atTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGa                                                                 <  1:435337/140‑1 (MQ=255)                                                                     atTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGa                                                                 <  2:324838/140‑1 (MQ=255)                                                                         tCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAAt                                                             >  2:771104/1‑140 (MQ=255)                                                                                ccGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCgcagc                                                      >  1:809155/1‑140 (MQ=255)                                                                                   tCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCagcagc                                                   >  1:513688/1‑140 (MQ=255)                                                                                          gCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGttt                                            <  1:771104/140‑1 (MQ=255)                                                                                           cGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGtttt                                           >  1:171289/1‑140 (MQ=255)                                                                                           cGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGtttt                                           >  1:176354/1‑140 (MQ=255)                                                                                           cGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGtttt                                           >  2:211952/1‑140 (MQ=255)                                                                                           cGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGtttt                                           >  1:275916/1‑140 (MQ=255)                                                                                           cGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGtttt                                           >  1:56895/1‑140 (MQ=255)                                                                                             gAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTcc                                         >  2:929771/1‑140 (MQ=255)                                                                                              aaCAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTccc                                        >  1:880599/1‑140 (MQ=255)                                                                                                   cgcgTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAAc                                   >  1:1002216/1‑140 (MQ=255)                                                                                                     cgTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCg                                 >  2:84051/1‑140 (MQ=255)                                                                                                     cgTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCg                                 >  2:360016/1‑140 (MQ=255)                                                                                                      gTAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCgg                                <  2:56895/140‑1 (MQ=255)                                                                                                       tAATGCTGGAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGa                               <  2:317926/140‑1 (MQ=255)                                                                                                              ggAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGcccc                        <  2:726510/140‑1 (MQ=255)                                                                                                               gAACACGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCa                       >  2:947299/1‑140 (MQ=255)                                                                                                               gAACACGAAACCG‑CTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAg                      >  2:352742/1‑140 (MQ=255)                                                                                                                    cGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACgg                            <  2:975953/130‑1 (MQ=255)                                                                                                                    cGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACgg                            <  2:183823/130‑1 (MQ=255)                                                                                                                    cGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACgg                            >  1:975953/1‑130 (MQ=255)                                                                                                                    cGAAACCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACgg                            >  1:183823/1‑130 (MQ=255)                                                                                                                       aaCCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAAcgcg               <  2:275916/140‑1 (MQ=255)                                                                                                                       aaCCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAAcgcg               <  2:176354/140‑1 (MQ=255)                                                                                                                       aaCCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAAcgcg               <  2:809155/140‑1 (MQ=255)                                                                                                                       aaCCGACTTTACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAAcgcg               >  1:99715/1‑140 (MQ=255)                                                                                                                              tttACGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCt        <  2:171289/140‑1 (MQ=255)                                                                                                                                 aCGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGAc     <  2:372803/140‑1 (MQ=255)                                                                                                                                 aCGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGAc     <  2:804717/140‑1 (MQ=255)                                                                                                                                 aCGATCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGAc     >  1:804717/1‑140 (MQ=255)                                                                                                                                    aTCGCGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTg  <  1:84051/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         ACTTTCGCTTTGATGGCTGCGGCATTCGGGCGCTCGCGACGCGGCAGCACCGTCAGGCCAAAAGCGATATTGTCGAACACCGTCATATGGCGGAACAGCGCGTAATGCTGGAACACGAAACCGACTTTACGATCACGTGCGTGCAGGCGGCTCACGTCGGTGCCGTGGAAGCGAATATGCCCGCTGGTTTGATGCTCCAGCCCGGCGATAATGCGCAGCAGCGTGGTTTTCCCGGAACCGGACGGCCCCAGCAACGCGACCATCTGACCTG  >  NC_000913/2540465‑2540735

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