[−][src]Function varlociraptor::utils::filter_by_threshold

pub fn filter_by_threshold<E: Event>(
    calls: &mut Reader, 
    threshold: Option<LogProb>, 
    out: &mut Writer, 
    events: &[E], 
    vartype: &VariantType
) -> Result<(), Box<dyn Error>>

Filter a VCF record stream by a minimum threshold on the sum of posterior probabilities of a given set of Events. The threshold should be an informative false discovery rate (FDR) threshold, e.g. determined with the varlociraptor FDR control functionality.

Arguments

  • calls - BCF reader with varlociraptor calls
  • threshold - minimum threshold for the sum of posterior probabilities of the set of Events considered
  • calls - BCF writer for the filtered varlociraptor calls
  • events - the set of Events to filter on
  • vartype - the variant type to consider