BRESEQ :: Evidence

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%	98.8 / NA	28	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 (14/14); total (14/14)

AGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGA  >  NC_000913/2533915‑2534130

agcctctGTGACTGACTCAAGGTACCGCCGGAACTTTCGCCACAGAAGATACAGACGAGCAGAAAACGGCTGCACAGCTGGTTATACTCATCGCACAACTCGAACAAATTCCCGGGCGCTTTTAAAAGTCAGTCACAAg                                                                               >  2:672273/4‑139 (MQ=255)

aGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAg                                                                               >  2:1063502/1‑139 (MQ=255)

  aCTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTcaca                                                                                 >  2:261536/1‑135 (MQ=255)

  aCTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTcaca                                                                                 <  1:261536/135‑1 (MQ=255)

     ctGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAgg                                                                          <  1:1066558/139‑1 (MQ=255)

         gCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAg                                                                               <  1:187645/130‑1 (MQ=255)

         gCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAg                                                                               >  2:187645/1‑130 (MQ=255)

              aCTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTAt                                                                 >  2:800535/1‑139 (MQ=255)

                     gTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCa                                                          <  1:1063502/139‑1 (MQ=255)

                         cGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTc                                                                                    <  2:920386/109‑1 (MQ=255)

                         cGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTc                                                                                    >  1:920386/1‑109 (MQ=255)

                          gTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTc                                                                                    >  2:672724/1‑108 (MQ=255)

                          gTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTc                                                                                    <  1:672724/108‑1 (MQ=255)

                               gACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAg                                                <  2:82615/139‑1 (MQ=255)

                                aCTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAgg                                                                          >  2:313060/1‑112 (MQ=255)

                                aCTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAgg                                                                          <  1:313060/112‑1 (MQ=255)

                                   aTCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATc                                            >  2:161941/1‑139 (MQ=255)

                                    tctcCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATcc                                           >  2:194046/1‑139 (MQ=255)

                                               ggCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCa                                                          >  1:611982/1‑113 (MQ=255)

                                               ggCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCa                                                          <  2:611982/113‑1 (MQ=255)

                                               ggCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGc                                >  2:1290088/1‑139 (MQ=255)

                                                gCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCt                               <  1:1091521/139‑1 (MQ=255)

                                                gCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCt                               >  1:552116/1‑139 (MQ=255)

                                                     gACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCgg                          <  1:800535/139‑1 (MQ=255)

                                                            aGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGct                   <  1:1290088/139‑1 (MQ=255)

                                                                              tGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAagaaga  <  2:63530/138‑1 (MQ=255)

                                                                              tGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAagaaga  >  1:63530/1‑138 (MQ=255)

                                                                                                            ttCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGc               <  2:376349/95‑1 (MQ=255)

                                                                                                            ttCCCGGGCTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGc               >  1:376349/1‑95 (MQ=255)

AGACTCTGGGCCTGACTCAAGGTACCGTTGTGACTATCTCCGCAGAAGGCGAAGACGAGCAGAAAGCGGTTGAACATCTGGTTAAACTGATGGCGGAACTCGAGTAATTTCCCGGGTTCTTTTAAAAATCAGTCACAAGTAAGGTAGGGTTATGATTTCAGGCATTTTAGCATCCCCGGGTATCGCTTTCGGTAAAGCTCTGCTTCTGAAAGAAGA  >  NC_000913/2533915‑2534130

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