Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,502,083	A→G	R801R (CGT→CGC)	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,083	0	A	G	100.0%	177.6 / NA	50	R801R (CGT→CGC)	fryA	putative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/33);  total (17/33)

GGCGTCAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTACGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATG  >  NC_000913/2501935‑2502223                                                                                                                                                     |                                                                                                                                             ggCGTCAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGcg                                                                                                                                            <  1:33856/151‑1 (MQ=255)    gTCAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTgcgcgc                                                                                                                                         >  2:572973/1‑151 (MQ=255)       aGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAg                                                                                                                                      <  1:108589/151‑1 (MQ=255)         tttGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCt                                                                                                                                    >  1:637033/1‑151 (MQ=255)          ttGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCt                                                                                                                                    <  2:340743/150‑1 (MQ=255)          ttGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCt                                                                                                                                    <  2:631020/150‑1 (MQ=255)                 gACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGagaa                                                                                                                             <  2:777171/150‑1 (MQ=255)                 gACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGagaa                                                                                                                             <  2:743403/150‑1 (MQ=255)                    gCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACt                                                                                                                         <  2:23192/151‑1 (MQ=255)                    gCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACt                                                                                                                         <  2:697555/151‑1 (MQ=255)                       gggAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTc                                                                                                                      <  2:504944/151‑1 (MQ=255)                         gagaGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTcac                                                                                                                    <  1:247898/151‑1 (MQ=255)                         gagaGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTcac                                                                                                                    >  1:310999/1‑151 (MQ=255)                         gagaGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTcac                                                                                                                    <  1:675473/151‑1 (MQ=255)                          agagTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTcaca                                                                                                                   <  2:596994/151‑1 (MQ=255)                           gagTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACAt                                                                                                                  >  2:55985/1‑151 (MQ=255)                            agTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATg                                                                                                                 <  1:339427/151‑1 (MQ=255)                             gTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGa                                                                                                                <  2:35210/151‑1 (MQ=255)                               aCAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGAt                                                                                                               <  1:271651/150‑1 (MQ=255)                                       cGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCtt                                                                                                      <  1:37486/151‑1 (MQ=255)                                           ttCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCTCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGtt                                                                                                  >  1:217563/1‑151 (MQ=255)                                            tCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGtt                                                                                                  <  1:194699/150‑1 (MQ=255)                                                 gCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCAcc                                                                                            <  2:251005/151‑1 (MQ=255)                                                  cTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCAccc                                                                                           >  2:322325/1‑151 (MQ=255)                                                         cAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCa                                                                                    <  1:724660/151‑1 (MQ=255)                                                          aGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTc                                                                                     >  1:197046/1‑149 (MQ=255)                                                                  aaTTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCa                                                                           <  2:185335/151‑1 (MQ=255)                                                                      cTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCAcc                                                                         >  1:744472/1‑149 (MQ=255)                                                                        gTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCa                                                                        <  1:461118/148‑1 (MQ=255)                                                                        gTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGtt                                                                     >  2:137828/1‑151 (MQ=255)                                                                         tttCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTc                                                                    >  2:44024/1‑151 (MQ=255)                                                                          ttCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGTCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTc                                                                    <  1:621512/150‑1 (MQ=255)                                                                            ccAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGGAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGa                                                                  <  2:724643/150‑1 (MQ=255)                                                                              agcagcGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACtt                                                               >  1:754513/1‑151 (MQ=255)                                                                                         gATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCt                                                    <  1:370889/151‑1 (MQ=255)                                                                                                   gCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTgccag                                           <  1:149949/150‑1 (MQ=255)                                                                                                    cATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTgccagcc                                         >  1:429588/1‑151 (MQ=255)                                                                                                          gcgcGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCg                                    <  1:432202/150‑1 (MQ=255)                                                                                                           cgcgGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGt                                  >  2:272409/1‑151 (MQ=255)                                                                                                            gcgGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGtt                                 <  2:268401/151‑1 (MQ=255)                                                                                                                  aacaacGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCa                           <  2:92656/151‑1 (MQ=255)                                                                                                                           tGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCa                   >  2:683067/1‑150 (MQ=255)                                                                                                                            ggCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGGCTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCa                   <  2:308002/149‑1 (MQ=255)                                                                                                                                gAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTg             <  1:741094/151‑1 (MQ=255)                                                                                                                                 aaGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTg             >  2:636376/1‑150 (MQ=255)                                                                                                                                   gTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCg                       >  2:331210/1‑138 (MQ=255)                                                                                                                                   gTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCg                       <  1:331177/138‑1 (MQ=255)                                                                                                                                   gTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGAt           <  1:546308/150‑1 (MQ=255)                                                                                                                                         aTTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGa     >  2:568198/1‑150 (MQ=255)                                                                                                                                           tCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATg  >  2:33067/1‑151 (MQ=255)                                                                                                                                            cACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATg  <  2:714959/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             GGCGTCAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTACGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATG  >  NC_000913/2501935‑2502223

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