已知转录本test.ID，根据biomart包批量获得ensembl数据库中对应序列信息，并以fa格式输出

该方法需要魔法，我已进行了调试，可以成功运行

转录本test.ID样式(ENSMUST00000040608.4)，以xlsx形式储存

library(Biostrings)
library(biomaRt)
library(openxlsx)
library(seqinr)rm(list = ls())
genelist <- read.xlsx("待测基因testID汇总.xlsx")
genelist <- genelist$test_id
#移除genelist的小数点
genelist<-unlist(lapply(genelist, FUN = function(x) {return(strsplit(x, split = ".",fixed = T)[[1]][1])}))mart <- useMart("ensembl", dataset="mmusculus_gene_ensembl") # 小鼠的数据库
Ensemble_to_seq <- function(x) {cat("Getting sequence for gene", x, "\n")  # 添加打印语句seq = biomaRt::getSequence(id = x, type = "ensembl_transcript_id", seqType = "cdna", mart = mart)seq<-as.data.frame(seq)seq$"Length"<-lapply(seq[,1],function(y){return(nchar(y))})return(seq)
}outTab <- data.frame()
for (i in genelist){outTab<-rbind(outTab,Ensemble_to_seq(i))} # 此步需要魔法# 保存为fa格式文件
for (i in 1:nrow(outTab)) {# 定义序列seq <- DNAStringSet(outTab$cdna[i])# 定义描述信息descr <- paste(outTab$ensembl_transcript_id[i], "length=", outTab1$Length[i], sep = "|")# 输出序列write.fasta(seq, names = outTab$ensembl_transcript_id[i], file.out = "output.fa", open = "a")
}

最后是输出结果展示

已知转录本test.ID，根据biomart包批量获得ensembl数据库中对应序列信息，并以fa格式输出相关推荐

1.已知本原多项式，利用Matlb中的simulink构成m序列产生器。2.已知任意本原多项式，利用matlb软件编程求解其对应的m序列以及m序列产生过程。
1. 已知本原多项式,利用Matlb中的simulink构成m序列产生器.2.已知任意本原多项式,利用matlb软件编程求解其对应的m序列以及m序列产生过程. m序列是最长线性反馈移位寄存器的简称,他 ...
已知三角形三边长怎么求面积_解三角形问题中的常见错解分析
解三角形问题是个难点,怎样才能突破这个难点呢? 只有正确理解三角形中的边角关系,即三角形中的边角等量关系.边角的不等关系及内角和关系,才能克服这个难点. 下面快和包sir一起对解三角形问题中的常见错误 ...
mysql查询已知两列第三列,MySQL查询从两列中查找出现的次数？
使用MySQL GROUP_BY从两列中查找出现的次数.让我们首先创建一个表-mysql> create table DemoTable -> ( -> Name1 varchar( ...
c语言已知一个班有10位同学,一个班10位同学参加期末c语言考试,先要输出全班同学的成绩,统计全班同学的平...
#include #define n 5 struct student { int num; char name[9]; float math; float english; float physic ...
爬取知网博硕士文献及中国专利存到mysql数据库中的代码及其注意事项
今天因为需要做了一个爬取知网博硕士论文及中国专利的爬虫,在制作的过程中遇到了不少坑,在网上查资料时都是很老的资源,在现在知网的反爬虫下不起作用,所以我来写这篇文章来供大家参考.(这篇文章主要介绍通过改 ...
ensembl数据库BioMart使用，查找转录本ID在ensembl与refseq中的对应关系
如何查找转录本ID在ensembl与refseq中的对应关系 Ensembl数据库地址: http://grch37.ensembl.org/index.html 进入BioMart 选择DataS ...
jQuery九宫格跑马灯抽奖（已知结果）
描述: 写两个抽奖,等级不同,找到一个原生写法的,原作者链接如下: 参考( https://www.cnblogs.com/yangqing22/p/14065754.html ) 略改了一下,换成了 ...
已知空间一点到另外两点直线的距离
转自:http://www.cnblogs.com/clarkustb/archive/2008/11/04/1326500.html 已知空间一点到另外两点直线的距离设空间中的三点为M1,M2,M ...
Win7安装VC++6.0已知的兼容性问题的解决方法
转自:http://www.veryhuo.com/a/view/53057.html 今天写程序让人生气的是写的所有的程序没一个对的,不能连接运行.原来经过好多次的纠结后来才知道是win7不兼容导致 ...

已知转录本test.ID，根据biomart包批量获得ensembl数据库中对应序列信息，并以fa格式输出

已知转录本test.ID，根据biomart包批量获得ensembl数据库中对应序列信息，并以fa格式输出相关推荐

最新文章

热门文章