bambu icon indicating copy to clipboard operation
bambu copied to clipboard

Visualizing gene annotation on top of reported/discovered isoforms

Open imilenkovic opened this issue 3 years ago • 2 comments

Hello, I was asking if there is a way to add the gene annotation with all the annotated exons on top of the plot showing all isoforms (annotated and unannotated) that is produced using plotBambu? This way it would be easy to visually compare the detected isoforms to the annotated gene. Thanks!

imilenkovic avatar Feb 07 '22 14:02 imilenkovic

Hi @imilenkovic , unfortunately we don't provide a visualization for gene annotation on top of transcripts for now. We do provide visualization for gene level annotation and expression, but currently this functionality is a bit buggy. We are currently working to fix that bug and let you know once that is fixed maybe.

cying111 avatar Feb 08 '22 00:02 cying111

Hi @imilenkovic , we have fixed the bug for visualization of gene level annotation and expression. Below are the codes that you can follow: `seGene <- transcriptToGeneExpression(se) # se is the output from bambu

plotBambu(seGene, group.variable = NULL, type = 'annotation', gene_id = gene_id_of_interest, transcript_id = NULL)`

Thank you and let us know if you have further questions!

cying111 avatar Feb 14 '22 06:02 cying111