首页 > 其他 > 详细

序列处理之seqinr(fetch)

时间:2018-09-17 12:13:49      阅读:470      评论:0      收藏:0      [点我收藏+]

缺点:需要联网,经常出错,不是操作问题而是因为网络问题

安装

if("seqinr" %in% rownames(installed.packages()) == FALSE) {source("http://bioconductor.org/biocLite.R");biocLite("seqinr")}
suppressMessages(library(seqinr))
ls(‘package:seqinr‘)

####Retrieving a sequence and write into FASTA file#########

1) 选择要去fetch序列的数据库(这里已genebank为例)

choosebank()   #查看有哪些数据库
choosebank(‘genbank‘)

技术分享图片

2)一旦选择好了数据库,用query信息进行收索

BRCA1<- query("BRCA1", "SP=Homo sapiens AND K=BRCA1")

3)查看query返回的对象所有属性

attributes(BRCA1)
mynames <- getName(BRCA1)  #查看所有搜索到的名称
length(mynames)             #查看共检索到多少,写入到文档的时候可以用来用来检查 209

 技术分享图片

技术分享图片

 4)查看所有收索到的序列所包含的属性

BARC1$req

 技术分享图片

5)提取上述特定的序列(第一条序列),及注释信息

myseq <- getSequence(BRCA1$req[[1]])
annots <- getAnnot(BRCA1$req[[1]]) myseq

 技术分享图片

6) 获取所有检索到的序列,并以fasta格式写入到文件中

all_myseqs <- getSequence(BRCA1)    #所有收索到的序列
write.fasta(all_myseqs, mynames, file.out = "MyBRCA.fasta") #将所有收索到的序列写入fasta格式文件。

6)关闭接口,防止打开多个接口

closebank()

除了上述之外,如果你知道ID号码也可以根据数据库AC属性 (AC attribute)提取搜索到的特定序列

U61268<-query("BRCA1", "SP=Homo sapiens AND AC=U61268")  #也可以根据特定的ID进行搜索
attributes(U61268)
U61268$req
U61268_seq <- getSequence(U61268$req[[1]])
U61268_annots <- getAnnot(U61268$req[[1]]) 

技术分享图片

 

序列处理之seqinr(fetch)

原文:https://www.cnblogs.com/djx571/p/9661347.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!