code: def run_vardict(align_bams, items, ref_file, assoc_files, region=None,
                  out_file=None):
      
      items = shared.add_highdepth_genome_exclusion(items)
      if vcfutils.is_paired_analysis(align_bams, items):
          call_file = _run_vardict_paired(align_bams, items, ref_file,
                                          assoc_files, region, out_file)
      else:
          vcfutils.check_paired_problems(items)
          call_file = _run_vardict_caller(align_bams, items, ref_file,
                                          assoc_files, region, out_file)
      return call_file
docs: Run VarDict variant calling.