BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,697,216	T→C	intergenic (‑164/+57)	uidR ← / ← hdhA	DNA‑binding transcriptional repressor UidR/7‑alpha‑hydroxysteroid dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,697,216	0	T	C	100.0%	227.3 / NA	56	intergenic (‑164/+57)	uidR/hdhA	DNA‑binding transcriptional repressor UidR/7‑alpha‑hydroxysteroid dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (30/26); total (30/26)

GTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGG  >  NC_000913/1697170‑1697260

gTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCtga                                           >  1:2606882/1‑50 (MQ=255)

   aaCCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCtgatg                                         <  1:2738045/49‑1 (MQ=255)

    aCCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAg                                       <  1:1656904/50‑1 (MQ=255)

      cAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGtt                                     >  1:2798776/1‑50 (MQ=255)

      cAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGtt                                     <  1:2842113/50‑1 (MQ=255)

       aaTAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGttt                                    <  1:2205146/50‑1 (MQ=255)

       aaTAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGtt                                     >  1:1144397/1‑49 (MQ=255)

       aaTAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGtt                                     >  1:1227371/1‑49 (MQ=255)

        aTAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTa                                   <  1:439620/50‑1 (MQ=255)

         tAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTAt                                  <  1:2929618/50‑1 (MQ=255)

         tAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTAt                                  <  1:1805591/50‑1 (MQ=255)

         tAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTAt                                  >  1:572472/1‑50 (MQ=255)

         tAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTAt                                  >  1:847833/1‑50 (MQ=255)

           gCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATcc                                <  1:1965092/50‑1 (MQ=255)

           gCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATcc                                <  1:524303/50‑1 (MQ=255)

           gCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGCTTATcc                                >  1:1740211/1‑50 (MQ=39)

            ccGTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATccc                               <  1:3046682/50‑1 (MQ=255)

              gTGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCgg                             >  1:104049/1‑50 (MQ=255)

               tGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCgg                             >  1:1126628/1‑49 (MQ=255)

               tGATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCgg                             >  1:3432722/1‑49 (MQ=255)

               tGATGTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGt                            <  1:2628932/50‑1 (MQ=39)

                gATTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGt                            >  1:3389851/1‑49 (MQ=255)

                 aTTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTa                           >  1:2528260/1‑49 (MQ=255)

                 aTTTGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTa                           >  1:3343759/1‑49 (MQ=255)

                    tGGCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAaac                        <  1:1630244/49‑1 (MQ=255)

                      gCTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAaacaa                      <  1:2762760/49‑1 (MQ=255)

                       cTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAacaaca                    <  1:620327/50‑1 (MQ=255)

                       cTGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAacaaca                    >  1:1589999/1‑50 (MQ=255)

                        tGCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAacac                   >  1:3776078/1‑50 (MQ=255)

                         gCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacg                  >  1:643975/1‑50 (MQ=255)

                         gCCAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacg                  >  1:3337049/1‑50 (MQ=255)

                          ccAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacgc                 <  1:1452051/50‑1 (MQ=255)

                          ccAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacgc                 <  1:3158797/50‑1 (MQ=255)

                          ccAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacg                  <  1:3467243/49‑1 (MQ=255)

                           cAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacgca                <  1:1930644/50‑1 (MQ=255)

                           cAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacgca                >  1:3867564/1‑50 (MQ=255)

                           cAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacgca                >  1:2766873/1‑50 (MQ=255)

                           cAGTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacgca                <  1:1547684/50‑1 (MQ=255)

                             gTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacgcacg              <  1:1276640/50‑1 (MQ=255)

                             gTTATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAAcacgcacg              <  1:371556/50‑1 (MQ=255)

                              ttATTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACgg             >  1:3297029/1‑50 (MQ=255)

                                aTTTAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgtt           >  1:402295/1‑50 (MQ=255)

                                 tttAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgttg          <  1:2927720/50‑1 (MQ=255)

                                 tttAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgtt           <  1:1304834/49‑1 (MQ=255)

                                  ttAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgttgt         >  1:2310878/1‑50 (MQ=255)

                                   tAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgttgtt        >  1:2946348/1‑50 (MQ=255)

                                   tAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgttgtt        >  1:2088531/1‑50 (MQ=255)

                                   tAGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgttgtt        >  1:1273187/1‑50 (MQ=255)

                                    aGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgttgtt        >  1:2472909/1‑49 (MQ=255)

                                    aGCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGgttgtt        >  1:359974/1‑49 (MQ=255)

                                     gCAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGttt       <  1:848410/49‑1 (MQ=255)

                                      cAGAGACGCTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTAc     <  1:2386329/50‑1 (MQ=255)

                                       agagACGCTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTAcc    >  1:3056509/1‑50 (MQ=255)

                                       agagACGCTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTAcc    >  1:2431103/1‑50 (MQ=255)

                                        gagaCGCTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCg   <  1:487335/50‑1 (MQ=255)

                                         agaCGCTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCgg  >  1:729010/1‑50 (MQ=255)

GTTAACCAATAGCCGTGATTTGGCTGCCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGG  >  NC_000913/1697170‑1697260

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