#从基因组的gff3文件创建txdb数据库
#BiocManager::install("GenomicFeatures")
#BiocManager::install("AnnotationHub")
library(GenomicFeatures)
library(AnnotationHub)
setwd("E:/bioinformation/")
gff3file <- "ZmB73V5.genome.gff3" #本地需要存在这个基因组的gff3文件
#从gff3文件创制txdb数据库
txdbV5 <- makeTxDbFromGFF(file=gff3file,
dataSource="Zeamays (maize) B73_V5",
organism="Zea mays",
taxonomyId = 4577) #这里是物种在NCBI分类数据库的编号
#保存txdb到本地
saveDb(txdbV5,"Zeamay_B73_V5.txdb")
#从本地加载txdb
txdb <- loadDb("./Zeamay_B73_V4.txdb")#加载V4版本的txdb
txdb <- loadDb("./Zeamay_B73_V5.txdb")#加载V5版本的txdb
#查看数据库的内容
genes(txdb)$gene_id
exons(txdb)
从基因组的gff3文件创建txdb数据库
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 截止2021年12月22日,从NCBI下载GRCh38注释文件gff3的方法有:方法一:访问网址 https://...
- 本文参考自Repeat Library Construction-Advanced,整体思路一致,但是所用软件有所...
- def reverse_complement(seq): ntComplement = {'A':'T','T'...
- GDSC(Genomics of Drug Sensitivity in Cancer)数据库由英国桑格研究院开发...
- 研究开始前搜集、整理已有研究文献是常规操作,这个步骤往往需要耗费大量时间和精力。小编今天介绍的这个数据库,可以让你...