I am trying to get the concatenated GRanges of ORFs for each gene.
I thought I could just get a list of cds from the txdb database and concatenate this list from the first to the last position of each cds. But I got stuck
What i do so far is:
cds <- cds(txdb, columns=c("TXNAME","EXONRANK"))
cds_grl <- multisplit(cds, cds$TXNAME)
wich returns a list of GRange object
Do you know how I could concatenate each GR element of the list to get the first and last possible CDS GRranges for each gene?
Maybe there is another method ?