资讯

精准传达 • 有效沟通

从品牌网站建设到网络营销策划,从策略到执行的一站式服务

R语言如何实现可视化展示gff3格式基因组注释文件

这篇文章主要介绍R语言如何实现可视化展示gff3格式基因组注释文件,文中介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们一定要看完!

隆化网站制作公司哪家好,找成都创新互联!从网页设计、网站建设、微信开发、APP开发、响应式网站等网站项目制作,到程序开发,运营维护。成都创新互联2013年开创至今到现在10年的时间,我们拥有了丰富的建站经验和运维经验,来保证我们的工作的顺利进行。专注于网站建设就选成都创新互联


使用R语言的 ggbio这个包可视化gff3格式的基因组注释文件

我用到的文件是NCBI下载的拟南芥注释文件,为了减小计算压力,我只用到了gff文件的前119行,两个基因。

 首先是读入gff文件

用到的函数是 GenomicFeaturesR包中的 **makeTxDbFromGFF()**函数

library(GenomicFeatures)
txdb<-makeTxDbFromGFF(file="practice.gff",format="gff3")
   可视化

用到的 ggbio这个包中的 **autoplot()**这个函数

library(ggbio)
autoplot(txdb,
        which=GRanges("CP002684.1", IRanges(100, 9000)),
        names.expr = "gene_id")+
 theme_bw()
 

结果R语言如何实现可视化展示gff3格式基因组注释文件

可以通过fill参数设置不同的颜色

autoplot(txdb,
        which=GRanges("CP002684.1", IRanges(100, 9000)),
        names.expr = "gene_id",fill="red")+
 theme_bw()
 
R语言如何实现可视化展示gff3格式基因组注释文件  
image.png

不同的基因填充不同的颜色

autoplot(txdb,
        which=GRanges("CP002684.1", IRanges(100, 9000)),
        names.expr = "gene_id",aes(fill=gene_id))+
 theme_bw()
 
R语言如何实现可视化展示gff3格式基因组注释文件  
image.png

现在还不知道如何给同一个基因不同的部分(utr,exon,intron)等填充不同的颜色 还有就是 makeTxDbFromGFF()函数读入的数据存储格式还没搞懂

开头提到的参考资料里有一幅图将 reads数量, 覆盖度的折线图,vcf文件的结果,gff可视化的结果画到了一起,做基因组重测序分析应该会用得到。这里暂时不重复了。等用到的时候再说。R语言如何实现可视化展示gff3格式基因组注释文件

以上是“R语言如何实现可视化展示gff3格式基因组注释文件”这篇文章的所有内容,感谢各位的阅读!希望分享的内容对大家有帮助,更多相关知识,欢迎关注创新互联行业资讯频道!


分享名称:R语言如何实现可视化展示gff3格式基因组注释文件
当前链接:http://cdkjz.cn/article/pescgh.html
多年建站经验

多一份参考,总有益处

联系快上网,免费获得专属《策划方案》及报价

咨询相关问题或预约面谈,可以通过以下方式与我们联系

大客户专线   成都:13518219792   座机:028-86922220