I'm trying out DEXSeq for the identification of differentially expressed/spliced transcripts in my experiments. Using DEXSeq it is easy to identify the differentially used exons for each gene and by running
plotDEXSeq(dxr, "gene-name", displayTranscripts=T, names=T, legend=T)
it also possible to identify the transcripts for each gene.
Now, theoretically I would be able to identify which transcripts are present in which condition for each gene, based on these plots. However, this is extremely time consuming (even for all smaller number of genes) and I have to do this for multiple samples and hundreds or maybe thousands of genes. Therefore I was wondering if it is possible to produce a list of transcripts per gene for each condition using DEXSeq. Something like this:
|COND A||COND B|
which gives an easy overview of the transcripts expressed in each condition.
Thanks in advance,