REDseq (version 1.18.0)

buildREmap: Build a genome wide cut site map for a Restriction Enzyme (RE)

Description

Build a genome-wide cut map for a Restriction Enzyme (RE)

Usage

buildREmap(REpatternFilePath, format = "fasta", BSgenomeName, outfile)

Arguments

REpatternFilePath

File path storing the recognition pattern of a RE

format

format of the pattern file, either "fasta" (the default) or "fastq

BSgenomeName

BSgenome object, please refer to available.genomes in BSgenome package for details

outfile

temporary output file for writing the matched chromosome location to

Value

Output REmap as a RangedData

Examples

	library(REDseq)
	REpatternFilePath = system.file("extdata", "examplePattern.fa", package="REDseq")
	library(BSgenome.Celegans.UCSC.ce2)
	buildREmap( REpatternFilePath, BSgenomeName=Celegans, outfile=tempfile())