code: def freebayes(in_file, ref_file, vrn_files, data):
      
      out_file = _freebayes_cutoff(in_file, data)
      #out_file = _freebayes_custom(in_file, ref_file, data)
      return out_file
docs: FreeBayes filters: cutoff-based soft filtering.