Extract the DNA sequences of all exons of transcript models and concatenate to one sequence per
transcript. This is basically a wrapper for
which makes sure that the exons are correctly sorted according to their position in the
transcript (3' to 5').
1 2 3 4 5 6 7
DNAStringSet object containing the transcript
GRangesList: Obtain transcript sequence from
GRanges: Obtain transcript sequence from
1 2 3 4 5 6 7 8 9 10 11
library(BSgenome) # protein-coding exons of transcripts within chr11 region data("chr11_genes") target_txs <- split(chr11_genes, f = chr11_genes$transcript_id) # human genome (hg38) BSgenome object (needs to be installed separately from Bioconductor) hg38 <- getBSgenome("BSgenome.Hsapiens.UCSC.hg38") # get sequences for all target transcripts on chr11 txs_seqs <- getTxsSeq(target_txs, genome = hg38)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.