Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,091,007	C→T	A454T (GCG→ACG)	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,091,007	0	C	T	100.0%	193.7 / NA	56	A454T (GCG→ACG)	frvB	putative PTS enzyme IIBC component FrvB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/26);  total (30/26)

TCGCCGCTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGC  >  NC_000913/4090873‑4091142                                                                                                                                       |                                                                                                                                        tCGCCGCTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGc                                                                                                                                    >  1:1079532/1‑140 (MQ=255)    ccgcTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCaa                                                                                                                                 >  1:449513/1‑140 (MQ=255)       cTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCg                                                                                                                              >  1:1183813/1‑140 (MQ=255)                    cTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATa                                                                                                                 >  2:188725/1‑140 (MQ=255)                                 ttGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCAcc                                                                                                    >  2:847243/1‑140 (MQ=255)                                   gCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCgg                                                                                                  <  2:1226325/140‑1 (MQ=255)                                   gCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCgg                                                                                                  <  1:593820/140‑1 (MQ=255)                                   gCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCgg                                                                                                  <  2:280786/140‑1 (MQ=255)                                    cAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGc                                                                                                 <  1:984606/140‑1 (MQ=255)                                      gTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCgg                                                                                                  >  1:222459/1‑137 (MQ=255)                                      gTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCgg                                                                                                  <  2:222459/137‑1 (MQ=255)                                      gTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGa                                                                                               >  1:947663/1‑140 (MQ=255)                                        aaCTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACa                                                                                             >  1:1019089/1‑140 (MQ=255)                                        aaCTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACa                                                                                             >  1:908550/1‑140 (MQ=255)                                         aCTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACAt                                                                                            >  2:886795/1‑140 (MQ=255)                                         aCTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACAt                                                                                            >  1:1191953/1‑140 (MQ=255)                                           tCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCa                                                                                          <  2:685115/140‑1 (MQ=255)                                           tCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCa                                                                                          >  2:712172/1‑140 (MQ=255)                                              aTGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGc                                                                                       >  2:739276/1‑140 (MQ=255)                                                gTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCg                                                                                     >  1:445129/1‑140 (MQ=255)                                                 ttCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGa                                                                                    >  2:652724/1‑140 (MQ=255)                                                    aTGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCa                                                                                 >  1:213549/1‑140 (MQ=255)                                                     tGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAg                                                                                >  1:491983/1‑140 (MQ=255)                                                         cAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAg                                                                                <  2:85146/136‑1 (MQ=255)                                                         cAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAg                                                                                >  1:85146/1‑136 (MQ=255)                                                          aaGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGt                                                                           <  2:491983/140‑1 (MQ=255)                                                          aaGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGt                                                                           <  1:847243/140‑1 (MQ=255)                                                            gCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCa                                                                         >  1:481438/1‑140 (MQ=255)                                                               ccTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAg                                                                      <  2:1191953/140‑1 (MQ=255)                                                               ccTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAg                                                                      <  2:947663/140‑1 (MQ=255)                                                               ccTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAg                                                                      <  2:1183813/140‑1 (MQ=255)                                                               ccTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAg                                                                      <  1:652724/140‑1 (MQ=255)                                                                      cgcgGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCa                                                               >  2:447673/1‑140 (MQ=255)                                                                      cgcgGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCa                                                               <  1:566668/140‑1 (MQ=255)                                                                        cgGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCaaa                                                             >  1:1131794/1‑140 (MQ=255)                                                                         ggCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAAc                                                            <  1:712172/140‑1 (MQ=255)                                                                                                 caATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCACACCAGGTGCCAGACGTTTGATTCCGa                                    <  2:445129/140‑1 (MQ=255)                                                                                                      cTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCg                               <  2:1019089/140‑1 (MQ=255)                                                                                                      cTTTCAGGACGATAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCg                               <  1:188725/140‑1 (MQ=255)                                                                                                                 aTAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCAc                    <  2:908550/140‑1 (MQ=255)                                                                                                                 aTAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCAc                    <  1:447673/140‑1 (MQ=255)                                                                                                                  tAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGAt                             >  1:592776/1‑130 (MQ=255)                                                                                                                  tAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGAt                             <  2:592776/130‑1 (MQ=255)                                                                                                                  tAATGAATGAAATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACg                   >  2:134412/1‑140 (MQ=255)                                                                                                                       aatgaaATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGt              >  2:743419/1‑140 (MQ=255)                                                                                                                          gaaATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCa                     >  2:780428/1‑130 (MQ=255)                                                                                                                          gaaATGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCa                     <  1:780428/130‑1 (MQ=255)                                                                                                                           aaaTGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATg          >  1:257158/1‑140 (MQ=255)                                                                                                                            aaTGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCa                        >  1:228884/1‑125 (MQ=255)                                                                                                                            aaTGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCa                        <  2:228884/125‑1 (MQ=255)                                                                                                                            aaTGTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATgc         >  1:1032541/1‑140 (MQ=255)                                                                                                                               gTTCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATgccgc      >  1:622160/1‑140 (MQ=255)                                                                                                                                 tCAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGccgccc    >  1:226043/1‑140 (MQ=255)                                                                                                                                  cAGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGccgcccg   <  2:257158/140‑1 (MQ=255)                                                                                                                                   aGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCCCCGGCGACATCAGGCCGATCAGCGGGTAAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGccgcccgc  <  1:1274982/140‑1 (MQ=255)                                                                                                                                   aGCGTGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGccgcccgc  <  2:574818/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        TCGCCGCTGACGGCGCTGGCTTCGCACAACTGTTGCAGTAACTCAATGTTCATGGACAAGCTCCTGTTGCGCGGCTTTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGC  >  NC_000913/4090873‑4091142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑