Why Default blast search result shows fractional scores at last?
1
0
Entering edit mode
9.6 years ago

I've run Blast normally with default parameters. Alignment scores are give in fractions like 84.3. But the BLOSUM 62 matrix contains no fractional scores for any amino acids. So I would like to know how those fractional score is arrived. http://www.ncbi.nlm.nih.gov/protein/9501247?report=fasta

http://www.ncbi.nlm.nih.gov/Class/FieldGuide/BLOSUM62.txt

blast • 1.4k views
ADD COMMENT
1
Entering edit mode
9.6 years ago
pld 5.1k

http://homepages.ulb.ac.be/~dgonze/TEACHING/stat_scores.pdf

http://www.ncbi.nlm.nih.gov/blast/tutorial/Altschul-1.html

You have to remember that BLAST is not a simple sequence aligner, it tries to find you the likely best hit(s) as quickly as possible. This means BLAST trades accuracy (not precision) for speed, which means it uses many heuristics to speed the search of the sequence space up.

The score isn't an actual alignment score, but a transformation on that score to normalize on the size of the search space and the scoring system used.

If you just want a raw score, you shouldn't be using BLAST.

ADD COMMENT

Login before adding your answer.

Traffic: 3001 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6