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%	143.4 / NA	40	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 (13/27);  total (13/27)

GTCAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTACGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATG  >  NC_000913/2501938‑2502223                                                                                                                                                  |                                                                                                                                             gTCAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTGTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTgcgcgc                                                                                                                                         <  2:79065/151‑1 (MQ=255)   cAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCa                                                                                                                                       <  1:110844/151‑1 (MQ=255)   cAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCa                                                                                                                                       <  2:260293/151‑1 (MQ=255)    aGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAg                                                                                                                                      <  2:819964/151‑1 (MQ=255)                cGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGAc                                                                                                                          >  2:445940/1‑151 (MQ=255)                            aCAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGAtt                                                                                                              <  2:762664/151‑1 (MQ=255)                            aCAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGAtt                                                                                                              <  1:561200/151‑1 (MQ=255)                              aGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCa                                                                                                            <  2:601040/151‑1 (MQ=255)                                cgcCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATc                                                                                                          <  2:765006/151‑1 (MQ=255)                                   ccGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCt                                                                                                       >  2:446630/1‑151 (MQ=255)                                      ttttCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTc                                                                                                     <  1:423962/150‑1 (MQ=255)                                        ttCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGtt                                                                                                  >  1:401844/1‑151 (MQ=255)                                          cAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTa                                                                                                 <  2:793091/150‑1 (MQ=255)                                             ggCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCAc                                                                                             >  1:599459/1‑151 (MQ=255)                                                                       ttCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTc                                                                    <  2:271232/150‑1 (MQ=255)                                                                         ccAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGa                                                                  <  2:715011/150‑1 (MQ=255)                                                                          cagcagCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACt                                                                >  1:207920/1‑151 (MQ=255)                                                                            gcagcGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTc                                                              <  1:766413/151‑1 (MQ=255)                                                                                   cAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCAtt                                                       <  2:239642/151‑1 (MQ=255)                                                                                    aGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCAttt                                                      >  2:838689/1‑151 (MQ=255)                                                                                     ggATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTc                                                     <  1:181673/151‑1 (MQ=255)                                                                                         aCTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGa                                                  >  1:817376/1‑150 (MQ=255)                                                                                              gtgCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTgcca                                            <  2:638404/151‑1 (MQ=255)                                                                                                gCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTgccag                                           <  1:432422/150‑1 (MQ=255)                                                                                                 cATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTgccagcc                                         >  2:244875/1‑151 (MQ=255)                                                                                                      tGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCCCATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCaa                                       <  1:518750/148‑1 (MQ=255)                                                                                                      tGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTgccagcca                                        <  1:353152/147‑1 (MQ=255)                                                                                                      tGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTgccagcca                                        <  1:244851/147‑1 (MQ=255)                                                                                                      tGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAAt                                      <  2:740958/149‑1 (MQ=255)                                                                                                        cgcgGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGt                                  <  2:51384/151‑1 (MQ=255)                                                                                                                acaacGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAg                          <  1:269185/151‑1 (MQ=255)                                                                                                                  aacGACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGcc                        >  2:246321/1‑151 (MQ=255)                                                                                                                     gACTGGCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCggg                     <  2:599498/151‑1 (MQ=255)                                                                                                                         ggCGGAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCaa                  <  2:156389/150‑1 (MQ=255)                                                                                                                             gAAGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTg             <  2:225731/151‑1 (MQ=255)                                                                                                                               aGTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGAt           >  1:696201/1‑151 (MQ=255)                                                                                                                                gTTTTTATTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGAt           >  1:358696/1‑150 (MQ=255)                                                                                                                                      aTTCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGa     >  2:101587/1‑150 (MQ=255)                                                                                                                                        tCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAAt   <  2:731764/150‑1 (MQ=255)                                                                                                                                        tCACCAGTTTGCGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATg  >  1:454372/1‑151 (MQ=255)                                                                                                                                        tCAACAATTTGCTCGCAAGCTGCTAGAAGACTTTAACATGATTCACTCCTTCGTTAACTCTCAGCGTCATCATGATCTCCAGTCCGACTTCTCATAATTCTGACTGCCAGGCAATCGGTTTCGCCAGCCGTGCAATAGTGATACTGGAAt   <  1:349238/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             GTCAGTTTGGTCTGACGCCGGGAGAGTACAGCGCCCGTTTTCAGGGCTAACGTCAGAAGGTTAATTCTGTTTCCAGCAGCGTCAGGATACTTTGTGCATCTTGCGCGGCAAACAACGACTGGCGGAAGTTTTTATTCACCAGTTTACGCGCCAGCTGCGAGAAGACTTTCACATGATTCATCCCTTCGTTAGCACCCAGCGTCAGCATGATCACCAGTTCGACTTCGCCCATTTCTGACTGCCAGCCAATCGGTTTCGCCAGCCGGGCAATGCTGATACTGGAATG  >  NC_000913/2501938‑2502223

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