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%	168.0 / 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 (24/22);  total (24/22)

CTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCTGCCGA  >  NC_000913/2533902‑2534166                                                                                                                                  |                                                                                                                                        cTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTaaa                                                                                                                                <  2:481227/139‑1 (MQ=255)              aGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTc                                                                                                                                       >  2:390226/1‑119 (MQ=255)                     gggCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTa                                                                                                            >  1:507162/1‑139 (MQ=255)                            aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCGGTc                                                                                                                        >  1:950268/1‑120 (MQ=255)                            aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCGGTc                                                                                                                        <  2:950268/120‑1 (MQ=255)                            aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGt                                                                                                                         >  2:550418/1‑119 (MQ=255)                            aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGt                                                                                                                         <  1:550418/119‑1 (MQ=255)                            aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAgg                                                                                                              <  2:801411/130‑1 (MQ=255)                            aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAgg                                                                                                              >  1:801411/1‑130 (MQ=255)                            aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTAt                                                                                                     >  2:813245/1‑139 (MQ=255)                                aaGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAg                                                                                                               >  1:324159/1‑125 (MQ=255)                                aaGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAg                                                                                                               <  2:324159/125‑1 (MQ=255)                                             gACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAg                                                                                    <  1:813245/139‑1 (MQ=255)                                              aCTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAgg                                                                                            >  1:483926/1‑130 (MQ=255)                                              aCTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAgg                                                                                            <  2:483926/130‑1 (MQ=255)                                              aCTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGc                                                                                   >  1:858455/1‑139 (MQ=255)                                                   ctcCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAAt                                                                                                                             >  2:730815/1‑92 (MQ=255)                                                   ctcCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAAt                                                                                                                             <  1:730815/92‑1 (MQ=255)                                                   ctcCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATccc                                                                              >  2:583302/1‑139 (MQ=255)                                                      cGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTa                                                                                                            >  2:867718/1‑106 (MQ=255)                                                      cGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTa                                                                                                            <  1:867718/106‑1 (MQ=255)                                                             ggCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGc                                                                    >  1:101263/1‑139 (MQ=255)                                                                gAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCttt                                                                 <  2:858455/139‑1 (MQ=255)                                                                    aCGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGt                                                             <  1:267185/139‑1 (MQ=255)                                                                       aGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCTAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTaaa                                                          >  1:714394/1‑139 (MQ=255)                                                                                  gTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCt                                               >  2:238784/1‑139 (MQ=255)                                                                                      aaCATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTa                                                                                                            <  2:499573/74‑1 (MQ=255)                                                                                      aaCATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTa                                                                                                            >  1:499573/1‑74 (MQ=255)                                                                                             ggTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCt                                                                   <  2:215982/108‑1 (MQ=255)                                                                                             ggTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCt                                                                   >  1:215982/1‑108 (MQ=255)                                                                                                       aTGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAg                                                         <  1:157912/108‑1 (MQ=255)                                                                                                       aTGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAg                                                         >  2:157912/1‑108 (MQ=255)                                                                                                        tGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAgg                                                                                            <  1:100266/72‑1 (MQ=255)                                                                                                        tGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAgg                                                                                            >  2:100266/1‑72 (MQ=255)                                                                                                        tGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCAtt                         >  2:901239/1‑139 (MQ=255)                                                                                                             gAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGt                                                                                                                         >  1:898792/1‑38 (MQ=255)                                                                                                             gAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGt                                                                                                                         <  2:898792/38‑1 (MQ=255)                                                                                                                    aGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGct                                                       >  1:300706/1‑97 (MQ=255)                                                                                                                    aGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGct                                                       <  2:300706/97‑1 (MQ=255)                                                                                                                    aGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGaaaaaa             >  2:137206/1‑139 (MQ=255)                                                                                                                           tCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGTGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCTg      <  1:583302/139‑1 (MQ=255)                                                                                                                             ccGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCa                                                                                  >  1:171395/1‑61 (MQ=255)                                                                                                                             ccGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCa                                                                                  <  2:171395/61‑1 (MQ=255)                                                                                                                               gggCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCTGCCGa  <  2:714394/139‑1 (MQ=255)                                                                                                                               gggCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCTGCCGa  <  1:238784/139‑1 (MQ=255)                                                                                                                               gggCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCTGCCGa  <  2:101263/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        CTGTTTAAACTGCAGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGACGAAATTGTCATTGACCGGAAAAAAATTTCTGCCGA  >  NC_000913/2533902‑2534166

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