On 06/01/2011 09:06 AM, Ivan Gregoretti wrote:
Hello IRanges connoisseurs,

Is there a function like subseq() but with the ability to accept a
vector of starts?


For instance, if I want to chop the beginning of the murine
mitochondrial DNA into 50 nucleotides, I'd do

library(BSgenome.Mmusculus.UCSC.mm9)

subseq(Mmusculus[["chrM"]], start=1,   width=50)
subseq(Mmusculus[["chrM"]], start=51,  width=50)
subseq(Mmusculus[["chrM"]], start=101, width=50)

...

however, it would be more convenient something like this

subseq(Mmusculus[["chrM"]], start=c(1,51,101), width=50)

Thank you

Ivan

_______________________________________________
Bioc-sig-sequencing mailing list
Bioc-sig-sequencing@r-project.org
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing


--
Computational Biology
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N. PO Box 19024 Seattle, WA 98109

Location: M1-B861
Telephone: 206 667-2793

_______________________________________________
Bioc-sig-sequencing mailing list
Bioc-sig-sequencing@r-project.org
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing

Reply via email to