Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,502,713	C→T	L591L (CTG→CTA)	fryA ←	putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,502,713	0	C	T	100.0%	182.3 / NA	49	L591L (CTG→CTA)	fryA	putative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/24);  total (25/24)

GAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAGCACACC  >  NC_000913/2502568‑2502857                                                                                                                                                  |                                                                                                                                                 gAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGa                                                                                                                                              >  2:419082/1‑150 (MQ=255)          ggCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGt                                                                                                                                    >  2:805074/1‑151 (MQ=255)              actcatactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCAGTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATc                                                                                                                                <  2:185601/151‑1 (MQ=255)                      tcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTtat                                                                                                                         >  1:395393/1‑150 (MQ=255)                       caGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTtatat                                                                                                                       <  1:524167/151‑1 (MQ=255)                        aGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTtatat                                                                                                                       >  2:589940/1‑150 (MQ=255)                            cATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgc                                                                                                                    <  1:654993/149‑1 (MQ=255)                              tCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgcggcg                                                                                                                >  1:722417/1‑151 (MQ=255)                               ccAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgcggcg                                                                                                                <  1:37953/150‑1 (MQ=255)                                cAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGAt                                                                                                              <  1:619302/151‑1 (MQ=255)                                 aGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGAt                                                                                                              >  1:147010/1‑150 (MQ=255)                                 aGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGAt                                                                                                              >  2:396255/1‑150 (MQ=255)                                    cccagcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATAcgc                                                                                                          >  1:148439/1‑151 (MQ=255)                                    cccagcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATAcgc                                                                                                          <  2:543725/151‑1 (MQ=255)                                       agcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATAcgcgcg                                                                                                       >  2:317305/1‑151 (MQ=255)                                        gcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCgg                                                                                                      <  1:622662/151‑1 (MQ=255)                                         cccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCggg                                                                                                     >  1:909212/1‑151 (MQ=255)                                          ccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCggg                                                                                                     <  1:29672/150‑1 (MQ=255)                                            aaGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGtta                                                                                                  <  1:891641/151‑1 (MQ=255)                                            aaGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGtt                                                                                                   >  2:453175/1‑150 (MQ=255)                                              gCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGttatt                                                                                                <  2:564597/151‑1 (MQ=255)                                                          gATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGc                                                                                    >  1:620408/1‑151 (MQ=255)                                                              aCGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgcggcg                                                                                                                >  1:67279/1‑119 (MQ=255)                                                              aCGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgcggcg                                                                                                                <  2:67257/119‑1 (MQ=255)                                                               cGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACa                                                                               <  1:72228/151‑1 (MQ=255)                                                                                    ccGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTgg                                                          >  2:376959/1‑151 (MQ=255)                                                                                     cGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGa                                                         <  2:18515/151‑1 (MQ=255)                                                                                      gCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTgg                                                          >  1:487674/1‑149 (MQ=255)                                                                                       cAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGc                                                       >  2:390941/1‑151 (MQ=255)                                                                                        aaaTGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGcc                                                      >  1:897618/1‑151 (MQ=255)                                                                                            gCCTACCCATTTGCCCCGCTGATGCGCAGTGGTATCTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGAt                                                   <  1:627709/150‑1 (MQ=255)                                                                                                aCCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTg                                               >  2:595151/1‑150 (MQ=255)                                                                                                    aTTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTgaagaa                                          <  2:249322/151‑1 (MQ=255)                                                                                                        gCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACGGGGTCATATCGTTGGAGCCGATACTGAAGAAATc                                       <  2:871363/150‑1 (MQ=255)                                                                                                          cccGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGAc                                    >  1:69213/1‑151 (MQ=255)                                                                                                               tGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCAt                               >  1:896448/1‑151 (MQ=255)                                                                                                                   gcgcAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGc                            <  2:783492/150‑1 (MQ=255)                                                                                                                            tAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGa                                     >  2:213536/1‑132 (MQ=255)                                                                                                                             aaCTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCg                 >  1:812159/1‑151 (MQ=255)                                                                                                                                tATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCgatg              >  1:532571/1‑151 (MQ=255)                                                                                                                                   ttGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCgatgat            <  2:170192/150‑1 (MQ=255)                                                                                                                                    tgctgcAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGatgatg           <  2:498496/150‑1 (MQ=255)                                                                                                                                     gctgcAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGt          <  2:413980/150‑1 (MQ=255)                                                                                                                                     gctgcAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGt          <  1:854386/150‑1 (MQ=255)                                                                                                                                     gctgcAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTa         <  1:691436/151‑1 (MQ=255)                                                                                                                                        gcAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAg        <  1:267635/149‑1 (MQ=255)                                                                                                                                        gcAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAGc       >  2:858840/1‑150 (MQ=255)                                                                                                                                         cAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAGcac     >  1:519797/1‑151 (MQ=255)                                                                                                                                            cATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAGCACAcc  <  1:673866/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 GAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAGCACACC  >  NC_000913/2502568‑2502857

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