如何用COBALT构建本地的多序列比对（Linux系统）

互联网2013-08-19

2861

COBALT 是一个蛋白的多序列比对工具，也用到RPS-Blast , BLASTP, 和PHI-BLAST等工具，并且也用到了conserved domain database (CDD) 和PROSITE protein-motif database来保证COBALT 比对结果的质量。关于在线COBALT的用法看另一篇文章” COBALT：NCBI在线蛋白多序列比对（比ClustalW还强大的工具） ”。

这里再讲解如何构建本地的COBALT多序列比对，在Linux 系统的配置步骤如下：

1．登陆Ftp：ftp://ftp.NCBI .nlm.nih.gov/pub/agarwala/cobalt/，下载cobalt里的所有文件

ncftp /pub/agarwala/cobalt > ls
cdd.aux   cdd.freq   cdd.phr  cdd.psd   cdd.psq    cobalt.linux     README
cdd.blocks  cdd.loo  dd.pin   cdd.psi  cdd.rps     patterns     README.cobalt

2．设置全局变量：setenv COBALT ，这里不设置也是可以的。下面直接用全路径就是

3．准备你所要比对的蛋白序列（Fasta 格式）

4．运行的代码如下（$COBALT也可以用全路径代替）：

       $COBALT/cobalt.linux -db $COBALT/cdd \
       -b $COBALT/cdd.blocks -i  \
       -p $COBALT/patterns -f $COBALT/cdd.freq

5．输出的结果是Fasta格式。

查看详细说明，用

$COBALT/cobalt.linux –help

更详细的说明文档看cobalt文件夹里的README.cobalt。

(责任编辑：admin)

相关产品推荐

cobalt protoporphyrin IX，阿拉丁

￥3999.90

KU-55933 是有效的 ATM 抑制剂，IC50 和 Ki 值分别为 12.9 和 2.2 nM；对 ATM 的选择性比对 DNA-PK，PI3K/PI4K，ATR 和 mTOR 高。

￥500

大肠杆菌蛋白表达| 大肠杆菌表达系统（重组蛋白表达）| 大肠杆菌重组蛋白表达| 原核(大肠杆菌) 蛋白表达与纯化服务

￥4000

α-糜蛋白酶来源于牛胰腺(TLCK处理,序列)，9004-07-3，阿拉丁

￥823.90

热喷涂粉末-粒度分布碳化钨/钴（针状）；Thermal Spray Powder - Particle Size Distribution Tungsten Carbide/Cobalt (Acicular)

￥650