简单画个基因结构吧

引言
有时候只想看看基因的结构，没有其他信号数据输入什么的。简化版的 transcript_plot 来干这个事情。参数大部分和 coverage_plot 一致。
devtools::install_github(“junjunlab/omicScope”)
library(omicScope)

gtf <- rtracklayer::import(“../../index-data/Mus_musculus.GRCm38.102.gtf.gz”)

transcript_plot(gtf_file = gtf,
target_gene = c(“Actb”,”Gapdh”))

添加 transcript_id:
transcript_plot(gtf_file = gtf,
target_gene = c(“Actb”,”Gapdh”),
add_gene_label = TRUE,
gene_label_aes = “transcript_id”,
gene_label_vjust = 0.25)

去掉矩形背景：
transcript_plot(gtf_file = gtf,
target_gene = c(“Actb”,”Gapdh”),
add_gene_label = TRUE,
gene_label_aes = “transcript_id”,
gene_label_vjust = 0.25,
add_backsqure = FALSE)

指定基因组区间绘图：
transcript_plot(gtf_file = gtf,
target_region = c(“1:127061844-128745904”,
“2:121127504-121354559”),
add_gene_label = TRUE,
gene_label_aes = “gene_name”,
gene_label_vjust = 0.25,
add_backsqure = FALSE)

结尾
路漫漫其修远兮,吾将上下而求索。
欢迎加入生信交流群。加我微信我也拉你进 微信群聊老俊俊生信交流群(微信交流群需收取 20 元入群费用,一旦交费,拒不退还!(防止骗子和便于管理)) 。