Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,504,391	T→C	E32G (GAA→GGA)	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,504,391	0	T	C	100.0%	170.5 / NA	49	E32G (GAA→GGA)	fryA	putative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/23);  total (26/23)

CATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTTCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCAGGAGATCATCTGAAATCCGT  >  NC_000913/2504272‑2504524                                                                                                                        |                                                                                                                                      cATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCa                                                                                                                    <  1:800352/139‑1 (MQ=255)  aTCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                   <  1:696907/139‑1 (MQ=255)  aTCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAAGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTGCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                   <  2:561504/139‑1 (MQ=255)                     aGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACt                                                                                                                           <  2:960701/112‑1 (MQ=255)                     aGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACt                                                                                                                           >  1:960701/1‑112 (MQ=255)                            aaCTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACg                                                                                         <  1:304229/139‑1 (MQ=255)                              cTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAg                                                                                       <  2:279245/139‑1 (MQ=255)                                 tCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGa                                                                                           <  1:444947/132‑1 (MQ=255)                                 tCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGa                                                                                           >  2:444947/1‑132 (MQ=255)                                   aTTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGt                                                                                  >  2:251314/1‑139 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAGTAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGc                                                                                      >  2:719269/1‑127 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                   <  2:511281/98‑1 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                   >  1:511281/1‑98 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAg                                                                                       >  1:474296/1‑126 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAg                                                                                       <  2:474296/126‑1 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGc                                                                                      <  1:719269/127‑1 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGtt                                                                          <  1:626107/139‑1 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGtt                                                                          >  2:328634/1‑139 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGtt                                                                          >  2:872292/1‑139 (MQ=255)                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGtt                                                                          >  2:69015/1‑139 (MQ=255)                                            ggtgccggtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGc                                                                                      <  2:138747/126‑1 (MQ=255)                                            ggtgccggtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGc                                                                                      >  1:138747/1‑126 (MQ=255)                                            ggtgccggtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCATTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTa                                                                         >  2:769434/1‑139 (MQ=255)                                             gtgccggtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCTACGCCGGACGAGCATGTAGACCGTTAg                                                                        >  2:111721/1‑139 (MQ=255)                                             gtgccggtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAg                                                                        >  2:340906/1‑139 (MQ=255)                                             gtgccggtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAg                                                                        >  2:252489/1‑139 (MQ=255)                                                ccggtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCa                                                                     >  2:180237/1‑139 (MQ=255)                                                ccggtcgcAATCATCGCCAGCGAGTGTTTCGGCTCTGCTTTTGCGTTCGGGCTATTGTTAATTAAAGTGATTCCTCTTTGCCGCTGGCTGTACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCa                                                                     <  2:896743/139‑1 (MQ=255)                                                   gtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTTATAAAAGTGATTCCGCt                                                                                                                                   <  1:660231/74‑1 (MQ=255)                                                   gtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCt                                                                                                                                   >  2:660231/1‑74 (MQ=255)                                                      ccAATCAGCTCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGAc                                                               >  1:707854/1‑139 (MQ=255)                                                      ccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGAc                                                               >  2:734381/1‑139 (MQ=255)                                                              cgcCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAAt                                                       >  2:196524/1‑139 (MQ=255)                                                                ccAGCGAGCTTTTCGCATCTTCCTTTGCGTTCTGGCGAGGGTTAATGAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCTTTAGGCAAAGGACAGAGAAATTg                                                     <  2:732248/139‑1 (MQ=255)                                                                ccAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTg                                                     <  1:584647/139‑1 (MQ=255)                                                                 cAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGa                                                    <  2:707854/139‑1 (MQ=255)                                                                      aGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGt                                               <  1:872292/139‑1 (MQ=255)                                                                        cTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTATTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTa                                             >  1:732248/1‑139 (MQ=255)                                                                                   tGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCt                                  <  1:83254/139‑1 (MQ=255)                                                                                     ccTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCa                                                                                                  >  2:246151/1‑73 (MQ=255)                                                                                     ccTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCa                                                                                                  <  1:246151/73‑1 (MQ=255)                                                                                         tGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGt                            <  1:196524/139‑1 (MQ=255)                                                                                         tGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGt                            <  1:2168/139‑1 (MQ=255)                                                                                               cTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCa                      >  2:951383/1‑139 (MQ=255)                                                                                               cTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCa                      >  1:903378/1‑139 (MQ=255)                                                                                                     aTGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACg                                                                                         >  2:636056/1‑66 (MQ=255)                                                                                                     aTGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACg                                                                                         <  1:636056/66‑1 (MQ=255)                                                                                                                  gTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACg                                                                                         >  2:195075/1‑53 (MQ=255)                                                                                                                  gTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACg                                                                                         <  1:195075/53‑1 (MQ=255)                                                                                                                   tGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCAGGAGATCATCTGAAATCCGt  >  2:914725/1‑139 (MQ=255)                                                                                                                        |                                                                                                                                      CATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTTCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCAGGAGATCATCTGAAATCCGT  >  NC_000913/2504272‑2504524

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