This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Compute Blastx E-Value

Does anyone know how blastx e-values are computed? Is it the same as blastp, except the query sequence is divided by 3? I've tried using E = (m/3)n2^(-S), but I get slightly different results than from the blastx output. I need to compute blastx e-value because my blast database is split up into several parts.

Thanks in advance

blast ncbi
• 3,263 views • link
updated by Biostar Biostar • written by snayfach • snayfach

2 answers

I believe the blastx e-value uses 2*m, not m/3, because all 6-frames are included.

• 0 views • link
updated by Ram Ram • written by Bill Pearson Bill Pearson

I don't understand why, but this formula actually got me pretty close (well under an order of magnitude difference):

E = (m/6)(n)(2^-S)

-m is the sequence length in nucleotides
-n is the database size in residues
-S is the bitscore

It should be noted that blast performs a "finite size correction" in which a value is subtracted from the query and database sequence length. It's unclear what this value actually is.

http://blast.ncbi.nlm.nih.gov/Blast.cgi?CMD=Web&PAGE_TYPE=BlastNews

• 0 views • link
updated by Ram Ram • written by snayfach • snayfach

Log in to answer this question.

More posts like this