Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,999,707	C→A	A227S (GCG→TCG)	tcyJ ←	cystine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,999,707	0	C	A	100.0%	105.7 / NA	31	A227S (GCG→TCG)	tcyJ	cystine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/15);  total (16/15)

ATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGCCACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATACGCCCTACGCGCAGAT  >  NC_000913/1999578‑1999832                                                                                                                                  |                                                                                                                              aTGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGcgacac                                                                                                                            >  1:372764/1‑133 (MQ=255)       attattTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGa                                                                                                                      <  1:87174/133‑1 (MQ=255)         tattTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACt                                                                                                                    >  1:567782/1‑133 (MQ=255)                 cacaTCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACgg                                                                                                            <  1:622012/133‑1 (MQ=255)                 cacaTCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGACCTCATTTACTTGACGGAGCGACACACCAGACTCCTGACgg                                                                                                            >  2:169138/1‑133 (MQ=255)                   caTCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGga                                                                                                          <  2:567782/133‑1 (MQ=255)                         cACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGc                                                                                                    >  1:6014/1‑133 (MQ=255)                         cACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTACTTTACGCAGCGACACGCCAGACTCCTGACGGGAGaa                                                                                                       >  2:610884/1‑130 (MQ=255)                         cACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTACTTTACGCAGCGACACGCCAGACTCCTGACGGGAGaa                                                                                                       <  1:610884/130‑1 (MQ=255)                              aaCCATTTTTTGGAAGGGGCTTGCGGAGAGCCCTCTTTTTGAATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCAc                                                                                               <  2:209503/133‑1 (MQ=255)                              aaCCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACGCGCCAGAGTCCTGATGGGAGAGTGCTTCAc                                                                                               >  1:209503/1‑133 (MQ=255)                              aaCCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCAc                                                                                               >  2:58657/1‑133 (MQ=255)                              aaCCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCAc                                                                                               >  1:570100/1‑133 (MQ=255)                              aaCCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCAc                                                                                               >  2:462594/1‑133 (MQ=255)                                                gCTTGCAGCGTGCCACCTATTTACATTTCCACAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTAt                                                                             <  1:193885/133‑1 (MQ=255)                                                      agagTGCCATCTTTTTGCATTTCCGCAATTGCATTATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGAGTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTgg                                                                       <  2:372764/133‑1 (MQ=255)                                                             cATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTttctt                                                                >  1:476356/1‑133 (MQ=255)                                                                cTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCAc                                                             <  1:58657/133‑1 (MQ=255)                                                                           tCCGCAATCGCATCAATCCCTGCTTTCAGCTGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAgcg                                                  <  1:169138/133‑1 (MQ=255)                                                                                  ttGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAg                                           <  2:476356/133‑1 (MQ=255)                                                                                    gCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGAc                                         <  2:108364/133‑1 (MQ=255)                                                                                                    tCAGCAGGTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGc                         >  1:155057/1‑133 (MQ=255)                                                                                                          ggTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGAt                   <  2:155057/133‑1 (MQ=255)                                                                                                           gTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATa                  >  1:131425/1‑133 (MQ=255)                                                                                                           gTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATa                  >  1:549117/1‑133 (MQ=255)                                                                                                           gTCCTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCCCCGCCGGCCGATCAACGAGGAAGGCATCGaca                  >  2:282226/1‑131 (MQ=255)                                                                                                              cTCATTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATACGc               >  2:315720/1‑133 (MQ=255)                                                                                                                 aTTTCCTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATACGCCCt            <  2:549117/133‑1 (MQ=255)                                                                                                                      cTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATACGCCCTAcgcg       <  2:64459/133‑1 (MQ=255)                                                                                                                      cTTTACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATACGCCCTAcgcg       >  1:138561/1‑133 (MQ=255)                                                                                                                       tttACGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATACGCCCTACgcgc      <  2:340415/133‑1 (MQ=255)                                                                                                                           cGCAGCGACACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATACGCCCTACGCGCAGAt  >  2:51959/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                              ATGCTGATTATTTGGTCACATCAGCACCAAACCATTTTTCGGAAAGGGCTTGCAGAGTGCCATCTTTTTGCATTTCCGCAATTGCATCATTCACTGCTTTCAGCAGGTCCTCATTTCCTTTACGCAGCGCCACGCCAGACTCCTGACGGGAGAATGCTTCACCGGTTACTGCCAGCGTATCGTTGGTTTTCTTCACCAGATCCAGCGCCGCCAGACGATCAACGAGGATCGCATCGATACGCCCTACGCGCAGAT  >  NC_000913/1999578‑1999832

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