BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,843,412	C→T	R12H (CGT→CAT)	hycH ←	formate hydrogenlyase assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,843,412	0	C	T	100.0%	199.5 / NA	50	R12H (CGT→CAT)	hycH	formate hydrogenlyase assembly protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (27/23); total (27/23)

CTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCG  >  NC_000913/2843366‑2843456

ctgctgCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATgg                                            >  1:2853930/1‑49 (MQ=255)

  gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                          >  1:2544266/1‑49 (MQ=255)

  gctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                          >  1:339679/1‑49 (MQ=255)

   ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                          >  1:2846970/1‑48 (MQ=255)

   ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                         <  1:2742134/49‑1 (MQ=255)

   ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTc                                         <  1:609245/49‑1 (MQ=255)

    tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                          >  1:966585/1‑47 (MQ=255)

    tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                          >  1:382517/1‑47 (MQ=255)

      cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                        <  1:1720435/47‑1 (MQ=255)

      cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGt                                      <  1:1659236/49‑1 (MQ=255)

       gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                     <  1:530059/49‑1 (MQ=255)

       gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                     <  1:1312480/49‑1 (MQ=255)

       gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                     <  1:650400/49‑1 (MQ=255)

       gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGCCTCTGtt                                     <  1:461710/49‑1 (MQ=255)

        ccTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                    <  1:2151055/49‑1 (MQ=255)

        ccTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCCCAGtt                                     >  1:2244920/1‑48 (MQ=255)

            ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                 >  1:112229/1‑48 (MQ=255)

            ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                <  1:2105937/49‑1 (MQ=255)

            ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                <  1:1502485/49‑1 (MQ=255)

            ggcgggcgTGGCATCGTTCTCATCAAAAAATTTATGGCTCAGTTGAc                                  <  1:512146/47‑1 (MQ=38)

             gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                >  1:2024176/1‑48 (MQ=255)

              cgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                              >  1:572167/1‑49 (MQ=255)

                ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                            <  1:242015/49‑1 (MQ=255)

                   gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                         >  1:502700/1‑49 (MQ=255)

                     ggCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                       >  1:417226/1‑49 (MQ=255)

                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                  <  1:151449/48‑1 (MQ=255)

                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                  <  1:1512993/48‑1 (MQ=255)

                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                  <  1:1613754/48‑1 (MQ=255)

                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                 >  1:1528523/1‑49 (MQ=255)

                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                 <  1:2842635/49‑1 (MQ=255)

                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                 >  1:328313/1‑49 (MQ=255)

                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                 >  1:380957/1‑49 (MQ=255)

                           gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                 >  1:387404/1‑49 (MQ=255)

                            ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                >  1:71742/1‑49 (MQ=255)

                             tctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt               >  1:1777487/1‑49 (MQ=255)

                              ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt               >  1:511256/1‑48 (MQ=255)

                               tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa             >  1:1465354/1‑49 (MQ=255)

                                catcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt            >  1:139518/1‑49 (MQ=255)

                                 atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc              <  1:1087506/46‑1 (MQ=255)

                                 atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATc           <  1:1256038/49‑1 (MQ=255)

                                   caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg          >  1:2221083/1‑48 (MQ=255)

                                    aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa        >  1:671281/1‑49 (MQ=255)

                                      tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa      >  1:2276724/1‑49 (MQ=255)

                                      tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa      >  1:2627912/1‑49 (MQ=255)

                                      tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa      >  1:1729689/1‑49 (MQ=255)

                                      tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa      <  1:288582/49‑1 (MQ=255)

                                       aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc     >  1:290057/1‑49 (MQ=255)

                                        aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc     <  1:2070026/48‑1 (MQ=255)

                                         aTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc   <  1:1944137/49‑1 (MQ=255)

                                          tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgcg  <  1:645076/49‑1 (MQ=255)

CTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCG  >  NC_000913/2843366‑2843456

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