View source: R/filter_homologs_core_set.R
filter_homologs_core_set | R Documentation |
This function aims to retrieve a core set of best blast hits for each query sequence that is shared across all species. In other words, only query sequences
that generated blast hits in all species in the input blast_tbl
were retained.
by filtering a blast_tbl
using the following criteria.
A best hit is defined as (fulfilling all three critaria):
max(alig_length): only the hit having the longest alignment length is retained.
qcovhsp >= min_qcovhsp
: only hits that have a query coverage of at least min_qcovhsp
are retained.
max(bit_score): only the hit having the highest bit-score is retained.
filter_homologs_core_set(blast_tbl, min_qcovhsp = 50)
blast_tbl |
a BLAST table generated with |
min_qcovhsp |
minimum query coverage of the hit in percent |
Hajk-Georg Drost
filter_best_hits
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.