Question: Blast raw score calculation
gravatar for cristianrohr768
5.1 years ago by
cristianrohr76830 wrote:


I search for the past two days, trying to find out how the blast raw score is calculated, i read i lot, read the answers here, but i can't found the right answer to my problem

I'm using a 240 aa query, blastp agains refseq, the first result is


Range 1: 140 to 379

Score Expect Method Identities Positives Gaps
462 bits(1190) 2e-160 Compositional matrix adjust. 239/240(99%) 240/240(100%) 0/240(0%)

Raw score 1190, with 1 mismatch

the protein is from H sapiens, and this result appears in the 4 place


Range 1: 1 to 240

Score Expect Method Identities Positives Gaps
450 bits(1157) 6e-158 Compositional matrix adjust. 240/240(100%) 240/240(100%) 0/240(0%)

The raw score is 1157,

how is this possible, that a perfect match has a lower raw score??


blast raw score calculation • 2.9k views
ADD COMMENTlink modified 4.4 years ago by Biostar ♦♦ 20 • written 5.1 years ago by cristianrohr76830

The aligned sequences are different: "Range 1: 140 to 379" and "Range 1: 1 to 240". Have a look at the diagonal of the BLOSUM62 matrix and you will understand the reason.

ADD REPLYlink written 4.4 years ago by lh332k
gravatar for Damian Kao
5.1 years ago by
Damian Kao15k
Damian Kao15k wrote:

Two possible reasons. 1) The raw score is not length normalized. Longer sequences could have bigger raw score values. 2) The raw score is also calculated from a substitution matrix. The purpose of the matrix is to give a score based on the likelihood of an amino acid being substituted for another. Some substitutions will have better score than others. So depending on the amino acid content of your sequence, you can have different scores.



By the way, if you want to learn how blast works. Ian Korf's BLAST essentials book is probably the most comprehensive resource:

ADD COMMENTlink modified 5.1 years ago • written 5.1 years ago by Damian Kao15k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1680 users visited in the last hour