Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,595,737	G→T	L223I (CTC→ATC)	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,737	0	G	T	100.0%	167.5 / NA	47	L223I (CTC→ATC)	livH	branched chain amino acid/phenylalanine ABC transporter membrane subunit LivH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (26/21);  total (26/21)

CGAGCACCGCCGCGGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGCGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAA  >  NC_000913/3595659‑3595814                                                                               |                                                                              cGAGCACCGCCGCGGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGatca                                                                             <  1:628171/81‑1 (MQ=255)    gCACCGCCGCGGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAtcagca                                                                          >  1:628229/1‑81 (MQ=255)     cACCGCCGCGGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAtcagca                                                                          >  1:628233/1‑80 (MQ=255)      aCCGCCGCGGTAAAGGCTTTCATCCCGGCCATAAGGCCGATGTAGGGGTTAATGACGCCGTAGAACTGCCCGATCAGcaca                                                                        <  2:628171/81‑1 (MQ=255)         gccgcGGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACAccc                                                                     >  2:628254/1‑81 (MQ=255)          ccgcGGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCg                                                                    <  2:628172/81‑1 (MQ=255)           cgcgGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGc                                                                   >  2:628243/1‑81 (MQ=255)           cgcgGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGc                                                                   >  1:628221/1‑81 (MQ=255)            gcgGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGcc                                                                  <  2:628173/81‑1 (MQ=255)             cgGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCa                                                                 <  1:628176/81‑1 (MQ=255)             cgGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCa                                                                 <  2:628174/81‑1 (MQ=255)             cgGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCa                                                                 >  1:628228/1‑81 (MQ=255)              ggTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCAc                                                                >  2:628219/1‑81 (MQ=255)                  aaGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCAccgcc                                                            >  1:628227/1‑81 (MQ=255)                  aaGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCAccgcc                                                            >  2:628222/1‑81 (MQ=255)                    ggCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACcgccgc                                                          >  2:628245/1‑81 (MQ=255)                       tttCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCAt                                                       >  2:628235/1‑81 (MQ=255)                         tCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATc                                                      <  1:628177/80‑1 (MQ=255)                           aTCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATcgcc                                                   >  2:628251/1‑81 (MQ=255)                           aTCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATcgc                                                    <  2:628178/80‑1 (MQ=255)                           aTCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATcgc                                                    <  2:628177/80‑1 (MQ=255)                            tCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATcgccg                                                  <  1:628180/81‑1 (MQ=255)                              ccGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATcgccga                                                 <  2:628180/80‑2 (MQ=255)                                ggCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGcc                                              <  1:628182/81‑1 (MQ=255)                                   cATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAAt                                           >  1:628231/1‑81 (MQ=255)                                    aTAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATc                                          <  2:628182/81‑1 (MQ=255)                                      aaaGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATca                                         >  2:628229/1‑80 (MQ=255)                                             aTGTAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTc                                 <  2:628183/81‑1 (MQ=255)                                                tAGGGGTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgc                              >  1:628237/1‑81 (MQ=255)                                                  ggggTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgcgc                            <  2:628184/81‑1 (MQ=255)                                                  ggggTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAgcgc                            <  2:628185/81‑1 (MQ=255)                                                     gTTGATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAAt                         >  2:628237/1‑81 (MQ=255)                                                     gTTAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAAt                         >  1:628235/1‑81 (MQ=255)                                                      ttAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAAt                         >  1:628242/1‑80 (MQ=255)                                                      ttAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATc                        <  1:628187/81‑1 (MQ=255)                                                       tAATGACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCa                       >  2:628256/1‑81 (MQ=255)                                                           gACGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcg                   >  1:628256/1‑81 (MQ=255)                                                             cGCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggt                 >  2:628253/1‑81 (MQ=255)                                                              gCCGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtc                >  1:628250/1‑81 (MQ=255)                                                               ccGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcg               <  2:628187/81‑1 (MQ=255)                                                                cGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcgg              >  1:628239/1‑81 (MQ=255)                                                                cGTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcgg              >  1:628243/1‑81 (MQ=255)                                                                 gTAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCcggtcggt             >  1:628234/1‑81 (MQ=255)                                                                  tAGAACTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGTGCAATCACCCGGTCGgtg            <  1:628189/81‑1 (MQ=255)                                                                     aaCTGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTa         <  2:628189/81‑1 (MQ=255)                                                                        tGACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATg      >  1:628240/1‑81 (MQ=255)                                                                         gACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGc     >  2:628250/1‑81 (MQ=255)                                                                         gACCGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGc     >  2:628240/1‑81 (MQ=255)                                                                           ccGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCa   <  1:628191/81‑1 (MQ=255)                                                                            cGATCAGCACACCCGCCACCGCCGCCATCGCCGAGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCaa  <  1:628192/81‑1 (MQ=255)                                                                               |                                                                              CGAGCACCGCCGCGGTAAAGGCTTTCATCCCGGCCATAAAGCCGATGTAGGGGTTAATGACGCCGTAGAACTGACCGAGCAGCACACCCGCCACCGCCGCCATCGCCGCGCCAATCACAAAGGTCAGCGCAATCACCCGGTCGGTGTTAATGCCAA  >  NC_000913/3595659‑3595814

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