【生信】自建rpsblast的蛋白结构域数据库

NCBI上BLAST的手册中对makeprofiledb的介绍：

在做rpsblast时，可以对自建的蛋白结构域数据库做比对。如：要对PKS/NRPS蛋白结构域做BLAST，如何自建一个PKS/NRPS蛋白结构域数据库？

一、流程

S1：挑出PKS/NRPS蛋白结构域的smp文件。

S2：放到一个文件夹里。

S3：写一个PKS.pn文件，里面是S1所挑的smp文件的名字，一名一行。PKS.pn同样放在这个文件夹里。

S4：参考NCBI上BLAST的手册中对makeprofiledb的介绍在该文件夹里输入

makeprofiledb -title PKS.v.1.0 -in PKS.pn -out PKS -threshold 9.82 -scale 100.0 -dbtype rps -index true

则建好了一个名为PKS的数据库

S5：输入以下命令行，则是用这个自建的PKS数据库，对20genome.fasta运行rpstblastn

rpstblastn -query 20genome.fasta -db ./domains/PKS -out op.out -evalue 0.00001 -outfmt 6

参考NCBI的BLAST手册

最后编辑于：2020.04.26 18:15:36

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。