How to find the BSgenome object that provided the seqinfo of a GRanges object ?
1
0
Entering edit mode
@charles-plessy-7857
Last seen 14 months ago
Japan

Hello, I construct GRanges objects and pass them the seqinfo of a relevant BSgenome so that I get chromosome sizes. Sometimes, I would like to retrieve the DNA sequence of the ranges programmatically.

Practically it means that, using the information in the GRanges object I need to find the original BSgenome object. For instance, the seqinfo of BSgenome.Scerevisiae.UCSC.sacCer3, leads me to sacCer3 name, and so will do the GRanges objects that are contructed with the same seqinfo.

> BSgenome.Scerevisiae.UCSC.sacCer3 |> genome() |> unique()
[1] "sacCer3"

Is there a Bioc function that, when given "sacCer3" as input would return either the BSgenome object or its name ?

BSgenome • 1.0k views
ADD COMMENT
2
Entering edit mode
@herve-pages-1542
Last seen 1 day ago
Seattle, WA, United States

Yup, getBSgenome().

ADD COMMENT

Login before adding your answer.

Traffic: 548 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6