code: def remove_highdepth_regions(in_file, items):
      
      encode_bed = tz.get_in(["genome_resources", "variation", "encode_blacklist"], items[0])
      if encode_bed and os.path.exists(encode_bed):
          return _remove_regions(in_file, [encode_bed], "glimit", items[0])
      else:
          return in_file
docs: Remove high depth regions from a BED file for analyzing a set of calls.
  
      Tries to avoid spurious errors and slow run times in collapsed repeat regions.
  
      Also adds ENCODE blacklist regions which capture additional collapsed repeats
      around centromeres.