code: def _apply_priority_filter(in_file, priority_file, data):
      
      out_file = "%s-priority%s" % utils.splitext_plus(in_file)
      if not utils.file_exists(out_file):
          with file_transaction(data, out_file) as tx_out_file:
              header = ('##INFO=<ID=EPR,Number=.,Type=String,'
                        'Description="Somatic prioritization based on external annotations, '
                        'identify as likely germline">')
              header_file = "%s-repeatheader.txt" % utils.splitext_plus(tx_out_file)[0]
              with open(header_file, "w") as out_handle:
                  out_handle.write(header)
              if "tumoronly_germline_filter" in dd.get_tools_on(data):
                  filter_cmd = ("bcftools filter -m '+' -s 'LowPriority' "
                                """-e "EPR[0] != 'pass'" |""")
              else:
                  filter_cmd = ""
              cmd = ("bcftools annotate -a {priority_file} -h {header_file} "
                     "-c CHROM,FROM,TO,REF,ALT,INFO/EPR {in_file} | "
                     "{filter_cmd} bgzip -c > {tx_out_file}")
              do.run(cmd.format(**locals()), "Run external annotation based prioritization filtering")
      vcfutils.bgzip_and_index(out_file, data["config"])
      return out_file
docs: Annotate variants with priority information and use to apply filters.