powered by
writeHits(seqname, matches, strand, file = "", append = FALSE)
library(REDseq) x <- DNAString("AAGCGCGATATG") m <- matchPattern("GCGC", x) REDseq:::writeHits(seqname="chr1", m, strand="+",file="exampleWriteHits.bed", append=FALSE)
Run the code above in your browser using DataLab