makeblastdb注意事项
来源:互联网 发布:商家如何联系淘宝客服 编辑:程序博客网 时间:2024/06/11 04:02
makeblastdb注意事项
本地blast,需要自行构建数据库。那么就要用到makeblastdb。这个命令很简单。
USAGE makeblastdb [-h] [-help] [-in input_file] [-input_type type] -dbtype molecule_type [-title database_title] [-parse_seqids] [-hash_index] [-mask_data mask_data_files] [-mask_id mask_algo_ids] [-mask_desc mask_algo_descriptions] [-gi_mask] [-gi_mask_name gi_based_mask_names] [-out database_name] [-max_file_sz number_of_bytes] [-logfile File_Name] [-taxid TaxID] [-taxid_map TaxIDMapFile] [-version]DESCRIPTION Application to create BLAST databases, version 2.4.0+Use '-help' to print detailed descriptions of command line arguments
建库命令如下:
makeblastdb -in IWGSC_v1.0_blastdb.fasta -dbtype nucl -title CS_v1.0_full -parse_seqids -out CS_v1.0_full
但是需要注意序列名字的问题,如果不注意,使用blastdbcmd调取目标序列时会出现问题。fasta格式中序列名字必须使用“|”隔开。比如序列原始名字是>chr1A ,则需要修改成>lcl|chr1A。这样就可以了,调取序列的命令是:
blastdbcmd -db CS_v1.0_full -entry chr1A -range 2-100 #因未写输出文件,默认输出到屏幕>chr1A TAAACCCTAAACCCTAAACCCTAAACCCTAAACCCCTAACCCTAAACCCTAACCAAAACCCTAAACCCTAAACCCCTAAACCCTAAACCCTAACCTAAA
如果不按上述修改名字,会出现 OID not found的错误。
当要调取的序列较大时可以将序列输出☞文件,见如下命令。然后
blastdbcmd -db /data2/Fshare/IWGSC_v1.0_Formatdb/CS_v1.0_merge -entry chr1A_part1 -range 100000-110000 -out result.fastascp mollab@192.168.3.253:/home/mollab/result.fasta ~/Desktop/ #此条命令需要在本机执行,不是在服务器端执行,并且需要输入密码,文件保存在你自己电脑的桌面上
0 0
- makeblastdb注意事项
- makeblastdb
- BLAST+中makeblastdb参数详解
- 串行的blast+2.2.25本地化: blastdbcmd,makeblastdb,blastn.简介
- 串行的blast+2.2.25本地化: blastdbcmd,makeblastdb,blastn.简介
- 注意事项!
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- 注意事项
- OSGI实战(3)-实战OSGI-翻译助手项目01
- springMVC中获取request和response
- javaWeb项目ssm框架中jdbc.properties文件加载不了或获取不到值的解决方法
- Helo World
- tomcat提示端口被占用,Tomcat v7.0 Server at localhost are already in use,tomcat端口已经被使用
- makeblastdb注意事项
- RabbitMQ(一)—RabbitMQ基础知识
- 【C#】登录成功后关闭本窗体打开主窗体
- 【PAT】1011. World Cup Betting
- 关于this的一些学习笔记
- 年薪5w和50w的职场人思考力有什么区别
- 彻底理解ThreadLocal
- vuejs教程1-简单实例
- 获得select下拉选择框的值和文本