Constructs CDS information on transcripts from custom annotation using reference-based approach.

buildCDS(object, verbose = FALSE)

Arguments

object

factRObject

verbose

Whether to print out messages (Default: FALSE)

Value

Updated factRObject. The GTF GenomicRanges object will now contain "CDS" type features.

The `cds` variable from the transcript metadata will return "yes" for coding transcripts and "no" for non-coding transcripts

Details

This function will firstly identify custom transcripts with identical exon structures as those in the reference annotation. If these reference transripts contain CDS segments, the coordinates will be passed to its paired- custom transcript. The putative start codon of the remaining transcripts will be the first ATG sequence that is in-frame with the coding sequence of reference transcripts and an in-frame stop codon will be determined.

See also

Examples

## Load sample data
data(factRsample)

## Build coding sequences
factRsample <- buildCDS(factRsample)