buildREmap {REDseq} | R Documentation |
Build a genome-wide cut map for a Restriction Enzyme (RE)
buildREmap(REpatternFilePath, format = "fasta", BSgenomeName, outfile)
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 |
Output REmap as a RangedData
Lihua Julie Zhu
library(REDseq) REpatternFilePath = system.file("extdata", "examplePattern.fa", package="REDseq") library(BSgenome.Celegans.UCSC.ce2) buildREmap( REpatternFilePath, BSgenomeName=Celegans, outfile=tempfile())