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%	171.7 / NA	44	R12H (CGT→CAT)	hycH	formate hydrogenlyase assembly protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (19/25); total (19/25)

GCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACG  >  NC_000913/2843368‑2843454

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

 ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCt                                        <  1:1289761/48‑1 (MQ=255)

    cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                     <  1:900356/48‑1 (MQ=255)

    cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                     <  1:670895/48‑1 (MQ=255)

     gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                   >  1:2251179/1‑49 (MQ=255)

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

     gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                   >  1:6325/1‑49 (MQ=255)

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

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

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

          ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                <  1:754273/47‑1 (MQ=255)

           gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                             <  1:1816527/49‑1 (MQ=255)

           gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                             <  1:2086797/49‑1 (MQ=255)

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

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

               gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                         <  1:344531/49‑1 (MQ=255)

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

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

                    gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                     <  1:338567/48‑1 (MQ=255)

                    gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                     <  1:624459/48‑1 (MQ=255)

                    gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                    >  1:253965/1‑49 (MQ=255)

                     cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                   >  1:354443/1‑49 (MQ=255)

                       tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                 <  1:121076/49‑1 (MQ=255)

                        cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                 >  1:445005/1‑48 (MQ=255)

                        cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                 >  1:374057/1‑48 (MQ=255)

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

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

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

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

                            ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt             <  1:623072/48‑1 (MQ=255)

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

                               atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCa           <  1:14663/47‑1 (MQ=255)

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

                                tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCg        <  1:1795254/49‑1 (MQ=255)

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

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

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

                                  aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg       <  1:846975/48‑1 (MQ=255)

                                   aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa      >  1:1908894/1‑48 (MQ=255)

                                   aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt     <  1:692722/49‑1 (MQ=255)

                                   aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt     <  1:1125282/49‑1 (MQ=255)

                                      aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg  >  1:1860375/1‑49 (MQ=255)

                                      aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg  <  1:1360692/49‑1 (MQ=255)

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

GCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACG  >  NC_000913/2843368‑2843454

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