blast+本地化中blastp操作(基于PDB库)—linux[通俗易懂]
大家好,又见面了,我是你们的朋友全栈君。
blast+本地化的构建对于流程化处理大量数据序列很方便,blast+是将blast模块化,分为了蛋白质序列比对蛋白数据库(blastp)、核酸序列比对核酸数据库(blastn)、核酸序列比对蛋白质数据库(blastx)、蛋白质比对翻译后的核酸数据库(tblastn)、
翻译后的核酸序列比对翻译后的核酸数据库(tblastx)
BLAST+ 链接地址: ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST
得到最新版本的blast+压缩包:ncbi-blast-2.9.0+-x64-win64.tar.gz
对安装包进行解压:
[wangh@master Softbacks]$ tar -zxvf ncbi-blast-2.9.0+-x64-linux.tar.gz
[wangh@master Softbacks]$ cd ncbi-blast-2.9.0+
[wangh@master ncbi-blast-2.9.0+]$ vim ~/.bashrc # 对blast+进行环境配置,进入变量配置环境中后,按i或者o切换到插入(编辑模式下)输入下列路径
# ncbi-blast
export PATH=/path/ncbi-blast-2.9.0+/bin/:$PATH
#######按Esc键退出编辑环境,再输入 :wq 命令进行写入保存(w)退出(q) 按回车键(Enter)退出
[wangh@master ncbi-blast-2.9.0+]$ source ~/.bashrc #刷新你的环境配置文件,使得系统识别你刚加入的环境变量(也可以关闭当前终端再次打开,系统自动刷新环境配置文件)
blastp 的用法:
数据库从PDB中下载: ftp://ftp.wwpdb.org/pub/pdb/derived_data/pdb_seqres.txt
用cd-hit处理去冗余(解压后加变量后就能用)
$: cd-hit -i /path/pdb_seqres_simp.fasta -o PDB_simp.fasta -c 0.9
格式化数据库:
$: makeblastdb –in PDB_simp.fasta –dbtype prot –out PDB
-in:待格式化处理的fasta文件(一般是从PDB/NCBI里下载所有的相关或者整个库中的序列);
-dbtype: 数据库类型,prot或者nucl;
-out: 输出的数据库名;
蛋白质比对蛋白数据库(blastp):
blastp –query name.fasta –out name.blast –db PDB –outfmt 6 –evalue 1e-5 –num_threads 11 –max_target_seqs 6
参数说明:
-query: 输入文件路径及文件名(.fasta格式);
-out: 输出文件路径及文件名(.blast);
-db: 格式化了的数据库路径及数据库名(数据库可以从PDB/NCBI里下载所有的相关/整个库中的序列);
-outfmt:输出的文件格式,6是tabular格式对应BLAST的m8格式;
-evalue: 设置输出结果的e-value值;
-max_target_seqs:找到最大的目标的数目,也可以用-num_descriptions,tabular格式输出结果的条数;
-num_threads :线程数,并行运行时,缩短计算时间。
发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/168185.html原文链接:https://javaforall.cn
相关文章
- Linux自动获取IP地址:轻松操作指南(linux自动获取ip地址)
- Linux必不可少的分区策略(linux必须的分区)
- 追忆:忘记Linux密码的日子(忘记linux密码)
- Linux编程指南:助你应对未来(linux编程培训)
- Linux下安装Nginx的步骤详解(linux下安装nginx)
- Linux下远程文件操作的实用命令(linux远程文件命令)
- Linux查看最大连接数的简单方法(linux查看最大连接数)
- 卸载Linux下OpenJDK:一步一步操作(linux卸载openjdk)
- Linux复制虚拟机:快速、简单的方式享受实时便利(linux复制虚拟机)
- 使用Linux编写设备驱动程序(linux驱动编写)
- 操作Linux线程实现原子化操作的研究(linux线程原子)
- 掌握Linux命令行之门——学习使用通配符(linux命令通配符)
- Linux解开注释之路:快速破解技巧(linux取消注释)
- Linux查看进程日志:一次简单的操作(linux查看进程日志)
- 解决Linux终端乱码问题(linux显示乱码)
- 文本编辑玩转Linux终端:vi文本编辑操作技巧(linux终端vi)
- Linux系统关闭防火墙操作指南(关闭linux的防火墙)
- 操作系统学习Linux的分支操作系统——实现灵活操作(linux的分支)
- 破坏性操作:关闭 Linux 防火墙(关闭linux的防火墙)
- Linux系统的防火墙关闭方法(关闭linux的防火墙)
- Linux的分支:探索云计算的领域(linux的分支)
- Linux跳过繁琐操作,提高工作效率!(linux跳过)
- Linux命令大全CHM:指令一键查,linux操作轻松搞定(linux命令大全chm)
- 轻松操作:打开Linux系统的步骤(linux系统怎么打开)
- 开源Linux游戏服务器:提升服务质量与游戏体验。(游戏服务端开源linux)
- Linux脚本操作:安全删除文件(linux脚本删除文件)
- 文件Linux下重命名文件的方法(linux下重命名)
- 开启Linux微信机器人,快捷操作无忧(linux微信机器人)
- 从C语言入门到Linux编程(c语言入门linux)
- 成为Linux架构师:如何获得这个证书?(linux架构师证书)
- Linux使用指南:替换指令简单教程(linux替换指令)
- Linux磁盘挂载:让你轻松管理存储空间(linux磁盘挂载)
- 交流Linux内核:实战交流与分享(linux内核 博客)
- 用Skia在Linux上实现多彩图像处理(skia linux)
- Linux创建用户与目录:快速实现操作(linux创建用户及目录)