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

Only 100 of 132 total aligned reads displayed.

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458

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

    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  1:4710015/50‑1 (MQ=255)

    ctgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCa                                          <  2:2479460/50‑1 (MQ=255)

     tgcGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAg                                         >  1:3766518/1‑50 (MQ=255)

       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  2:938626/50‑1 (MQ=255)

       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:5140394/50‑1 (MQ=255)

       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  1:5269309/50‑1 (MQ=255)

       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       <  2:3684317/50‑1 (MQ=255)

       cgcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGtt                                       >  2:4012362/1‑50 (MQ=255)

        gcCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTg                                      <  1:4328564/50‑1 (MQ=255)

         ccTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGa                                     <  2:4127530/50‑1 (MQ=255)

          cTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGAc                                    >  2:3189517/1‑50 (MQ=255)

          cTCGGCGGGCGTGGCATCGGTCTCCTCAAAAAATTTATGGCGCCGTTGAc                                    >  2:3887283/1‑50 (MQ=255)

            cggcggGCGTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   >  1:2781519/1‑49 (MQ=255)

             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:4244625/48‑1 (MQ=255)

             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  2:3164343/48‑1 (MQ=255)

             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  2:4573134/48‑1 (MQ=255)

             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACt                                   <  1:2033988/48‑1 (MQ=255)

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

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

             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTg                                  >  1:89881/1‑49 (MQ=255)

             ggcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGa                                 <  2:2381941/50‑1 (MQ=255)

              gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                >  2:1526646/1‑50 (MQ=255)

              gcgggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGaa                                <  2:4555944/50‑1 (MQ=255)

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

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

                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  2:5146460/50‑1 (MQ=255)

                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  2:5685992/50‑1 (MQ=255)

                gggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAaca                              <  1:4031275/50‑1 (MQ=255)

                 ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  2:727991/50‑1 (MQ=255)

                 ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  1:1079957/50‑1 (MQ=255)

                 ggcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             <  1:2299458/50‑1 (MQ=255)

                  gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  1:1032287/1‑49 (MQ=255)

                  gcgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAacac                             >  1:2370552/1‑49 (MQ=255)

                   cgTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcacca                           <  1:4272397/50‑1 (MQ=255)

                    gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          >  2:172316/1‑50 (MQ=255)

                    gTGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAAcaccac                          <  2:1447572/50‑1 (MQ=255)

                     tGGCATCGTTCTCATCACTAAATTTATGGCTCAGTTGCCTGAACaccacc                         <  2:430948/50‑1 (MQ=255)

                     tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         <  1:499620/50‑1 (MQ=255)

                     tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         <  1:2576119/50‑1 (MQ=255)

                     tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         <  1:5443907/50‑1 (MQ=255)

                     tGGCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACaccacc                         <  2:1119532/50‑1 (MQ=255)

                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       <  2:4933687/50‑1 (MQ=255)

                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       <  2:368632/50‑1 (MQ=255)

                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtt                       >  1:5060659/1‑50 (MQ=255)

                       gCATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCt                        <  1:1778689/49‑1 (MQ=255)

                        cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      >  1:4351031/1‑50 (MQ=255)

                        cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  2:3105133/50‑1 (MQ=255)

                        cATCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCttt                      <  1:5305103/50‑1 (MQ=255)

                         aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:1549569/1‑50 (MQ=255)

                         aTCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     <  1:2605353/50‑1 (MQ=255)

                          tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCtttt                     >  1:2416092/1‑49 (MQ=255)

                          tCGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTc                    <  2:4032357/50‑1 (MQ=255)

                           cGTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCa                   >  2:1776568/1‑50 (MQ=255)

                            gTTCTCATCAATAAATTTATTGCTCAGTTGACTGAACACCACCTTTTCAc                  <  1:1918184/50‑1 (MQ=39)

                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCTCCTTTTCAc                  <  1:2963766/50‑1 (MQ=255)

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

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

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

                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:1830386/1‑50 (MQ=255)

                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  <  1:3509600/50‑1 (MQ=255)

                            gTTCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCAc                  >  1:4655607/1‑50 (MQ=255)

                             ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  1:342116/1‑50 (MQ=255)

                             ttCTCATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACt                 >  2:2431135/1‑50 (MQ=255)

                               ctcATCAATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTc                <  2:1419224/49‑1 (MQ=255)

                                tcatcaATAAATTTATGTCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  2:4704845/1‑50 (MQ=255)

                                tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  2:4253469/1‑50 (MQ=255)

                                tcatcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              >  1:759919/1‑50 (MQ=255)

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

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

                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:5020299/48‑1 (MQ=255)

                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  2:931524/48‑1 (MQ=255)

                                  atcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCAt              <  1:2097403/48‑1 (MQ=255)

                                   tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg           <  1:255020/50‑1 (MQ=255)

                                   tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg           <  1:745057/50‑1 (MQ=255)

                                   tcaATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCgg           <  1:1239676/50‑1 (MQ=255)

                                    caATAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          >  2:5068366/1‑50 (MQ=255)

                                     aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGa          <  2:380357/49‑1 (MQ=255)

                                     aaTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGAt         <  2:1353933/50‑1 (MQ=255)

                                      aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        <  2:508980/50‑1 (MQ=255)

                                      aTAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATa        <  2:216530/50‑1 (MQ=255)

                                       tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:2171696/50‑1 (MQ=255)

                                       tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  1:3993939/50‑1 (MQ=255)

                                       tAAATTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAc       <  2:576897/50‑1 (MQ=255)

                                        aaaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcg      >  2:4770665/1‑50 (MQ=255)

                                        aaaTTTATGGCTCAGATGACTGAACACCACCTTTTCACTCATCGGATAcg      >  2:5503135/1‑50 (MQ=255)

                                         aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     >  1:3820316/1‑50 (MQ=255)

                                         aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  1:1358696/50‑1 (MQ=255)

                                         aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     <  2:2817339/50‑1 (MQ=255)

                                         aaTTTATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATAcgc     >  2:3511215/1‑50 (MQ=255)

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

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

                                         aaTTTATGGCTCAGTTGACTGAACACCAACTTTTCACTCATCGGATAcgc     <  1:1001538/50‑1 (MQ=39)

                                          aTTTATGGCTCAGTTGACTGAACACCACCCTTTCACTCATCGGATAcgcg    >  2:816725/1‑50 (MQ=255)

                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   <  2:2677045/50‑1 (MQ=255)

                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   >  2:5432024/1‑50 (MQ=255)

                                           tttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACgcgc   <  2:366877/50‑1 (MQ=255)

                                            ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  1:2972610/1‑50 (MQ=255)

                                            ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  2:4818608/1‑50 (MQ=255)

                                            ttATGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGcc  >  2:5237172/1‑50 (MQ=255)

CCTGCTGCGCCTCGGCGGGCGTGGCATCGTTCTCATCAATAAATTTACGGCTCAGTTGACTGAACACCACCTTTTCACTCATCGGATACGCGCC  >  NC_000913/2843365‑2843458

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