get_enzyme_cutsites: get_enzyme_cutsites

In mervesa/HiCDCPlus: Hi-C Direct Caller Plus

get_enzyme_cutsites

Description

This function finds all restriction enzyme cutsites of a given genome, genome version, and set of cut patterns

Usage

get_enzyme_cutsites(sig, gen = "Hsapiens", gen_ver = "hg19", chrs = NULL)

Arguments

sig	a set of restriction enzyme cut patterns (e.g., 'GATC' or c('GATC','GANTC'))
gen	name of the species: e.g., default 'Hsapiens'
gen_ver	genomic assembly version: e.g., default 'hg19'
chrs	a subset of chromosomes' e.g., c('chr21','chr22'). Defaults to all chromosomes (except Y and M) in the genome specified by gen and gen_ver.

Value

list of chromosomes.

Examples

get_enzyme_cutsites(gen='Hsapiens',gen_ver='hg19',
sig=c('GATC','GANTC'),chrs=c('chr22'))

mervesa/HiCDCPlus documentation built on June 8, 2022, 3:43 a.m.