1.将 gtf 文件下载好之后上传到RStudio-Server。
-
或者在自己的电脑上就直接读取就行。
gtf
如图所示我下载的是v43版本的基因组注释文件，现在最新版应该已经是v44了。
2.进行读取
经过搜索教程，发现读取gtf文件有一个专门的R包：rtracklayer
代码如下：
library(rtracklayer)
gtf <- rtracklayer::import('./gencode.v43.annotation.gtf')
gtf <- as.data.frame(gtf)
gtf$seqnames %>% table()
3.结果如下
gtf读取.png
