Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,534,031	T→C	intergenic (+10/‑35)	ptsH → / → ptsI	phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)/PEP‑protein phosphotransferase of PTS system (enzyme I)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,534,031	0	T	C	100.0%	163.3 / NA	46	intergenic (+10/‑35)	ptsH/ptsI	phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)/PEP‑protein phosphotransferase of PTS system (enzyme I)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/26);  total (20/26)

TTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCT  >  NC_000913/2533905‑2534161                                                                                                                               |                                                                                                                                   tttAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAAt                                                                                                                        <  2:213247/139‑1 (MQ=255)     aaCTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGt                                                                                                                    <  1:1169035/139‑1 (MQ=255)     aaCTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGt                                                                                                                    <  1:1015496/139‑1 (MQ=255)      aCTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTc                                                                                                                   <  1:1108976/139‑1 (MQ=255)           aGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAg                                                                                                              >  2:882950/1‑139 (MQ=255)           aGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAg                                                                                                              <  2:524378/139‑1 (MQ=255)               tctGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAg                                                                                                          <  1:882950/139‑1 (MQ=255)                ctGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAgg                                                                                                         <  2:520212/139‑1 (MQ=255)                       tGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGtt                                                                                                  >  1:1070944/1‑139 (MQ=255)                         aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGt                                                                                                                    >  2:621353/1‑119 (MQ=255)                         aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGt                                                                                                                    <  1:621353/119‑1 (MQ=255)                         aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTcac                                                                                                                 <  2:755905/122‑1 (MQ=255)                         aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTcac                                                                                                                 >  1:755905/1‑122 (MQ=255)                         aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAg                                                                                                              >  1:344524/1‑125 (MQ=255)                         aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAg                                                                                                              <  2:344524/125‑1 (MQ=255)                         aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTAt                                                                                                >  1:259791/1‑139 (MQ=255)                                gTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCa                                                                                         <  2:1070944/139‑1 (MQ=255)                                 tACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAg                                                                                        <  2:259791/139‑1 (MQ=255)                                   ccGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGc                                                                                      >  2:82934/1‑139 (MQ=255)                                           aCTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGc                                                                              >  1:1077699/1‑139 (MQ=255)                                                ctcCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCAtt                                                                                   <  1:30710/129‑1 (MQ=255)                                                ctcCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCAtt                                                                                   >  2:30710/1‑129 (MQ=255)                                                   cGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACaa                                                                                                               <  1:861086/98‑1 (MQ=255)                                                   cGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACaa                                                                                                               >  2:861086/1‑98 (MQ=255)                                                     cAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAgg                                                                                       <  2:908297/120‑1 (MQ=255)                                                     cAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAgg                                                                                       <  2:158062/120‑1 (MQ=255)                                                     cAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAgg                                                                                       >  1:158062/1‑120 (MQ=255)                                                     cAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAgg                                                                                       >  1:908297/1‑120 (MQ=255)                                                     cAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCATTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGt                                                                    >  2:197411/1‑139 (MQ=255)                                                      aGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCATTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTa                                                                   <  1:197411/139‑1 (MQ=255)                                                          ggCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAAt                                                                                                                        >  2:1288927/1‑82 (MQ=255)                                                          ggCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAAt                                                                                                                        <  1:1288927/82‑1 (MQ=255)                                                                  cGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTa                                                       <  1:82237/139‑1 (MQ=255)                                                                  cGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTa                                                       <  2:363230/139‑1 (MQ=255)                                                                      cAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATcc                                                                          <  1:1218945/116‑1 (MQ=255)                                                                      cAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATcc                                                                          >  2:1218945/1‑116 (MQ=255)                                                                         aaaGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATcc                                                                          <  2:61187/113‑1 (MQ=255)                                                                         aaaGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATcc                                                                          >  1:61187/1‑113 (MQ=255)                                                                              ggTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTc                                           >  2:828772/1‑139 (MQ=255)                                                                                          ggTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAAt                                                                                                                        >  1:688915/1‑50 (MQ=255)                                                                                          ggTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAAt                                                                                                                        <  2:688915/50‑1 (MQ=255)                                                                                             tAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCgg                                                                      >  2:807427/1‑97 (MQ=255)                                                                                             tAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCgg                                                                      <  1:807427/97‑1 (MQ=255)                                                                                                            aCTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGa             >  1:264852/1‑139 (MQ=255)                                                                                                                       ttCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCt  <  2:264852/139‑1 (MQ=255)                                                                                                                       ttCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCt  <  2:1077699/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                   TTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCT  >  NC_000913/2533905‑2534161

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