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%	181.5 / NA	51	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 (24/27);  total (24/27)

CGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTTCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCAGGAGATCATCTGAAATCCGTCAG  >  NC_000913/2504256‑2504527                                                                                                                                        |                                                                                                                                         cgccgcgccTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGc                                                                                                                                       >  2:277343/1‑139 (MQ=255)  gccgcgccTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCt                                                                                                                                      <  1:409490/139‑1 (MQ=255)       gccTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGcc                                                                                                                                 <  2:592133/139‑1 (MQ=255)         cTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCAc                                                                                                                               <  2:391424/139‑1 (MQ=255)            ctctTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTgg                                                                                                                            >  1:474535/1‑139 (MQ=255)                  aTCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                      <  2:474535/139‑1 (MQ=255)                      cTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                      >  1:101362/1‑135 (MQ=255)                      cTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                      <  2:101362/135‑1 (MQ=255)                      cTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTggt                                                                                                                  >  1:445983/1‑137 (MQ=255)                      cTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGtt                                                                                                                  <  2:427997/139‑1 (MQ=255)                                     aGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACg                                                                                                   >  2:559735/1‑139 (MQ=255)                                       gCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTAGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGcc                                                                                                 >  2:737096/1‑139 (MQ=255)                                        cTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCg                                                                                                <  1:737096/139‑1 (MQ=255)                                            aaCTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACg                                                                                            <  1:168666/139‑1 (MQ=255)                                                 tCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCAt                                                                                       <  1:559735/139‑1 (MQ=255)                                                    ttAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTa                                                                                    <  1:193849/139‑1 (MQ=255)                                                      aaaTAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGGATCTGCCTTTGCGTTCTGGCGATGGTTATTAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGa                                                                                  >  1:255600/1‑139 (MQ=255)                                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCt                                                                                                                >  1:849610/1‑104 (MQ=255)                                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCt                                                                                                                <  2:849610/104‑1 (MQ=255)                                                           gggTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGtt                                                                             >  1:19109/1‑139 (MQ=255)                                                               gccggtgccAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGc                                                                         >  1:184129/1‑139 (MQ=255)                                                                       cAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACa                                                                 <  2:19109/139‑1 (MQ=255)                                                                              cgcCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACt                                                                                                                     <  1:857958/80‑1 (MQ=255)                                                                              cgcCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACt                                                                                                                     >  2:857958/1‑80 (MQ=255)                                                                              cgcCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAAt                                                          >  2:1215810/1‑139 (MQ=255)                                                                               gcCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAAtt                                                         <  2:184129/139‑1 (MQ=255)                                                                                    cGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACt                                                                                                                     >  1:587106/1‑74 (MQ=255)                                                                                    cGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACt                                                                                                                     <  2:587106/74‑1 (MQ=255)                                                                                    cGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATc                                                    <  2:846185/139‑1 (MQ=255)                                                                                        cTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTa                                                >  2:473993/1‑139 (MQ=255)                                                                                              gCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATaa                                          <  1:473993/139‑1 (MQ=255)                                                                                                aTCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAAcc                                        <  2:255600/139‑1 (MQ=255)                                                                                                   tGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCa                                                                                                                       >  2:414592/1‑57 (MQ=255)                                                                                                   tGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCa                                                                                                                       <  1:414592/57‑1 (MQ=255)                                                                                                     ccTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                      >  2:60588/1‑56 (MQ=255)                                                                                                     ccTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                      <  1:60588/56‑1 (MQ=255)                                                                                                        ttGCGTTCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCggg                                <  1:1031234/139‑1 (MQ=255)                                                                                                             ttCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTa                                                                                                             <  2:1008157/57‑1 (MQ=255)                                                                                                             ttCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTa                                                                                                             >  1:1008157/1‑57 (MQ=255)                                                                                                              tCTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATc                          <  1:333215/139‑1 (MQ=255)                                                                                                               cTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGt                               >  1:851889/1‑133 (MQ=255)                                                                                                               cTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGt                               <  2:851889/133‑1 (MQ=255)                                                                                                               cTGGCGATGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCa                         >  2:721584/1‑139 (MQ=255)                                                                                                                     aTGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCa                                                                                                                       >  2:212270/1‑39 (MQ=255)                                                                                                                     aTGGTTAATAAAAGTGATTCCGCTTTGCCACTGGCTGCa                                                                                                                       <  1:212270/39‑1 (MQ=255)                                                                                                                          taataaAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                      <  2:400083/35‑1 (MQ=255)                                                                                                                          taataaAAGTGATTCCGCTTTGCCACTGGCTGCAc                                                                                                                      >  1:400083/1‑35 (MQ=255)                                                                                                                                   tGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACg                                                                                            <  1:868769/52‑1 (MQ=255)                                                                                                                                   tGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACg                                                                                            >  2:868769/1‑52 (MQ=255)                                                                                                                                   tGATTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCAGGAGATCATCTGAAATCCGt     >  1:77429/1‑139 (MQ=255)                                                                                                                                     aTTCCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATATCCGCTCCGGGTTAATCAGGAG‑TCATCTGAAATCCGTCAg  >  2:852889/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                         CGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTCGCATCTGCCTTTGCGTTCTGGCGATGGTTAATAAAAGTGATTTCGCTTTGCCACTGGCTGCACTGTTCTTTAAGTTCCCACGCCGGACGAGCATGTAGACCGTTAGGCAAAGGACAGAGAAATTGAATCGTTAACATAACCGCTCCGGGTTAATCAGGAGATCATCTGAAATCCGTCAG  >  NC_000913/2504256‑2504527

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