Utility functions related to RNA/DNA sequences, such as extracting
RNA/DNA sequences for features defined in
1 2 3 4 5
(In alphabetic order)
object with the genomic DNA sequence.
object with the genome sequence.
FaFile-class (defined in
Rsamtools) with the genomic sequence of the genome build
matching the Ensembl version of the
The file is retrieved using the
thus, at least for the first invocation, an internet connection is
required to locate and download the file; subsequent calls will
load the cached file instead.
If no fasta file for the actual Ensembl version is available the
function tries to identify a file matching the species and genome
build version of the closest Ensembl release and returns that
See the vignette for an example to work with such files.
TwoBitFile-class (defined in the
rtracklayer package) with the genomeic sequence of the
genome build matching the Ensembl version of the
object. The file is retrieved from
AnnotationHub and hence
requires (at least for the first query) an active internet
connection to download the respective resource. If no DNA sequence
matching the Ensembl version of
x is available, the
function tries to find the genomic sequence of the best matching
genome build (closest Ensembl release) and returns that.
ensembldb vignette for details.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
## Loading an EnsDb for Ensembl version 86 (genome GRCh38): library(EnsDb.Hsapiens.v86) edb <- EnsDb.Hsapiens.v86 ## Not run: ## Retrieve a TwoBitFile with the gneomic DNA sequence matching the organism, ## genome release version and, if possible, the Ensembl version of the ## EnsDb object. Dna <- getGenomeTwoBitFile(edb) ## Extract the transcript sequence for all transcripts encoded on chromosome ## Y. ##extractTranscriptSeqs(Dna, edb, filter=SeqNameFilter("Y")) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.