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%	167.0 / NA	45	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 (24/21);  total (24/21)

GAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAGC  >  NC_000913/2502568‑2502852                                                                                                                                                  |                                                                                                                                            gAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGaa                                                                                                                                        <  2:516369/151‑1 (MQ=255) gAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGATCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAgg                                                                                                                                          >  2:765883/1‑149 (MQ=255)  aaTACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGATGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAAc                                                                                                                                       <  2:429639/151‑1 (MQ=255)     aCGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGa                                                                                                                                     <  2:338432/150‑1 (MQ=255)      cgcgGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATg                                                                                                                                   >  1:537159/1‑151 (MQ=255)         gggCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCg                                                                                                                                >  2:364582/1‑151 (MQ=255)              actcatactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAAt                                                                                                                            <  2:128502/150‑1 (MQ=255)              actcatactcaGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAAt                                                                                                                            <  2:351751/150‑1 (MQ=255)                       caGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTtatat                                                                                                                  <  1:137795/151‑1 (MQ=255)                       caGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTtata                                                                                                                   <  2:217595/150‑1 (MQ=255)                         gCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAg                                                                                                                <  2:239563/151‑1 (MQ=255)                             aTCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgcggc                                                                                                            >  1:17418/1‑151 (MQ=255)                              tCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgcggcg                                                                                                           >  2:658713/1‑151 (MQ=255)                              tCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgcggcg                                                                                                           >  2:14511/1‑151 (MQ=255)                               ccAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAgcggcg                                                                                                           <  1:568144/150‑1 (MQ=255)                                     ccagcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATAcgcg                                                                                                    <  2:619916/151‑1 (MQ=255)                                       agcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATAcgcgcg                                                                                                  >  2:360085/1‑151 (MQ=255)                                        gcccaAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCgg                                                                                                 <  2:17225/151‑1 (MQ=255)                                            aaGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGtta                                                                                             <  1:756452/151‑1 (MQ=255)                                             aGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGttat                                                                                            >  1:278180/1‑151 (MQ=255)                                              gCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGttatt                                                                                           <  1:322292/151‑1 (MQ=255)                                               cAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGttatta                                                                                          <  1:43993/151‑1 (MQ=255)                                                   aGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGAt                                                                                      <  1:449834/151‑1 (MQ=255)                                                         aGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCg                                                                                >  2:78756/1‑151 (MQ=255)                                                           aTAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGc                                                                               >  1:61813/1‑150 (MQ=255)                                                            tAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCCGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCAt                                                                             >  2:520666/1‑151 (MQ=255)                                                                  cTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATAc                                                                           >  1:666299/1‑147 (MQ=255)                                                                      cACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTg                                                                   >  2:124463/1‑151 (MQ=255)                                                                             cAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCat                                                             >  2:158077/1‑150 (MQ=255)                                                                                  cACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGtt                                                       <  2:561444/151‑1 (MQ=255)                                                                                        aaaTGCCTACCCATTTGCTCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGcc                                                 <  1:33036/151‑1 (MQ=255)                                                                                         aaTGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCg                                                <  2:458728/151‑1 (MQ=255)                                                                                           tGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGAGGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGAt                                              >  2:314724/1‑151 (MQ=255)                                                                                             ccTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATAc                                            <  1:714831/151‑1 (MQ=255)                                                                                              cTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGCTAc                                            >  2:176943/1‑150 (MQ=255)                                                                                                      ttGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAAt                                   >  2:535149/1‑151 (MQ=255)                                                                                                       tGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGaaa                                    <  1:216591/149‑1 (MQ=255)                                                                                                           ccGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGAcc                              >  2:307237/1‑151 (MQ=255)                                                                                                             gCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTc                            >  2:70019/1‑151 (MQ=255)                                                                                                                gATGCGCAGTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATc                         >  2:586126/1‑151 (MQ=255)                                                                                                                        gTGGTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGt                 >  1:659888/1‑151 (MQ=255)                                                                                                                           gTAACTATTTGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGt              >  2:680162/1‑151 (MQ=255)                                                                                                                               ctctttGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGGTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCCTCTCAGCAGTGGTCgat          >  1:661452/4‑151 (MQ=255)                                                                                                                                   ttGCTGCAACATGCGTAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCgatgat       <  1:482358/150‑1 (MQ=255)                                                                                                                                        gcAACATGCGTAGGAACGATGGTGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAGc  >  2:774143/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                            GAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCACCGCCCAGTTCACCGCAAATGCCTACCCATTTGCCCCGCTGATGCGCAGTGGTAACTATTTGCTGCAACATGCGCAGGAACGATGGCGTAATCGGGTTATATAGCGGCGATACGCGCGGGTTATTACGATCGACCGCATACAGATACTGGGTCATATCGTTGGAGCCGATACTGAAGAAATCGACCTCATCGCAGAAGTGGTCGATGATGTAGC  >  NC_000913/2502568‑2502852

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