Relative position of the start and stop codon along the transcript

Share:

Description

Relative position of the start and stop codon along the transcript

Usage

1
orfRelativePos(cdsTransc, exonGRanges)

Arguments

cdsTransc

a GRangesList. It contains the CDS coordinates grouped by transcript.

exonGRanges

a GRangesList. It contains the exon coordinates grouped by transcript.

Value

a list. A list of relative positions of the start and end of ORFs.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
#make a txdb object containing the annotations for the specified species.
#In this case hg19.
txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene::TxDb.Hsapiens.UCSC.hg19.knownGene

#get all CDSs by transcript
cds <- GenomicFeatures::cdsBy(txdb, by="tx", use.names=TRUE)

#get all exons by transcript
exonGRanges <- GenomicFeatures::exonsBy(txdb, by="tx", use.names=TRUE)

#retrieve the positions of start and end codons relative to the transcript
cdsPosTransc <- orfRelativePos(cds, exonGRanges)