Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,090,946	G→A	A474V (GCC→GTC)	frvB ←	putative PTS enzyme IIBC component FrvB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,090,946	0	G	A	100.0%	191.6 / NA	53	A474V (GCC→GTC)	frvB	putative PTS enzyme IIBC component FrvB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/25);  total (28/25)

TCATCGCCGCTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAG  >  NC_000913/4090870‑4091022                                                                             |                                                                             tCATCGCCGCTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACttt                                                                          >  2:826521/1‑81 (MQ=255)   aTCGCCGCTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTc                                                                         <  1:826460/80‑1 (MQ=255)   aTCGCCGCTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCg                                                                        >  2:826461/1‑81 (MQ=255)       ccgcTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCttt                                                                    <  2:826463/81‑1 (MQ=255)       ccgcTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCttt                                                                    <  2:826462/81‑1 (MQ=255)        cgcTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTg                                                                   >  2:826520/1‑81 (MQ=255)         gcTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgc                                                                  >  2:826545/1‑81 (MQ=255)         gcTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgc                                                                  >  1:826518/1‑81 (MQ=255)         gcTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgc                                                                  >  1:826522/1‑81 (MQ=255)          cTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgcg                                                                 >  2:826554/1‑81 (MQ=255)          cTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgcg                                                                 <  2:826464/81‑1 (MQ=255)            gACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgcggc                                                               >  2:826523/1‑81 (MQ=255)            gACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgcggc                                                               >  2:826543/1‑81 (MQ=255)             aCGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgcggcg                                                              >  2:826504/1‑81 (MQ=255)             aCGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTgcgg                                                                <  2:826465/79‑1 (MQ=255)              cGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGt                                                             <  2:826466/81‑1 (MQ=255)                       cTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAAt                                                    >  1:826507/1‑81 (MQ=255)                        ttCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATc                                                   >  1:826555/1‑81 (MQ=255)                         tCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATcc                                                  >  1:826538/1‑81 (MQ=255)                         tCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATcc                                                  >  2:826557/1‑81 (MQ=255)                            cacaACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCttt                                               >  2:826532/1‑81 (MQ=255)                            cacaACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCttt                                               <  2:826467/81‑1 (MQ=255)                             acaACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTc                                              <  2:826469/81‑1 (MQ=255)                             acaACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTc                                              <  2:826468/81‑1 (MQ=255)                                aCTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAg                                            >  1:826551/1‑80 (MQ=255)                                    ttGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGAc                                         <  1:826470/79‑1 (MQ=255)                                      gCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGAt                                      >  1:826564/1‑80 (MQ=255)                                      gCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGAt                                      <  2:826471/80‑1 (MQ=255)                                       cAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATa                                     <  1:826473/80‑1 (MQ=255)                                       cAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATa                                     <  1:826472/80‑1 (MQ=255)                                           aaCTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATaatgaa                                >  2:826535/1‑81 (MQ=255)                                            aCTCAATGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATaatgaat                               <  2:826474/81‑1 (MQ=255)                                                aaTGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAAtgaatgaa                            <  2:826475/80‑1 (MQ=255)                                                 aTGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAAt                          >  2:826566/1‑81 (MQ=255)                                                 aTGTTCATGGACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAAt                          <  2:826476/81‑1 (MQ=255)                                                          gACAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAgcgc                 <  2:826477/81‑1 (MQ=255)                                                            cAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAgcgc                 >  2:826548/1‑79 (MQ=255)                                                            cAAGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGa               <  1:826478/81‑1 (MQ=255)                                                              aGCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAgg             <  1:826479/81‑1 (MQ=255)                                                               gCTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGc            >  1:826556/1‑81 (MQ=255)                                                                cTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGcc           <  2:826481/81‑1 (MQ=255)                                                                cTCCTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGcc           <  2:826480/81‑1 (MQ=255)                                                                  ccTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCaa         >  2:826563/1‑81 (MQ=255)                                                                   cTGTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAAt        <  2:826482/81‑1 (MQ=255)                                                                     gTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCg      >  2:826565/1‑81 (MQ=255)                                                                     gTTGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCg      >  2:826559/1‑81 (MQ=255)                                                                      ttGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATc       <  2:826483/79‑1 (MQ=255)                                                                      ttGCGCGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGc     <  1:826484/81‑1 (MQ=255)                                                                        gcgcGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGcc    <  2:826485/80‑1 (MQ=255)                                                                        gcgcGACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGc     >  2:826549/1‑79 (MQ=255)                                                                          gcgACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGATTGAAATGTTCAGCGCGAGGCCAATCGCCAg  >  2:826525/1‑80 (MQ=255)                                                                          gcgACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAg  >  2:826561/1‑80 (MQ=255)                                                                          gcgACTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAg  >  2:826539/1‑80 (MQ=255)                                                                             |                                                                             TCATCGCCGCTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAG  >  NC_000913/4090870‑4091022

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