This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Blast E-Value To Database Size

Hi,

I'm thinking of splitting the database to smaller chunks. And, blast my sequences against them each on a separate process. My only concern is the results (which I will merge later).

Would the resulting e-value be affected by database content when smaller subsets are used? I have a hunch that it would not matter when all the subset results later becomes concatenated. Please correct me if I'm wrong.

blast statistics
• 17,642 views • link
updated by Neilfws Neilfws • written by Hranjeev Hranjeev

1 answer

The statistics of BLAST scores are described in this article. It's quite mathematics-heavy, but also quite readable; just take your time and re-read several times.

The short answer is that yes, e-values are dependent on database size. If you think about it intuitively, there's a higher probability of finding a match in a large database than in a smaller database.

That said, it is possible to re-calculate e-values by combining the results when the database is split. This is implemented in, for example, mpiBLAST. It would be a good idea to study their website, code and publication to see how they handle the problem.

See also the discussion of recalculating e-value in this paper or do a quick web search for "BLAST split database e-value calculate" - it's quite a widely-discussed issue.

• 139 views • link
written by Neilfws Neilfws

thanks much appreciated

• 0 views • link
written by Sequer • Sequer

Totally agree with the answer but you can set manually the database size using the parameter "-z". On that way, you can split the db file into smaller pieces, make your queries and then merge results.

• 0 views • link
written by scapella scapella

I think you also need to set the number of sequences in the database (N) to calculate the edge adjustment parameter (l or "ell"). The adjustment is done for you if you use NOBLAST

• 0 views • link
written by colinDotAIBN • colinDotAIBN

Oops my hunch was wrong. Anyways any easy way to count the number of letters, N (total letters) of a database?

• 0 views • link
written by Hranjeev Hranjeev

Please validate if the nr database atm is 5784003470 letters in size

• 0 views • link
written by Hranjeev Hranjeev

try grep -v '^>' something.fasta | grep -o [ACTGNactg] | wc -l for fasta files before building database

• 0 views • link
written by Biomonika (Noolean) Biomonika (Noolean)

When I said "database size" I refered to the total number of sequences in your database, I didn't refer to the total number of residues on it.

• 0 views • link
written by scapella scapella

Log in to answer this question.

More posts like this