Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,595,767	C→A	A213S (GCG→TCG)	livH ←	branched chain amino acid/phenylalanine ABC transporter membrane subunit LivH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,595,767	0	C	A	100.0%	113.9 / NA	33	A213S (GCG→TCG)	livH	branched chain amino acid/phenylalanine ABC transporter membrane subunit LivH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/17);  total (16/17)

CGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGCGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTCCGCG  >  NC_000913/3595689‑3595843                                                                               |                                                                             cGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCgagc                                                                            >  1:694719/1‑81 (MQ=255)   gCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGcc                                                                           >  2:694713/1‑80 (MQ=255)   gCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGcc                                                                           >  1:694710/1‑80 (MQ=255)             cGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAgg                                                                >  1:694728/1‑81 (MQ=255)                 gTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAg                                                            >  2:694737/1‑81 (MQ=255)                  tAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgc                                                           <  1:694675/81‑1 (MQ=255)                   aGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgcg                                                          >  2:694733/1‑81 (MQ=255)                    ggggTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgcgc                                                         <  2:694676/81‑1 (MQ=255)                            tGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCAccc                                                 <  1:694677/81‑1 (MQ=255)                             gACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcg                                                >  1:694723/1‑81 (MQ=255)                              aCGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcgg                                               >  2:694717/1‑81 (MQ=255)                                gCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtc                                             <  1:694678/81‑1 (MQ=255)                                  cGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcgg                                           >  2:694725/1‑81 (MQ=255)                                     aGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGgtgt                                        >  1:694731/1‑81 (MQ=255)                                      gAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCGCAAAGGTCAGCGCAATCACCCGGTCGGTGtt                                       <  1:694679/81‑1 (MQ=255)                                            aCCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATg                                   <  1:694682/79‑1 (MQ=255)                                            aCCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGcc                                 <  1:694681/81‑1 (MQ=255)                                            aCCGAGCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGc                                  <  1:694680/80‑1 (MQ=255)                                                agcagcACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAg                              >  2:694730/1‑80 (MQ=255)                                                  cagcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAg                              >  2:694727/1‑78 (MQ=255)                                                  cagcaCACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAg                           <  2:694683/81‑1 (MQ=255)                                                       cacCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCg                      <  1:694684/81‑1 (MQ=255)                                                           cGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCAt                  <  1:694685/81‑1 (MQ=255)                                                            gCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCAtt                 >  1:694735/1‑81 (MQ=255)                                                             ccACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCAttt                <  1:694686/81‑1 (MQ=255)                                                               aCCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCAtttt               >  1:694745/1‑80 (MQ=255)                                                               aCCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTc              <  2:694687/81‑1 (MQ=255)                                                                ccgccgCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCa             <  1:694688/81‑1 (MQ=255)                                                                  gccgccATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCAtttt               >  2:694736/1‑77 (MQ=255)                                                                      ccATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCtt       >  2:694763/1‑81 (MQ=255)                                                                      ccATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCtt       <  1:694689/81‑1 (MQ=255)                                                                       cATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTc      >  2:694724/1‑81 (MQ=255)                                                                        aTCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTcc     <  2:694691/81‑1 (MQ=255)                                                                        aTCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTcc     <  1:694690/81‑1 (MQ=255)                                                                           gccgaGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTCcgcg  >  2:694739/1‑81 (MQ=255)                                                                               |                                                                             CGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGCGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAAGCAGACTCGCCATTTTCAGATCTTCCGCG  >  NC_000913/3595689‑3595843

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