I wanted to use the 'circlize' package to make some genomic graphs. I wanted to follow the example given in the vignette/manual (https://jokergoo.github.io/circlize_book/book/initialize-genomic-plot.html - section 9.3)
tp_family = readRDS(system.file(package = "circlize", "extdata", "tp_family_df.rds")) head(tp_family) ## gene start end transcript exon ## 1 TP53 7565097 7565332 ENST00000413465.2 7 ## 2 TP53 7577499 7577608 ENST00000413465.2 6 circos.genomicInitialize(tp_family)
But this sets the beginning and end of the plot to the beginning and end of the exon coordinates. If I add a dummy exon to the data above, it will extend the plot, but will also show it as an exon.
- How can I modify the sample data to incude the 5' and 3' UTRs?
- Also, how can I then extend the range for the gene such that the plot shows +/- 2kbp from the 5' and 3' ends?
thanks for your help!