Hi, I wondering if anybody can help me to generate masked (by RepeatMasker for instance) sequences.
I'm currently using Bsgenome to extract sequence from a BED file such as : library(BSgenome.Hsapiens.UCSC.hg18) genome<-Hsapiens FastaSeq<-getSeq(genome,"chr1",start=1000,end=1200, as.character=FALSE) I know that Bsgenome contains masks that can be apply by using : chr1 <- genome$chr1 active(masks(chr1)) <- TRUE So, I'm trying to use it to change the masks of the genome object. But I cannot modify it : active(masks(genome$chr1)) <- TRUE Error in `$<-`(`*tmp*`, "chr1", value = <S4 object of class "MaskedDNAString">) : no method for assigning subsets of this S4 class Is there a way get the masked sequence with the getSeq function ? Thanks. Arnaud. _______________________________________________ Bioc-sig-sequencing mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
