blastpgp生成pssm矩阵时出现错误blosum62错误

最新推荐文章于 2023-04-05 19:04:29 发布

Drone_xjw

最新推荐文章于 2023-04-05 19:04:29 发布

阅读量881

点赞数

分类专栏：生物信息学文章标签：机器学习

本文链接：https://blog.csdn.net/xjw9602/article/details/106251983

版权

生物信息学专栏收录该内容

22 篇文章 30 订阅

订阅专栏

环境

ubuntu
blast-2.2.26

错误
这个问题困扰了我一个多小时，因为我昨晚是可以跑生成pssm矩阵的blast命令

blastpgp -a 30 -i train_negative/new_negative_274_.fasta -d /home/xujunwei/dbs/uniref50 -h 0.001 -e 0.0001 -j 3 -Q pssm_files_train_negative/new_negative_274_.pssm

但是早上过来就报下面的错误。

解决办法
在blast-2.2.26的文件夹中，找到data文件夹下的blosum62文件，复制到当前执行路径下，再执行blastpgp命令就可以生成pssm矩阵了。

确定要放弃本次机会？

福利倒计时

: :

立减 ¥

普通VIP年卡可用

立即使用

Drone_xjw

关注关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

专栏目录

BLAST（2004版）本地化安装与使用，生成PSSM打分矩阵

u013313168的专栏

01-02

3677

BLAST（2004版）本地化安装与使用，生成PSSM打分矩阵一、软件安装 1、软件下载：该软件目前已更新到2.3.0+（ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/2.3.0/），我使用的是师姐给的2004年2.2.9版（ftp://ftp.ncbi.nlm.nih.gov/blast/executables/rel

psiblast生成pssm矩阵及psiblast参数详解

阿斯达克

03-28

1万+

PSI-BLAST : 敏感度更高的蛋白序列与蛋白序列之间的比对：一般用来生成pssm打分矩阵其常用命令有USAGE psiblast [-h] [-help] [-import_search_strategy filename] [-export_search_strategy filename] [-db database_name] [-dbsize num_letters] [-gili...

参与评论您还未登录，请先登录后发表或查看评论

BLOSUM62文件

05-22

POSSUM62文件--蛋白质序列打分矩阵文件。BLOSUM打分矩阵是一种在生物信息学中用于序列对比的氨基酸替换打分矩阵。BLOSUM 是“blocks substitution matrix”的缩写。它是目前常用的一种氨基酸替换打分矩阵。BLOSUM打分矩阵最早由 Steven Henikoff. 和 J.G Henikoff在他们的论文中被提出。其中，他们从BLOCKS数据库中对那些在高度保守序列中的蛋白质家族进行观察测量进而整理出了氨基酸替换的概率。他们继续使用对数胜算来计算矩阵中的分值。与PAM打分矩阵相比，BLOSUM打分矩阵的内容皆由观察得出。[1]在实际运用中，BLOSUM矩阵通常能获得更好的结果。

BLOSUM62的知识

最新发布

JEREMY的博客

04-05

3936

BLOSUM 是“blocks substitution matrix”的缩写。BLOSUM打分矩阵最早由 Steven Henikoff. 和 J.G Henikoff在他们的论文中被提出。其中，他们从BLOCKS数据库中对那些在高度保守序列中的蛋白质家族进行观察测量进而整理出了氨基酸替换的概率。他们继续使用对数胜算来计算矩阵中的分值。与PAM打分矩阵相比，BLOSUM打分矩阵的内容皆由观察得出。在实际运用中，BLOSUM矩阵通常能获得更好的结果。

用c写出程序，不打表，求出blosum62

y632481222的博客

03-21

473

Blosum62矩阵是通过对大量已知的蛋白质序列进行比对，统计它们之间的相似性得到的。具体来说，Blosum62矩阵中每个元素的值表示两个氨基酸之间的相似性得分，得分越高表示两个氨基酸之间的相似性越高。Blosum62是一种常用的蛋白质序列比对矩阵，可以用于计算两个蛋白质序列之间的相似性。要计算Blosum62矩阵，可以先定义一个二维数组来存储矩阵中每个元素的值。然后，可以通过对已知的蛋白质序列进行比对，统计它们之间的相似性得分来填充矩阵中的每个元素。

基础62 矩阵问题

qq_39338091的博客

04-23

344

62 矩阵问题作者: Turbo时间限制: 1S章节: 二维数组问题描述 : 从键盘上输入一个整数N，按以下规律输出一个(2N+1)*(2N+1)的矩阵：对角线的值为1，上半三角区域值为2，下半三角区域值为3，左半三角区域的值为4，右半三角区域的值为5。输入说明 : 你的程序需要从标准输入设备（通常为键盘）中读入多组测试数据。每组测试数据数据占一行，大于零，小于等于8。输出说明 ...

PSSM矩阵的两种生成方式

xiaobai1_1的博客

11-11

5213

一、本地生成方式——PSSM矩阵的生成（ncbi-blast-2.9.0±win64）本地生成需要自己下载ncbi-blast，还要自己动手安装及环境变量的配置；然后还要自己进行数据库构建，而一般我们需要拿来比对的数据库都比较大，下载很浪费时间。具体操作方式可以参考以下文章： https://blog.csdn.net/xiaobai1_1/article/details/103224665 二、POSSUM服务器服务器在线生成http://possum.erc.monash.edu/ POSSUM服务

linux下用psiblast批量生成pssm矩阵

褚骏逸的学习之路

01-28

2979

目录在linux上安装psiblast下载并编译用于比对的大型蛋白质数据库生成pssm批量生成方法文章中看到的手动生成pssm的方法在linux上安装psiblast 最好新建一个python环境，因为我发现conda安装blast默认的是python==3.6.11，可能会不小心把你的python版本改掉…然后你写好的代码全die了…… conda create -n blast python==3.6.11 source activate blast conda install -c bioconda

PSI-Blast最新版单机安装,批量生成Pssm打分矩阵

u013313168的专栏

08-23

6978

一、PSI-Blast在windows10系统下本地安装 1. PSI-blast最新版下载,版本：ncbi-blast-2.6.0+-win64.exe 网址：ftp://ftp.cnbi.nlm.nih.gov/blast/executables/blast+/LATEST/ 2.安装流程 2.1 点击.exe直接安装在非系统盘即可，在b

使用Blast本地数据库获得PSSM特征矩阵

Javis486的专栏

04-08

4595

1.下载数据db文件这次主要用到2中类型的蛋白质数据nr和nt，下载地址：ftp://ftp.ncbi.nlm.nih.gov/blast/db，由于nr和nt比较大，建议使用迅雷下载分割的小文件(nr64个,nt45个)，如nr.00.tar.gz 2.解压所有的.gz文件可以使用通配符解压(不知为何中途出错)， tar -xzvf *.gz 也可以使用脚本命令解压，下面给出pyt

blast常见错误|[blastpgp] WARNING: Unable to open BLOSUM62

GUET_DM_LQ的博客

05-22

938

解决方法：将BLOSUM62文件放到当前工作目录下

蛋白质功能预测中PSSM矩阵的生成

CBB_FT的博客

05-07

1万+

生成PSSM矩阵

本地BLAST的使用

duanhebo的博客

03-27

9933

本地BLAST的使用数据库的获取最简单的方法是直接到NCBI或别的网站去下载也可以将自己的序列，或与自己工作相关的序列进行整理构建成一个小型的数据库注意：以上文件格式一般可存为fasta格式构建BLAST用的数据库将已构建好的数据拷贝到你所设定的数据库所在文件夹运行cmd命令在cmd环境中输入如下所示命令 formatdb –i inseqs.fa –p F –o T –n db_...

本地blast详解

amw5181360的博客

11-21

2231

BLASTALL 用法 a.格式化序列数据库格式化序列数据库— —formatdb formatdb简单介绍: formatdb处理的都是格式为 ASN.1和 FASTA，而且不论是核苷酸序列数据库，还是蛋白质序列数据库；不论是使用Blastall ，还是Blastpgp，Mega Blast应用程序，这一步都是不可少的。 formatdb命令行参数: ...

生物信息学_序列比对（二）替换计分矩阵

金石就好的博客

07-12

3997

在序列比对时，因为空格插入不同产生了多种序列排序方案，可以通过替换计分矩阵这一计分规则来筛选出比对最大相似度的方案。核酸和蛋白质分别有不同的替换计分矩阵。核酸等价矩阵相同核苷酸之间的匹配得分为1，不同核苷酸间的替换为0。转换-颠换矩阵核酸的碱基按照结构可分为两类：嘌呤（两个环）嘧啶（一个环）不同结构之间替换称为颠换，相同结构替换为转换在进化中转换频率远大于颠换，所以转换得分为-1，颠换为-5。 BLAST矩阵经过大量实际比较得，核苷酸相同的+5，反之为-4，则效果较好，该

最长连续子序列偏移_生物信息学初识篇——第二章：序列比对（1）

weixin_39827315的博客

10-21

1065

点击查看往期内容回顾生物信息学初识篇——第一章：生物数据库第二章：序列比对(1)生物序列有自己的书写格式，而且格式很多。不同的处理软件会用到不同的格式，但是最常用的，大多数软件都识别的格式是 FASTA 格式。我们用一致度(identity)和相似度(similarity)这两个指标来定量描述序列有多相似。一.替换记分矩阵替换记分矩阵是反映残基之间相互替换率的矩阵。也就是说，它描述了残基...

PSI-BLAST参数解释