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

登录

该问题已被锁定！

2: 关注

3106: 浏览

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评论

关于作者

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

0: 回答

0: 文章

1: 问题

问题动态

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

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

关注人数: 2 人关注

相关问题

请问一下，使用conda install deeptools之后，出现如下报错，如何解决呢？: 2291 浏览 1 关注 1 回答 0 评论

使用MAS5.0处理过的芯片数据，能不能直接取对数后用limma包找差异表达基因，为什么不能？: 1918 浏览 2 关注 1 回答 0 评论

如何针对affy表达谱芯片计算结果分析不同样本之间是否表达差异: 1454 浏览 2 关注 1 回答 0 评论

运行roary软件不出结果，命令行用了这种: 2861 浏览 1 关注 2 回答 0 评论

使用cutadapt去掉adaptor后，fastqc结果的sequence distribution报错: 2366 浏览 2 关注 2 回答 0 评论

GAPIT包导出的GWAS结果如何添加新的阈值线？以及GAPIT的结果文件中的nobs、H&B.P.Value、Effect分别是什么意思？: 2567 浏览 2 关注 1 回答 0 评论

linux系统中使用fastqc报错: 3036 浏览 3 关注 3 回答 0 评论

blastp后，我想提取score的序列，请问使用python应该怎么写这个脚本: 1716 浏览 2 关注 1 回答 0 评论

SMRT三代测序Blasr结果coverage的分析以及可视化？: 1790 浏览 2 关注 1 回答 0 评论

linux下使用convert出现报错，可能是什么原因？如何解决？: 2317 浏览 2 关注 1 回答 0 评论

推荐内容

问关于基因间的相关性分析: 2306 浏览 2 关注 2 回答 0 评论

问 rna-seq数据校正: 2098 浏览 2 关注 2 回答 0 评论

问关于GO分析的问题: 1536 浏览 2 关注 1 回答 0 评论

问 3D基因组里compartment里一般是包含好多TAD的，但这图为什么compartment数量比TAD多这么多呢: 2238 浏览 2 关注 2 回答 0 评论

问关于ceRNA网络构建的后续分析有哪些？: 1535 浏览 2 关注 1 回答 0 评论

问关于生存分析的问题: 1063 浏览 1 关注 0 回答 0 评论

问如何按一个列表对基因型文件进行过滤，剔除不需要的样本？: 2264 浏览 2 关注 1 回答 0 评论

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

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

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025