Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,541,456	A→G	I75T (ATT→ACT)	cysW ←	sulfate/thiosulfate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,541,456	0	A	G	100.0%	143.9 / NA	41	I75T (ATT→ACT)	cysW	sulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/20);  total (21/20)

CCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGC  >  NC_000913/2541361‑2541551                                                                                                |                                                                                                ccAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCAAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAg                                                                                             >  1:1573755/1‑100 (MQ=255)       gCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGat                                                                                       <  2:1573745/100‑1 (MQ=255)        cGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCCGTCAGCGCGatc                                                                                      <  1:1037791/100‑1 (MQ=255)           cagcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCCCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCCCCatcatc                                                                                   >  2:1265221/1‑100 (MQ=255)           cagcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatcatc                                                                                   >  2:18117/1‑100 (MQ=255)           cagcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatcatc                                                                                   >  2:1293385/1‑100 (MQ=255)                gTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCTGTACGGCAGTCAGCGCGATCATCACCg                                                                               <  1:478955/99‑1 (MQ=255)                       gCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAg                                                                       >  1:1128247/1‑100 (MQ=255)                          tCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATg                                                                    <  2:720618/100‑1 (MQ=255)                          tCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATg                                                                    >  1:490713/1‑100 (MQ=255)                               ggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTg                                                               >  1:382915/1‑100 (MQ=255)                                gAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgc                                                              >  1:43125/1‑100 (MQ=255)                                 aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                             >  2:235801/1‑100 (MQ=255)                                 aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                             >  2:235740/1‑100 (MQ=255)                                     tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                         >  1:1093413/1‑100 (MQ=255)                                      aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                                        <  1:1086512/100‑1 (MQ=255)                                          cgcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCg                                                    >  1:1055408/1‑100 (MQ=255)                                           gcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCgg                                                   <  2:1553861/100‑1 (MQ=255)                                            cgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGa                                                  <  1:235740/100‑1 (MQ=255)                                            cgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGa                                                  <  1:235801/100‑1 (MQ=255)                                                aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                              >  1:723305/1‑100 (MQ=255)                                                  cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                            >  2:276311/1‑100 (MQ=255)                                                    gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAg                                          >  1:469796/1‑100 (MQ=255)                                                         gccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCt                                     <  1:391557/100‑1 (MQ=255)                                                         gccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCt                                     <  1:658175/100‑1 (MQ=255)                                                          ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                    <  2:618285/100‑1 (MQ=255)                                                                gAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAc                              >  2:43077/1‑100 (MQ=255)                                                                    gCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                          <  1:963834/100‑1 (MQ=255)                                                                    gCCGAACACCAGGTTAACCAGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGAGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                          <  2:1340902/100‑1 (MQ=255)                                                                     ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                         >  1:1489700/1‑100 (MQ=255)                                                                              aGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCt                >  2:657934/1‑100 (MQ=255)                                                                                  ttACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCt            <  2:345439/100‑1 (MQ=255)                                                                                   tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg           >  1:1139645/1‑100 (MQ=255)                                                                                     ccGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGaa         <  2:1055402/100‑1 (MQ=255)                                                                                      cGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAAt        <  1:276311/100‑1 (MQ=255)                                                                                       ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg       <  1:18117/100‑1 (MQ=255)                                                                                       ggTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg       <  2:1046805/100‑1 (MQ=255)                                                                                         tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc     >  2:1096312/1‑100 (MQ=255)                                                                                         tACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc     <  1:647684/100‑1 (MQ=255)                                                                                           cGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGAACGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg   <  2:128589/100‑1 (MQ=255)                                                                                            ggCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGc  >  1:1337221/1‑100 (MQ=255)                                                                                                |                                                                                                CCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGC  >  NC_000913/2541361‑2541551

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑