Question: Understanding Best-Hits filtering algorithm in BLAST+ applications
0
gravatar for Denis
16 months ago by
Denis40
Denis40 wrote:

Hi! I not completely understand Best-Hits filtering algorithm implemented in BLAST+ aplications. As far i can understand basic idea is to filter out short HSPs with lower E-value and save only relatively long HSPs. Am i right? From the BLAST documentation i have :

"The Best-Hit filtering algorithm is designed for use in applications that are searching for only the best matches for each query region reporting matches. Its -best_hit_overhang parameter, H, controls when an HSP is considered short enough to be filtered due to presence of another HSP. For each HSP A that is filtered, there exists another HSP B such that the query region of HSP A extends each end of the query region of HSP B by at most H times the length of the query region for B. Additional requirements that must also be met in order to filter A on account of B are:

i. evalue(A) >= evalue(B)

ii. score(A)/length(A) < (1.0 – score_edge) * score(B)/length(B)

We consider 0.1 to 0.25 to be an acceptable range for the -best_hit_overhang parameter and 0.05 to 0.25 to be an acceptable range for the -best_hit_score_edge parameter."

But i have several questions:

  1. HSP A and HSP B. Obviously they must related to the same query. But must both HSPs represent only two different hits? I mean, may HSP A and HSP B originate from the same hit sequence?

  2. E-value of HSP A must be greater those of HSP B, and opposite situation for bit score (bit score of HSP A< bit score HSP B). Right? I'm wondering what is about length of HSP A and HSP B?

  3. I can not understand the following sentence from the documentation: "For each HSP A that is filtered, there exists another HSP B such that the query region of HSP A extends each end of the query region of HSP B by at most H times the length of the query region for B." Is there any other explanation (may be graphical)?

blast alignment • 720 views
ADD COMMENTlink modified 16 months ago by RamRS21k • written 16 months ago by Denis40
Please log in to add an answer.

Help
Access

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