Hi, To begin I would like to thanks this community for all the work done and the spirit of sharing and helping each other it is very pleasant and rare nowadays. I am a french internship student in pathology and I am doing a work using restriction enzymes. I would like to see the distribution of the RE sites along the different chromosomes in order to see if there is some regions undercovered by one enzyme and which combinaison of two enzymes would fix this. In order to do that I found the REDSeq package, done by Lihua Julie Zhu. I worked on it yesterday but i'm facing a problem that i can't solve by myself.
My command is " library(REDseq) REpatternFilePath = myFile <- file.choose() library(BSgenome.Hsapiens.UCSC.hg38) buildREmap( REpatternFilePath, BSgenomeName= BSgenome.Hsapiens.UCSC.hg38, outfile= tempfile()) " My question is how to use the REmap created to export the data in a .csv file for example ? I don't know how to manage the outfile function.
Thanks a lot for any kind of help.