1. makeblastdb出错

你可能makeblastdb频繁报错，说是出现重复条目……
如果使用的是NCBI来源的nt.fa，请不要添加-hash_index这个参数。

2. 创建好的数据库不能使用

blastn或者blastp报了DB相关的错，网上有人提示你blastn和makeblastdb的版本不兼容。
如果你是通过conda安装的blast，那么恭喜你，你中了另一个天坑。赶紧删除这个conda环境，去NCBI的FTP服务器上下载个即开即用的ncbi-blast-版本号±x64-linux.tar.gz。conda的包有问题。

3. makeblastdb出来的数据库跑哪去了？

先cd到你想要保存的目录，然后再makeblastdb.

4. blastn或者blastp明明在命令行能运行，但是写在shell脚本里就报错offending value？

shell脚本的老问题，碰到写得稀碎的blast，就会爆炸。推荐使用python的subprocess.Popen，利用这个包传递命令不会报错，保证 -evalue 1e-6这种参数不会别识别为其他的。
网上论坛各种人让你检查是否输入的en dash还是em dash…… 没用，别耽误时间。

【天坑】BLAST比对序列时遭受的痛苦相关推荐

避开QQ聊天时遭受攻击的安全技巧(转)
避开QQ聊天时遭受攻击的安全技巧(转) QQ的密码.个人资料和聊天记录能否安全成为至关重要的问题,为了有效地防止聊天记录等本地信息的丢失和被,可以采取以下措施: 1.设置本地消息口令首先按下鼠标右键 ...
本地化BLAST+进行目标序列比对
一.在NCBI 中搜索BLAST 点击这个,然后来到这个页面,再点来到了这个下载页面,选择合适的,开始下载二.开始安装建议不要安装在C盘. 安装结束后,然后设置环境变量,在path后,添加一句 ...
用python编写程序来检查序列时否是回文序列【例如1234321就是回文序列】
a=input('请输入字符') b=a[::-1] if a==b:print("是回文序列") else:print('不是回文序列') 扩展: a = [0,1,2,3,4, ...
利用Bio.SeqIO提取基因组中特定序列时注意事项
例如: from Bio import SeqIO # 读取基因组文件 fa_seq = SeqIO.read(fasta_t_path+'\\genome.fasta', "fasta&q ...
python元祖迭代_python *式语法在迭代一个元组序列时的逻辑?
九州编程我猜想你的困惑可能是为什么:bar hello结果没有出现?因为你的records第二项写错了: ('batr', 'hello') 应该是 ('bar', 'hello'),In [1]: ...
《海豚湾》网友评论转载
http://movie.douban.com/review/2985674/ 2010-02-15 00:29:12 philip (电影是分享梦想的片刻) 我继续补充一下楼主的逻辑错误:) 1,纪 ...
第二章序列比对——Blast局部比对
第二章序列比对--Blast局部比对阅读量: 330 主要为基因组测序比对相关知识,部分内容作笔记自查使用.如有错误或遗漏还请海涵,可评论或邮箱联系. 最后修改时间:2020-04-16 16: ...
在你休息时，你的大脑运动皮层中重放习得的神经放电序列
导读以前在非人类动物中观察到的唤醒过程背后的神经激发模式的离线"重播"被认为是记忆巩固的一种机制.布朗大学(Brown University),麻省总医院(Massachuset ...
linux比对nt数据库,linux BLAST序列比对（nt/nr库）
安装本地blast序列比对软件,我们可以搜索一个查询序列定制数据库,例如想研究一个新测序的基因组,或者感兴趣的一组蛋白质序列.有时我们希望把程序插入到一个流程中,例如搜索一个大量的查询序列,例如你的测 ...

【天坑】BLAST比对序列时遭受的痛苦

1. makeblastdb出错

2. 创建好的数据库不能使用

3. makeblastdb出来的数据库跑哪去了？

4. blastn或者blastp明明在命令行能运行，但是写在shell脚本里就报错offending value？

【天坑】BLAST比对序列时遭受的痛苦相关推荐

最新文章

热门文章