Skip to contents

Given a BLAST output, filters query-subject pairs such that only HSPs with a coverage greater than mncvrg (specified in the pipeline parameters) remain. Filters both: query-subject and subject-query pairs, if one of the coverages is insufficient. HSP coverage is obtained from the BLAST column qcovs.

Usage

blast_filter(blast_res, ps, lvl = 3)

Arguments

blast_res

BLAST results

ps

Parameters list, generated with parameters()

lvl

Integer, number of message indentations indicating code depth.

Value

data.frame blast res