HHblits使用相同的数据库，得出的结果明显差于HHpred网页版的结果。该如何使用HHsuite命令行得到与网页端类似的效果 - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

1387: 浏览

HHblits使用相同的数据库，得出的结果明显差于HHpred网页版的结果。该如何使用HHsuite命令行得到与网页端类似的效果

生信问题

如题。

同样一条序列，使用HHblits在PDB70数据库中比对，无法得出高质量的比对，使用hhsearch也得出同样的结果。

但使用网页端的HHpred就可以注释出来。

并且我手动确认了，网页端HHpred注释出的domain id在我本地运行hhblits所用的pdb70数据库中存在。

但是它就是比对不出来。查看github官方文档也没有提到具体该如何使用HHsuite去复现HHpred的结果。

故来此请教该如何使用命令行版本的HHsuite，实现与网页端HHpred相同或类似的比对效果。

好问题 0 评论收藏举报

2 回答

chatGPT机器人机器人用户来自于: IANA
2023-05-26 18:31

首先，需要确认使用HHpred网页版得到的结果中使用的是哪个数据库。HHpred网页版默认使用了UniRef90数据库进行比对，而HHblits则使用PDB70数据库。这可能是造成差异的原因之一。如果需要得到与HHpred类似的结果，可以尝试使用UniRef90数据库进行比对。其次，可以尝试调整HHblits的一些参数，以改善比对效果。例如，可以尝试调整迭代次数、E-value阈值、gap扩展惩罚等参数，以找到最适合的参数组合。可以使用HHblits的--help命令查看各参数的含义和默认值。另外，HHpred网页版还使用了PSI-BLAST算法进行预处理，以进一步提高比对效果。可以尝试在HHblits之前使用PSI-BLAST进行预处理，以得到更好的结果。最后，如果以上方法都无法改善比对效果，可以考虑使用其他的比对工具，例如HMMER或者BLAST，以找到更适合的工具和参数组合。总之，要得到与HHpred类似的比对效果，需要综合考虑数据库、算法、参数等因素，并进行适当的调整和优化。

赞同 0 1评论

孟浩巍超级管理员用户来自于: 北京市
2023-05-26 19:57

我有几个建议：

去读一下官网的mannual，然后看看是否有不同；
提交的时候，看看有没有暴露后端提交的代码，因为网页版本质上还是命令行计算，看看他们的参数怎么设置的。
实在不行，写邮件给作者。

赞同 0 0评论

关于作者

: 聪明的生信汪 注册会员
这家伙很懒，还没有设置简介

0: 回答

0: 文章

1: 问题

问题动态

发布时间: 2023-05-26 18:21

更新时间: 2023-05-26 19:57

关注人数: 2 人关注

相关问题

两次运行同一个FQ文件产生的bam文件分析结果存在差异: 903 浏览 2 关注 2 回答 0 评论

使用bwa index 时遇到问题: 853 浏览 2 关注 1 回答 2 评论

在单细胞的细胞通讯分析软件cellphonedb的结果解读: 826 浏览 2 关注 1 回答 0 评论

菌侵染植物的转录组分析中，若菌的一个基因在当前时间没有引起植物基因的变化，过了三个小时后才引起植物基因的变化，如何分析这种现象或如何将这一类基因关联起来，有没有类似的文献推荐？: 846 浏览 2 关注 1 回答 0 评论

如何在Ubuntu 14.04 server版上安装Rstudio，并在网页上实现绘图功能？: 1152 浏览 4 关注 3 回答 0 评论

当LIMMA计算计算结果中的调整P值都大于0.05时，可否使用p值<0.05+fc>2来筛选差异表达基因？: 1030 浏览 2 关注 1 回答 0 评论

qpAdm结果解读: 944 浏览 2 关注 1 回答 0 评论

使用deeptools生成ChIP-seq信号热图与谱图: 2358 浏览 3 关注 2 回答 0 评论

seqtk的使用技巧，可以处理那些序列问题？: 869 浏览 1 关注 1 回答 0 评论

htseq使用出现故障: 699 浏览 2 关注 1 回答 0 评论

推荐内容

问 grep 能不能指定匹配第多少列？: 1002 浏览 3 关注 2 回答 1 评论

问 de novo mutation检测和判断: 622 浏览 2 关注 1 回答 0 评论

问 SPSS报错: 988 浏览 2 关注 1 回答 0 评论

问使用bedtools提取序列时的deyond length问题；bed负值: 1423 浏览 2 关注 2 回答 0 评论

问不同样品中寻找特异的OTU: 1014 浏览 2 关注 2 回答 0 评论

问染色体号是罗马数字怎么写sh循环: 766 浏览 1 关注 1 回答 0 评论

问 RNA Seq中什么样的reads是chimeric reads？是不是多数的chimeric reads都是circRNA的？: 803 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024