#!/usr/bin/env bash

# bcftools_annotate add_sgdemux
# 
# This wrapper script is auto-generated by viash 0.9.0 and is thus a derivative
# work thereof. This software comes with ABSOLUTELY NO WARRANTY from Data
# Intuitive.
# 
# The component may contain files which fall under a different license. The
# authors of this component should specify the license in the header of such
# files, or include a separate license file detailing the licenses of all included
# files.
# 
# Component authors:
#  * Theodoro Gasperin Terra Camargo (author)

set -e

if [ -z "$VIASH_TEMP" ]; then
  VIASH_TEMP=${VIASH_TEMP:-$VIASH_TMPDIR}
  VIASH_TEMP=${VIASH_TEMP:-$VIASH_TEMPDIR}
  VIASH_TEMP=${VIASH_TEMP:-$VIASH_TMP}
  VIASH_TEMP=${VIASH_TEMP:-$TMPDIR}
  VIASH_TEMP=${VIASH_TEMP:-$TMP}
  VIASH_TEMP=${VIASH_TEMP:-$TEMPDIR}
  VIASH_TEMP=${VIASH_TEMP:-$TEMP}
  VIASH_TEMP=${VIASH_TEMP:-/tmp}
fi

# define helper functions
# ViashQuote: put quotes around non flag values
# $1     : unquoted string
# return : possibly quoted string
# examples:
#   ViashQuote --foo      # returns --foo
#   ViashQuote bar        # returns 'bar'
#   Viashquote --foo=bar  # returns --foo='bar'
function ViashQuote {
  if [[ "$1" =~ ^-+[a-zA-Z0-9_\-]+=.+$ ]]; then
    echo "$1" | sed "s#=\(.*\)#='\1'#"
  elif [[ "$1" =~ ^-+[a-zA-Z0-9_\-]+$ ]]; then
    echo "$1"
  else
    echo "'$1'"
  fi
}
# ViashRemoveFlags: Remove leading flag
# $1     : string with a possible leading flag
# return : string without possible leading flag
# examples:
#   ViashRemoveFlags --foo=bar  # returns bar
function ViashRemoveFlags {
  echo "$1" | sed 's/^--*[a-zA-Z0-9_\-]*=//'
}
# ViashSourceDir: return the path of a bash file, following symlinks
# usage   : ViashSourceDir ${BASH_SOURCE[0]}
# $1      : Should always be set to ${BASH_SOURCE[0]}
# returns : The absolute path of the bash file
function ViashSourceDir {
  local source="$1"
  while [ -h "$source" ]; do
    local dir="$( cd -P "$( dirname "$source" )" >/dev/null 2>&1 && pwd )"
    source="$(readlink "$source")"
    [[ $source != /* ]] && source="$dir/$source"
  done
  cd -P "$( dirname "$source" )" >/dev/null 2>&1 && pwd
}
# ViashFindTargetDir: return the path of the '.build.yaml' file, following symlinks
# usage   : ViashFindTargetDir 'ScriptPath'
# $1      : The location from where to start the upward search
# returns : The absolute path of the '.build.yaml' file
function ViashFindTargetDir {
  local source="$1"
  while [[ "$source" != "" && ! -e "$source/.build.yaml" ]]; do
    source=${source%/*}
  done
  echo $source
}
# see https://en.wikipedia.org/wiki/Syslog#Severity_level
VIASH_LOGCODE_EMERGENCY=0
VIASH_LOGCODE_ALERT=1
VIASH_LOGCODE_CRITICAL=2
VIASH_LOGCODE_ERROR=3
VIASH_LOGCODE_WARNING=4
VIASH_LOGCODE_NOTICE=5
VIASH_LOGCODE_INFO=6
VIASH_LOGCODE_DEBUG=7
VIASH_VERBOSITY=$VIASH_LOGCODE_NOTICE

# ViashLog: Log events depending on the verbosity level
# usage: ViashLog 1 alert Oh no something went wrong!
# $1: required verbosity level
# $2: display tag
# $3+: messages to display
# stdout: Your input, prepended by '[$2] '.
function ViashLog {
  local required_level="$1"
  local display_tag="$2"
  shift 2
  if [ $VIASH_VERBOSITY -ge $required_level ]; then
    >&2 echo "[$display_tag]" "$@"
  fi
}

# ViashEmergency: log events when the system is unstable
# usage: ViashEmergency Oh no something went wrong.
# stdout: Your input, prepended by '[emergency] '.
function ViashEmergency {
  ViashLog $VIASH_LOGCODE_EMERGENCY emergency "$@"
}

# ViashAlert: log events when actions must be taken immediately (e.g. corrupted system database)
# usage: ViashAlert Oh no something went wrong.
# stdout: Your input, prepended by '[alert] '.
function ViashAlert {
  ViashLog $VIASH_LOGCODE_ALERT alert "$@"
}

# ViashCritical: log events when a critical condition occurs
# usage: ViashCritical Oh no something went wrong.
# stdout: Your input, prepended by '[critical] '.
function ViashCritical {
  ViashLog $VIASH_LOGCODE_CRITICAL critical "$@"
}

# ViashError: log events when an error condition occurs
# usage: ViashError Oh no something went wrong.
# stdout: Your input, prepended by '[error] '.
function ViashError {
  ViashLog $VIASH_LOGCODE_ERROR error "$@"
}

# ViashWarning: log potentially abnormal events
# usage: ViashWarning Something may have gone wrong.
# stdout: Your input, prepended by '[warning] '.
function ViashWarning {
  ViashLog $VIASH_LOGCODE_WARNING warning "$@"
}

# ViashNotice: log significant but normal events
# usage: ViashNotice This just happened.
# stdout: Your input, prepended by '[notice] '.
function ViashNotice {
  ViashLog $VIASH_LOGCODE_NOTICE notice "$@"
}

# ViashInfo: log normal events
# usage: ViashInfo This just happened.
# stdout: Your input, prepended by '[info] '.
function ViashInfo {
  ViashLog $VIASH_LOGCODE_INFO info "$@"
}

# ViashDebug: log all events, for debugging purposes
# usage: ViashDebug This just happened.
# stdout: Your input, prepended by '[debug] '.
function ViashDebug {
  ViashLog $VIASH_LOGCODE_DEBUG debug "$@"
}

# find source folder of this component
VIASH_META_RESOURCES_DIR=`ViashSourceDir ${BASH_SOURCE[0]}`

# find the root of the built components & dependencies
VIASH_TARGET_DIR=`ViashFindTargetDir $VIASH_META_RESOURCES_DIR`

# define meta fields
VIASH_META_NAME="bcftools_annotate"
VIASH_META_FUNCTIONALITY_NAME="bcftools_annotate"
VIASH_META_EXECUTABLE="$VIASH_META_RESOURCES_DIR/$VIASH_META_NAME"
VIASH_META_CONFIG="$VIASH_META_RESOURCES_DIR/.config.vsh.yaml"
VIASH_META_TEMP_DIR="$VIASH_TEMP"


# ViashHelp: Display helpful explanation about this executable
function ViashHelp {
  echo "bcftools_annotate add_sgdemux"
  echo ""
  echo "Add or remove annotations from a VCF/BCF file."
  echo ""
  echo "Inputs:"
  echo "    -i, --input"
  echo "        type: file, required parameter, multiple values allowed, file must exist"
  echo "        Input VCF/BCF file."
  echo ""
  echo "Outputs:"
  echo "    -o, --output"
  echo "        type: file, required parameter, output, file must exist"
  echo "        Output annotated file."
  echo ""
  echo "Options:"
  echo "    For examples on how to use use bcftools annotate see"
  echo "    http://samtools.github.io/bcftools/howtos/annotate.html."
  echo "    For more details on the options see"
  echo "    https://samtools.github.io/bcftools/bcftools.html#annotate."
  echo ""
  echo "    --a, --annotations"
  echo "        type: file, file must exist"
  echo "        VCF file or tabix-indexed FILE with annotations: CHR\\tPOS[\\tVALUE]+ ."
  echo ""
  echo "    --c, --columns"
  echo "        type: string"
  echo "        List of columns in the annotation file, e.g."
  echo "        CHROM,POS,REF,ALT,-,INFO/TAG."
  echo "        See man page for details."
  echo ""
  echo "    --C, --columns_file"
  echo "        type: file, file must exist"
  echo "        Read -c columns from FILE, one name per row, with optional --merge_logic"
  echo "        TYPE: NAME[ TYPE]."
  echo ""
  echo "    --e, --exclude"
  echo "        type: string"
  echo "        example: QUAL >= 30 && DP >= 10"
  echo "        Exclude sites for which the expression is true."
  echo "        See https://samtools.github.io/bcftools/bcftools.html#expressions for"
  echo "        details."
  echo ""
  echo "    --force"
  echo "        type: boolean_true"
  echo "        continue even when parsing errors, such as undefined tags, are"
  echo "        encountered."
  echo "        Note this can be an unsafe operation and can result in corrupted BCF"
  echo "        files."
  echo "        If this option is used, make sure to sanity check the result thoroughly."
  echo ""
  echo "    --H, --header_line"
  echo "        type: string"
  echo "        Header line which should be appended to the VCF header, can be given"
  echo "        multiple times."
  echo ""
  echo "    --h, --header_lines"
  echo "        type: file, file must exist"
  echo "        File with header lines to append to the VCF header."
  echo "        For example:"
  echo "          ##INFO=<ID=NUMERIC_TAG,Number=1,Type=Integer,Description=\"Example"
  echo "        header line\">"
  echo "          ##INFO=<ID=STRING_TAG,Number=1,Type=String,Description=\"Yet another"
  echo "        header line\">"
  echo ""
  echo "    --I, --set_id"
  echo "        type: string"
  echo "        Set ID column using a \`bcftools query\`-like expression, see man page for"
  echo "        details."
  echo ""
  echo "    --include"
  echo "        type: string"
  echo "        example: QUAL >= 30 && DP >= 10"
  echo "        Select sites for which the expression is true."
  echo "        See https://samtools.github.io/bcftools/bcftools.html#expressions for"
  echo "        details."
  echo ""
  echo "    --k, --keep_sites"
  echo "        type: boolean_true"
  echo "        Leave --include/--exclude sites unchanged instead of discarding them."
  echo ""
  echo "    --l, --merge_logic"
  echo "        type: string"
  echo "        When multiple regions overlap a single record, this option defines how"
  echo "        to treat multiple annotation values."
  echo "        See man page for more details."
  echo ""
  echo "    --m, --mark_sites"
  echo "        type: string"
  echo "        Annotate sites which are present (\"+\") or absent (\"-\") in the -a file"
  echo "        with a new INFO/TAG flag."
  echo ""
  echo "    --min_overlap"
  echo "        type: string"
  echo "        Minimum overlap required as a fraction of the variant in the annotation"
  echo "        -a file (ANN),"
  echo "        in the target VCF file (:VCF), or both for reciprocal overlap (ANN:VCF)."
  echo "        By default overlaps of arbitrary length are sufficient."
  echo "        The option can be used only with the tab-delimited annotation -a file"
  echo "        and with BEG and END columns present."
  echo ""
  echo "    --no_version"
  echo "        type: boolean_true"
  echo "        Do not append version and command line information to the output VCF"
  echo "        header."
  echo ""
  echo "    --O, --output_type"
  echo "        type: string"
  echo "        choices: [ u, z, b, v ]"
  echo "        Output type:"
  echo "          u: uncompressed BCF"
  echo "          z: compressed VCF"
  echo "          b: compressed BCF"
  echo "          v: uncompressed VCF"
  echo ""
  echo "    --pair_logic"
  echo "        type: string"
  echo "        choices: [ snps, indels, both, all, some, exact ]"
  echo "        Controls how to match records from the annotation file to the target"
  echo "        VCF."
  echo "        Effective only when -a is a VCF or BCF file."
  echo "        The option replaces the former uninuitive --collapse."
  echo "        See Common Options for more."
  echo ""
  echo "    --r, --regions"
  echo "        type: string"
  echo "        example: 20:1000000-2000000"
  echo "        Restrict to comma-separated list of regions."
  echo "        Following formats are supported: chr|chr:pos|chr:beg-end|chr:beg-[,…​]."
  echo ""
  echo "    --R, --regions_file"
  echo "        type: file, file must exist"
  echo "        Restrict to regions listed in a file."
  echo "        Regions can be specified either on a VCF, BED, or tab-delimited file"
  echo "        (the default)."
  echo "        For more information check manual."
  echo ""
  echo "    --regions_overlap"
  echo "        type: string"
  echo "        choices: [ pos, record, variant, 0, 1, 2 ]"
  echo "        This option controls how overlapping records are determined:"
  echo "        set to 'pos' or '0' if the VCF record has to have POS inside a region"
  echo "        (this corresponds to the default behavior of -t/-T);"
  echo "        set to 'record' or '1' if also overlapping records with POS outside a"
  echo "        region should be included (this is the default behavior of -r/-R,"
  echo "        and includes indels with POS at the end of a region, which are"
  echo "        technically outside the region);"
  echo "        or set to 'variant' or '2' to include only true overlapping variation"
  echo "        (compare the full VCF representation \"TA>T-\" vs the true sequence"
  echo "        variation \"A>-\")."
  echo ""
  echo "    --rename_annotations"
  echo "        type: file, file must exist"
  echo "        Rename annotations: TYPE/old\\tnew, where TYPE is one of"
  echo "        FILTER,INFO,FORMAT."
  echo ""
  echo "    --rename_chromosomes"
  echo "        type: file, file must exist"
  echo "        Rename chromosomes according to the map in file, with \"old_name"
  echo "        new_name\\n\" pairs"
  echo "        separated by whitespaces, each on a separate line."
  echo ""
  echo "    --samples"
  echo "        type: string"
  echo "        Subset of samples to annotate."
  echo "        See also"
  echo "        https://samtools.github.io/bcftools/bcftools.html#common_options."
  echo ""
  echo "    --samples_file"
  echo "        type: file, file must exist"
  echo "        Subset of samples to annotate in file format."
  echo "        See also"
  echo "        https://samtools.github.io/bcftools/bcftools.html#common_options."
  echo ""
  echo "    --single_overlaps"
  echo "        type: boolean_true"
  echo "        Use this option to keep memory requirements low with very large"
  echo "        annotation files."
  echo "        Note, however, that this comes at a cost, only single overlapping"
  echo "        intervals are considered in this mode."
  echo "        This was the default mode until the commit af6f0c9 (Feb 24 2019)."
  echo ""
  echo "    --x, --remove"
  echo "        type: string"
  echo "        List of annotations to remove."
  echo "        Use \"FILTER\" to remove all filters or \"FILTER/SomeFilter\" to remove a"
  echo "        specific filter."
  echo "        Similarly, \"INFO\" can be used to remove all INFO tags and \"FORMAT\" to"
  echo "        remove all FORMAT tags except GT."
  echo "        To remove all INFO tags except \"FOO\" and \"BAR\", use \"^INFO/FOO,INFO/BAR\""
  echo "        (and similarly for FORMAT and FILTER)."
  echo "        \"INFO\" can be abbreviated to \"INF\" and \"FORMAT\" to \"FMT\"."
}

# initialise variables
VIASH_MODE='run'
VIASH_ENGINE_ID='docker'

######## Helper functions for setting up Docker images for viash ########
# expects: ViashDockerBuild

# ViashDockerInstallationCheck: check whether Docker is installed correctly
#
# examples:
#   ViashDockerInstallationCheck
function ViashDockerInstallationCheck {
  ViashDebug "Checking whether Docker is installed"
  if [ ! command -v docker &> /dev/null ]; then
    ViashCritical "Docker doesn't seem to be installed. See 'https://docs.docker.com/get-docker/' for instructions."
    exit 1
  fi

  ViashDebug "Checking whether the Docker daemon is running"
  local save=$-; set +e
  local docker_version=$(docker version --format '{{.Client.APIVersion}}' 2> /dev/null)
  local out=$?
  [[ $save =~ e ]] && set -e
  if [ $out -ne 0 ]; then
    ViashCritical "Docker daemon does not seem to be running. Try one of the following:"
    ViashCritical "- Try running 'dockerd' in the command line"
    ViashCritical "- See https://docs.docker.com/config/daemon/"
    exit 1
  fi
}

# ViashDockerRemoteTagCheck: check whether a Docker image is available 
# on a remote. Assumes `docker login` has been performed, if relevant.
#
# $1                  : image identifier with format `[registry/]image[:tag]`
# exit code $?        : whether or not the image was found
# examples:
#   ViashDockerRemoteTagCheck python:latest
#   echo $?                                     # returns '0'
#   ViashDockerRemoteTagCheck sdaizudceahifu
#   echo $?                                     # returns '1'
function ViashDockerRemoteTagCheck {
  docker manifest inspect $1 > /dev/null 2> /dev/null
}

# ViashDockerLocalTagCheck: check whether a Docker image is available locally
#
# $1                  : image identifier with format `[registry/]image[:tag]`
# exit code $?        : whether or not the image was found
# examples:
#   docker pull python:latest
#   ViashDockerLocalTagCheck python:latest
#   echo $?                                     # returns '0'
#   ViashDockerLocalTagCheck sdaizudceahifu
#   echo $?                                     # returns '1'
function ViashDockerLocalTagCheck {
  [ -n "$(docker images -q $1)" ]
}

# ViashDockerPull: pull a Docker image
#
# $1                  : image identifier with format `[registry/]image[:tag]`
# exit code $?        : whether or not the image was found
# examples:
#   ViashDockerPull python:latest
#   echo $?                                     # returns '0'
#   ViashDockerPull sdaizudceahifu
#   echo $?                                     # returns '1'
function ViashDockerPull {
  ViashNotice "Checking if Docker image is available at '$1'"
  if [ $VIASH_VERBOSITY -ge $VIASH_LOGCODE_INFO ]; then
    docker pull $1 && return 0 || return 1
  else
    local save=$-; set +e
    docker pull $1 2> /dev/null > /dev/null
    local out=$?
    [[ $save =~ e ]] && set -e
    if [ $out -ne 0 ]; then
      ViashWarning "Could not pull from '$1'. Docker image doesn't exist or is not accessible."
    fi
    return $out
  fi
}

# ViashDockerPush: push a Docker image
#
# $1                  : image identifier with format `[registry/]image[:tag]`
# exit code $?        : whether or not the image was found
# examples:
#   ViashDockerPush python:latest
#   echo $?                                     # returns '0'
#   ViashDockerPush sdaizudceahifu
#   echo $?                                     # returns '1'
function ViashDockerPush {
  ViashNotice "Pushing image to '$1'"
  local save=$-; set +e
  local out
  if [ $VIASH_VERBOSITY -ge $VIASH_LOGCODE_INFO ]; then
    docker push $1
    out=$?
  else
    docker push $1 2> /dev/null > /dev/null
    out=$?
  fi
  [[ $save =~ e ]] && set -e
  if [ $out -eq 0 ]; then
    ViashNotice "Container '$1' push succeeded."
  else
    ViashError "Container '$1' push errored. You might not be logged in or have the necessary permissions."
  fi
  return $out
}

# ViashDockerPullElseBuild: pull a Docker image, else build it
#
# $1                  : image identifier with format `[registry/]image[:tag]`
# ViashDockerBuild    : a Bash function which builds a docker image, takes image identifier as argument.
# examples:
#   ViashDockerPullElseBuild mynewcomponent
function ViashDockerPullElseBuild {
  local save=$-; set +e
  ViashDockerPull $1
  local out=$?
  [[ $save =~ e ]] && set -e
  if [ $out -ne 0 ]; then
    ViashDockerBuild $@
  fi
}

# ViashDockerSetup: create a Docker image, according to specified docker setup strategy
#
# $1          : image identifier with format `[registry/]image[:tag]`
# $2          : docker setup strategy, see DockerSetupStrategy.scala
# examples:
#   ViashDockerSetup mynewcomponent alwaysbuild
function ViashDockerSetup {
  local image_id="$1"
  local setup_strategy="$2"
  if [ "$setup_strategy" == "alwaysbuild" -o "$setup_strategy" == "build" -o "$setup_strategy" == "b" ]; then
    ViashDockerBuild $image_id --no-cache $(ViashDockerBuildArgs "$engine_id")
  elif [ "$setup_strategy" == "alwayspull" -o "$setup_strategy" == "pull" -o "$setup_strategy" == "p" ]; then
    ViashDockerPull $image_id
  elif [ "$setup_strategy" == "alwayspullelsebuild" -o "$setup_strategy" == "pullelsebuild" ]; then
    ViashDockerPullElseBuild $image_id --no-cache $(ViashDockerBuildArgs "$engine_id")
  elif [ "$setup_strategy" == "alwayspullelsecachedbuild" -o "$setup_strategy" == "pullelsecachedbuild" ]; then
    ViashDockerPullElseBuild $image_id $(ViashDockerBuildArgs "$engine_id")
  elif [ "$setup_strategy" == "alwayscachedbuild" -o "$setup_strategy" == "cachedbuild" -o "$setup_strategy" == "cb" ]; then
    ViashDockerBuild $image_id $(ViashDockerBuildArgs "$engine_id")
  elif [[ "$setup_strategy" =~ ^ifneedbe ]]; then
    local save=$-; set +e
    ViashDockerLocalTagCheck $image_id
    local outCheck=$?
    [[ $save =~ e ]] && set -e
    if [ $outCheck -eq 0 ]; then
      ViashInfo "Image $image_id already exists"
    elif [ "$setup_strategy" == "ifneedbebuild" ]; then
      ViashDockerBuild $image_id --no-cache $(ViashDockerBuildArgs "$engine_id")
    elif [ "$setup_strategy" == "ifneedbecachedbuild" ]; then
      ViashDockerBuild $image_id $(ViashDockerBuildArgs "$engine_id")
    elif [ "$setup_strategy" == "ifneedbepull" ]; then
      ViashDockerPull $image_id
    elif [ "$setup_strategy" == "ifneedbepullelsebuild" ]; then
      ViashDockerPullElseBuild $image_id --no-cache $(ViashDockerBuildArgs "$engine_id")
    elif [ "$setup_strategy" == "ifneedbepullelsecachedbuild" ]; then
      ViashDockerPullElseBuild $image_id $(ViashDockerBuildArgs "$engine_id")
    else
      ViashError "Unrecognised Docker strategy: $setup_strategy"
      exit 1
    fi
  elif [ "$setup_strategy" == "push" -o "$setup_strategy" == "forcepush" -o "$setup_strategy" == "alwayspush" ]; then
    ViashDockerPush "$image_id"
  elif [ "$setup_strategy" == "pushifnotpresent" -o "$setup_strategy" == "gentlepush" -o "$setup_strategy" == "maybepush" ]; then
    local save=$-; set +e
    ViashDockerRemoteTagCheck $image_id
    local outCheck=$?
    [[ $save =~ e ]] && set -e
    if [ $outCheck -eq 0 ]; then
      ViashNotice "Container '$image_id' exists, doing nothing."
    else
      ViashNotice "Container '$image_id' does not yet exist."
      ViashDockerPush "$image_id"
    fi
  elif [ "$setup_strategy" == "donothing" -o "$setup_strategy" == "meh" ]; then
    ViashNotice "Skipping setup."
  else
    ViashError "Unrecognised Docker strategy: $setup_strategy"
    exit 1
  fi
}

# ViashDockerCheckCommands: Check whether a docker container has the required commands
#
# $1                  : image identifier with format `[registry/]image[:tag]`
# $@                  : commands to verify being present
# examples:
#   ViashDockerCheckCommands bash:4.0 bash ps foo
function ViashDockerCheckCommands {
  local image_id="$1"
  shift 1
  local commands="$@"
  local save=$-; set +e
  local missing # mark 'missing' as local in advance, otherwise the exit code of the command will be missing and always be '0'
  missing=$(docker run --rm --entrypoint=sh "$image_id" -c "for command in $commands; do command -v \$command >/dev/null 2>&1; if [ \$? -ne 0 ]; then echo \$command; exit 1; fi; done")
  local outCheck=$?
  [[ $save =~ e ]] && set -e
  if [ $outCheck -ne 0 ]; then
  	ViashError "Docker container '$image_id' does not contain command '$missing'."
  	exit 1
  fi
}

# ViashDockerBuild: build a docker image
# $1                               : image identifier with format `[registry/]image[:tag]`
# $...                             : additional arguments to pass to docker build
# $VIASH_META_TEMP_DIR             : temporary directory to store dockerfile & optional resources in
# $VIASH_META_NAME                 : name of the component
# $VIASH_META_RESOURCES_DIR        : directory containing the resources
# $VIASH_VERBOSITY                 : verbosity level
# exit code $?                     : whether or not the image was built successfully
function ViashDockerBuild {
  local image_id="$1"
  shift 1

  # create temporary directory to store dockerfile & optional resources in
  local tmpdir=$(mktemp -d "$VIASH_META_TEMP_DIR/dockerbuild-$VIASH_META_NAME-XXXXXX")
  local dockerfile="$tmpdir/Dockerfile"
  function clean_up {
    rm -rf "$tmpdir"
  }
  trap clean_up EXIT

  # store dockerfile and resources
  ViashDockerfile "$VIASH_ENGINE_ID" > "$dockerfile"

  # generate the build command
  local docker_build_cmd="docker build -t '$image_id' $@ '$VIASH_META_RESOURCES_DIR' -f '$dockerfile'"

  # build the container
  ViashNotice "Building container '$image_id' with Dockerfile"
  ViashInfo "$docker_build_cmd"
  local save=$-; set +e
  if [ $VIASH_VERBOSITY -ge $VIASH_LOGCODE_INFO ]; then
    eval $docker_build_cmd
  else
    eval $docker_build_cmd &> "$tmpdir/docker_build.log"
  fi

  # check exit code
  local out=$?
  [[ $save =~ e ]] && set -e
  if [ $out -ne 0 ]; then
    ViashError "Error occurred while building container '$image_id'"
    if [ $VIASH_VERBOSITY -lt $VIASH_LOGCODE_INFO ]; then
      ViashError "Transcript: --------------------------------"
      cat "$tmpdir/docker_build.log"
      ViashError "End of transcript --------------------------"
    fi
    exit 1
  fi
}

######## End of helper functions for setting up Docker images for viash ########

# ViashDockerFile: print the dockerfile to stdout
# $1    : engine identifier
# return : dockerfile required to run this component
# examples:
#   ViashDockerFile
function ViashDockerfile {
  local engine_id="$1"

  if [[ "$engine_id" == "docker" ]]; then
    cat << 'VIASHDOCKER'
FROM debian:stable-slim
ENTRYPOINT []
RUN apt-get update && \
  DEBIAN_FRONTEND=noninteractive apt-get install -y bcftools procps && \
  rm -rf /var/lib/apt/lists/*

RUN echo "bcftools: \"$(bcftools --version | grep 'bcftools' | sed -n 's/^bcftools //p')\"" > /var/software_versions.txt

LABEL org.opencontainers.image.authors="Theodoro Gasperin Terra Camargo"
LABEL org.opencontainers.image.description="Companion container for running component bcftools bcftools_annotate"
LABEL org.opencontainers.image.created="2024-11-14T12:50:48Z"
LABEL org.opencontainers.image.source="https://github.com/samtools/bcftools"
LABEL org.opencontainers.image.revision="402ad6a813c0d0b141bbf8f9a3aab9c270a74812"
LABEL org.opencontainers.image.version="add_sgdemux"

VIASHDOCKER
  fi
}

# ViashDockerBuildArgs: return the arguments to pass to docker build
# $1    : engine identifier
# return : arguments to pass to docker build
function ViashDockerBuildArgs {
  local engine_id="$1"

  if [[ "$engine_id" == "docker" ]]; then
    echo ""
  fi
}

# ViashAbsolutePath: generate absolute path from relative path
# borrowed from https://stackoverflow.com/a/21951256
# $1     : relative filename
# return : absolute path
# examples:
#   ViashAbsolutePath some_file.txt   # returns /path/to/some_file.txt
#   ViashAbsolutePath /foo/bar/..     # returns /foo
function ViashAbsolutePath {
  local thePath
  local parr
  local outp
  local len
  if [[ ! "$1" =~ ^/ ]]; then
    thePath="$PWD/$1"
  else
    thePath="$1"
  fi
  echo "$thePath" | (
    IFS=/
    read -a parr
    declare -a outp
    for i in "${parr[@]}"; do
      case "$i" in
      ''|.) continue ;;
      ..)
        len=${#outp[@]}
        if ((len==0)); then
          continue
        else
          unset outp[$((len-1))]
        fi
        ;;
      *)
        len=${#outp[@]}
        outp[$len]="$i"
      ;;
      esac
    done
    echo /"${outp[*]}"
  )
}
# ViashDockerAutodetectMount: auto configuring docker mounts from parameters
# $1                             : The parameter value
# returns                        : New parameter
# $VIASH_DIRECTORY_MOUNTS        : Added another parameter to be passed to docker
# $VIASH_DOCKER_AUTOMOUNT_PREFIX : The prefix to be used for the automounts
# examples:
#   ViashDockerAutodetectMount /path/to/bar      # returns '/viash_automount/path/to/bar'
#   ViashDockerAutodetectMountArg /path/to/bar   # returns '--volume="/path/to:/viash_automount/path/to"'
function ViashDockerAutodetectMount {
  local abs_path=$(ViashAbsolutePath "$1")
  local mount_source
  local base_name
  if [ -d "$abs_path" ]; then
    mount_source="$abs_path"
    base_name=""
  else
    mount_source=`dirname "$abs_path"`
    base_name=`basename "$abs_path"`
  fi
  local mount_target="$VIASH_DOCKER_AUTOMOUNT_PREFIX$mount_source"
  if [ -z "$base_name" ]; then
    echo "$mount_target"
  else
    echo "$mount_target/$base_name"
  fi
}
function ViashDockerAutodetectMountArg {
  local abs_path=$(ViashAbsolutePath "$1")
  local mount_source
  local base_name
  if [ -d "$abs_path" ]; then
    mount_source="$abs_path"
    base_name=""
  else
    mount_source=`dirname "$abs_path"`
    base_name=`basename "$abs_path"`
  fi
  local mount_target="$VIASH_DOCKER_AUTOMOUNT_PREFIX$mount_source"
  ViashDebug "ViashDockerAutodetectMountArg $1 -> $mount_source -> $mount_target"
  echo "--volume=\"$mount_source:$mount_target\""
}
function ViashDockerStripAutomount {
  local abs_path=$(ViashAbsolutePath "$1")
  echo "${abs_path#$VIASH_DOCKER_AUTOMOUNT_PREFIX}"
}
# initialise variables
VIASH_DIRECTORY_MOUNTS=()

# configure default docker automount prefix if it is unset
if [ -z "${VIASH_DOCKER_AUTOMOUNT_PREFIX+x}" ]; then
  VIASH_DOCKER_AUTOMOUNT_PREFIX="/viash_automount"
fi

# initialise docker variables
VIASH_DOCKER_RUN_ARGS=(-i --rm)

# initialise array
VIASH_POSITIONAL_ARGS=''

while [[ $# -gt 0 ]]; do
    case "$1" in
        -h|--help)
            ViashHelp
            exit
            ;;
        ---v|---verbose)
            let "VIASH_VERBOSITY=VIASH_VERBOSITY+1"
            shift 1
            ;;
        ---verbosity)
            VIASH_VERBOSITY="$2"
            shift 2
            ;;
        ---verbosity=*)
            VIASH_VERBOSITY="$(ViashRemoveFlags "$1")"
            shift 1
            ;;
        --version)
            echo "bcftools_annotate add_sgdemux"
            exit
            ;;
        --input)
            if [ -z "$VIASH_PAR_INPUT" ]; then
              VIASH_PAR_INPUT="$2"
            else
              VIASH_PAR_INPUT="$VIASH_PAR_INPUT;""$2"
            fi
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --input. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --input=*)
            if [ -z "$VIASH_PAR_INPUT" ]; then
              VIASH_PAR_INPUT=$(ViashRemoveFlags "$1")
            else
              VIASH_PAR_INPUT="$VIASH_PAR_INPUT;"$(ViashRemoveFlags "$1")
            fi
            shift 1
            ;;
        -i)
            if [ -z "$VIASH_PAR_INPUT" ]; then
              VIASH_PAR_INPUT="$2"
            else
              VIASH_PAR_INPUT="$VIASH_PAR_INPUT;""$2"
            fi
            [ $# -lt 2 ] && ViashError Not enough arguments passed to -i. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --output)
            [ -n "$VIASH_PAR_OUTPUT" ] && ViashError Bad arguments for option \'--output\': \'$VIASH_PAR_OUTPUT\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_OUTPUT="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --output. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --output=*)
            [ -n "$VIASH_PAR_OUTPUT" ] && ViashError Bad arguments for option \'--output=*\': \'$VIASH_PAR_OUTPUT\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_OUTPUT=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        -o)
            [ -n "$VIASH_PAR_OUTPUT" ] && ViashError Bad arguments for option \'-o\': \'$VIASH_PAR_OUTPUT\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_OUTPUT="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to -o. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --annotations)
            [ -n "$VIASH_PAR_ANNOTATIONS" ] && ViashError Bad arguments for option \'--annotations\': \'$VIASH_PAR_ANNOTATIONS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_ANNOTATIONS="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --annotations. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --annotations=*)
            [ -n "$VIASH_PAR_ANNOTATIONS" ] && ViashError Bad arguments for option \'--annotations=*\': \'$VIASH_PAR_ANNOTATIONS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_ANNOTATIONS=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --a)
            [ -n "$VIASH_PAR_ANNOTATIONS" ] && ViashError Bad arguments for option \'--a\': \'$VIASH_PAR_ANNOTATIONS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_ANNOTATIONS="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --a. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --columns)
            [ -n "$VIASH_PAR_COLUMNS" ] && ViashError Bad arguments for option \'--columns\': \'$VIASH_PAR_COLUMNS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_COLUMNS="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --columns. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --columns=*)
            [ -n "$VIASH_PAR_COLUMNS" ] && ViashError Bad arguments for option \'--columns=*\': \'$VIASH_PAR_COLUMNS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_COLUMNS=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --c)
            [ -n "$VIASH_PAR_COLUMNS" ] && ViashError Bad arguments for option \'--c\': \'$VIASH_PAR_COLUMNS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_COLUMNS="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --c. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --columns_file)
            [ -n "$VIASH_PAR_COLUMNS_FILE" ] && ViashError Bad arguments for option \'--columns_file\': \'$VIASH_PAR_COLUMNS_FILE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_COLUMNS_FILE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --columns_file. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --columns_file=*)
            [ -n "$VIASH_PAR_COLUMNS_FILE" ] && ViashError Bad arguments for option \'--columns_file=*\': \'$VIASH_PAR_COLUMNS_FILE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_COLUMNS_FILE=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --C)
            [ -n "$VIASH_PAR_COLUMNS_FILE" ] && ViashError Bad arguments for option \'--C\': \'$VIASH_PAR_COLUMNS_FILE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_COLUMNS_FILE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --C. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --exclude)
            [ -n "$VIASH_PAR_EXCLUDE" ] && ViashError Bad arguments for option \'--exclude\': \'$VIASH_PAR_EXCLUDE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_EXCLUDE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --exclude. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --exclude=*)
            [ -n "$VIASH_PAR_EXCLUDE" ] && ViashError Bad arguments for option \'--exclude=*\': \'$VIASH_PAR_EXCLUDE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_EXCLUDE=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --e)
            [ -n "$VIASH_PAR_EXCLUDE" ] && ViashError Bad arguments for option \'--e\': \'$VIASH_PAR_EXCLUDE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_EXCLUDE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --e. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --force)
            [ -n "$VIASH_PAR_FORCE" ] && ViashError Bad arguments for option \'--force\': \'$VIASH_PAR_FORCE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_FORCE=true
            shift 1
            ;;
        --header_line)
            [ -n "$VIASH_PAR_HEADER_LINE" ] && ViashError Bad arguments for option \'--header_line\': \'$VIASH_PAR_HEADER_LINE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_HEADER_LINE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --header_line. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --header_line=*)
            [ -n "$VIASH_PAR_HEADER_LINE" ] && ViashError Bad arguments for option \'--header_line=*\': \'$VIASH_PAR_HEADER_LINE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_HEADER_LINE=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --H)
            [ -n "$VIASH_PAR_HEADER_LINE" ] && ViashError Bad arguments for option \'--H\': \'$VIASH_PAR_HEADER_LINE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_HEADER_LINE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --H. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --header_lines)
            [ -n "$VIASH_PAR_HEADER_LINES" ] && ViashError Bad arguments for option \'--header_lines\': \'$VIASH_PAR_HEADER_LINES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_HEADER_LINES="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --header_lines. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --header_lines=*)
            [ -n "$VIASH_PAR_HEADER_LINES" ] && ViashError Bad arguments for option \'--header_lines=*\': \'$VIASH_PAR_HEADER_LINES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_HEADER_LINES=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --h)
            [ -n "$VIASH_PAR_HEADER_LINES" ] && ViashError Bad arguments for option \'--h\': \'$VIASH_PAR_HEADER_LINES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_HEADER_LINES="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --h. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --set_id)
            [ -n "$VIASH_PAR_SET_ID" ] && ViashError Bad arguments for option \'--set_id\': \'$VIASH_PAR_SET_ID\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_SET_ID="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --set_id. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --set_id=*)
            [ -n "$VIASH_PAR_SET_ID" ] && ViashError Bad arguments for option \'--set_id=*\': \'$VIASH_PAR_SET_ID\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_SET_ID=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --I)
            [ -n "$VIASH_PAR_SET_ID" ] && ViashError Bad arguments for option \'--I\': \'$VIASH_PAR_SET_ID\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_SET_ID="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --I. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --include)
            [ -n "$VIASH_PAR_INCLUDE" ] && ViashError Bad arguments for option \'--include\': \'$VIASH_PAR_INCLUDE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_INCLUDE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --include. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --include=*)
            [ -n "$VIASH_PAR_INCLUDE" ] && ViashError Bad arguments for option \'--include=*\': \'$VIASH_PAR_INCLUDE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_INCLUDE=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --keep_sites)
            [ -n "$VIASH_PAR_KEEP_SITES" ] && ViashError Bad arguments for option \'--keep_sites\': \'$VIASH_PAR_KEEP_SITES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_KEEP_SITES=true
            shift 1
            ;;
        --k)
            [ -n "$VIASH_PAR_KEEP_SITES" ] && ViashError Bad arguments for option \'--k\': \'$VIASH_PAR_KEEP_SITES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_KEEP_SITES=true
            shift 1
            ;;
        --merge_logic)
            [ -n "$VIASH_PAR_MERGE_LOGIC" ] && ViashError Bad arguments for option \'--merge_logic\': \'$VIASH_PAR_MERGE_LOGIC\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_MERGE_LOGIC="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --merge_logic. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --merge_logic=*)
            [ -n "$VIASH_PAR_MERGE_LOGIC" ] && ViashError Bad arguments for option \'--merge_logic=*\': \'$VIASH_PAR_MERGE_LOGIC\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_MERGE_LOGIC=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --l)
            [ -n "$VIASH_PAR_MERGE_LOGIC" ] && ViashError Bad arguments for option \'--l\': \'$VIASH_PAR_MERGE_LOGIC\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_MERGE_LOGIC="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --l. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --mark_sites)
            [ -n "$VIASH_PAR_MARK_SITES" ] && ViashError Bad arguments for option \'--mark_sites\': \'$VIASH_PAR_MARK_SITES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_MARK_SITES="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --mark_sites. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --mark_sites=*)
            [ -n "$VIASH_PAR_MARK_SITES" ] && ViashError Bad arguments for option \'--mark_sites=*\': \'$VIASH_PAR_MARK_SITES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_MARK_SITES=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --m)
            [ -n "$VIASH_PAR_MARK_SITES" ] && ViashError Bad arguments for option \'--m\': \'$VIASH_PAR_MARK_SITES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_MARK_SITES="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --m. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --min_overlap)
            [ -n "$VIASH_PAR_MIN_OVERLAP" ] && ViashError Bad arguments for option \'--min_overlap\': \'$VIASH_PAR_MIN_OVERLAP\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_MIN_OVERLAP="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --min_overlap. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --min_overlap=*)
            [ -n "$VIASH_PAR_MIN_OVERLAP" ] && ViashError Bad arguments for option \'--min_overlap=*\': \'$VIASH_PAR_MIN_OVERLAP\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_MIN_OVERLAP=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --no_version)
            [ -n "$VIASH_PAR_NO_VERSION" ] && ViashError Bad arguments for option \'--no_version\': \'$VIASH_PAR_NO_VERSION\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_NO_VERSION=true
            shift 1
            ;;
        --output_type)
            [ -n "$VIASH_PAR_OUTPUT_TYPE" ] && ViashError Bad arguments for option \'--output_type\': \'$VIASH_PAR_OUTPUT_TYPE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_OUTPUT_TYPE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --output_type. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --output_type=*)
            [ -n "$VIASH_PAR_OUTPUT_TYPE" ] && ViashError Bad arguments for option \'--output_type=*\': \'$VIASH_PAR_OUTPUT_TYPE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_OUTPUT_TYPE=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --O)
            [ -n "$VIASH_PAR_OUTPUT_TYPE" ] && ViashError Bad arguments for option \'--O\': \'$VIASH_PAR_OUTPUT_TYPE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_OUTPUT_TYPE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --O. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --pair_logic)
            [ -n "$VIASH_PAR_PAIR_LOGIC" ] && ViashError Bad arguments for option \'--pair_logic\': \'$VIASH_PAR_PAIR_LOGIC\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_PAIR_LOGIC="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --pair_logic. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --pair_logic=*)
            [ -n "$VIASH_PAR_PAIR_LOGIC" ] && ViashError Bad arguments for option \'--pair_logic=*\': \'$VIASH_PAR_PAIR_LOGIC\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_PAIR_LOGIC=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --regions)
            [ -n "$VIASH_PAR_REGIONS" ] && ViashError Bad arguments for option \'--regions\': \'$VIASH_PAR_REGIONS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REGIONS="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --regions. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --regions=*)
            [ -n "$VIASH_PAR_REGIONS" ] && ViashError Bad arguments for option \'--regions=*\': \'$VIASH_PAR_REGIONS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REGIONS=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --r)
            [ -n "$VIASH_PAR_REGIONS" ] && ViashError Bad arguments for option \'--r\': \'$VIASH_PAR_REGIONS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REGIONS="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --r. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --regions_file)
            [ -n "$VIASH_PAR_REGIONS_FILE" ] && ViashError Bad arguments for option \'--regions_file\': \'$VIASH_PAR_REGIONS_FILE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REGIONS_FILE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --regions_file. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --regions_file=*)
            [ -n "$VIASH_PAR_REGIONS_FILE" ] && ViashError Bad arguments for option \'--regions_file=*\': \'$VIASH_PAR_REGIONS_FILE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REGIONS_FILE=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --R)
            [ -n "$VIASH_PAR_REGIONS_FILE" ] && ViashError Bad arguments for option \'--R\': \'$VIASH_PAR_REGIONS_FILE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REGIONS_FILE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --R. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --regions_overlap)
            [ -n "$VIASH_PAR_REGIONS_OVERLAP" ] && ViashError Bad arguments for option \'--regions_overlap\': \'$VIASH_PAR_REGIONS_OVERLAP\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REGIONS_OVERLAP="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --regions_overlap. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --regions_overlap=*)
            [ -n "$VIASH_PAR_REGIONS_OVERLAP" ] && ViashError Bad arguments for option \'--regions_overlap=*\': \'$VIASH_PAR_REGIONS_OVERLAP\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REGIONS_OVERLAP=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --rename_annotations)
            [ -n "$VIASH_PAR_RENAME_ANNOTATIONS" ] && ViashError Bad arguments for option \'--rename_annotations\': \'$VIASH_PAR_RENAME_ANNOTATIONS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_RENAME_ANNOTATIONS="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --rename_annotations. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --rename_annotations=*)
            [ -n "$VIASH_PAR_RENAME_ANNOTATIONS" ] && ViashError Bad arguments for option \'--rename_annotations=*\': \'$VIASH_PAR_RENAME_ANNOTATIONS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_RENAME_ANNOTATIONS=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --rename_chromosomes)
            [ -n "$VIASH_PAR_RENAME_CHROMOSOMES" ] && ViashError Bad arguments for option \'--rename_chromosomes\': \'$VIASH_PAR_RENAME_CHROMOSOMES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_RENAME_CHROMOSOMES="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --rename_chromosomes. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --rename_chromosomes=*)
            [ -n "$VIASH_PAR_RENAME_CHROMOSOMES" ] && ViashError Bad arguments for option \'--rename_chromosomes=*\': \'$VIASH_PAR_RENAME_CHROMOSOMES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_RENAME_CHROMOSOMES=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --samples)
            [ -n "$VIASH_PAR_SAMPLES" ] && ViashError Bad arguments for option \'--samples\': \'$VIASH_PAR_SAMPLES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_SAMPLES="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --samples. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --samples=*)
            [ -n "$VIASH_PAR_SAMPLES" ] && ViashError Bad arguments for option \'--samples=*\': \'$VIASH_PAR_SAMPLES\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_SAMPLES=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --samples_file)
            [ -n "$VIASH_PAR_SAMPLES_FILE" ] && ViashError Bad arguments for option \'--samples_file\': \'$VIASH_PAR_SAMPLES_FILE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_SAMPLES_FILE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --samples_file. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --samples_file=*)
            [ -n "$VIASH_PAR_SAMPLES_FILE" ] && ViashError Bad arguments for option \'--samples_file=*\': \'$VIASH_PAR_SAMPLES_FILE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_SAMPLES_FILE=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --single_overlaps)
            [ -n "$VIASH_PAR_SINGLE_OVERLAPS" ] && ViashError Bad arguments for option \'--single_overlaps\': \'$VIASH_PAR_SINGLE_OVERLAPS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_SINGLE_OVERLAPS=true
            shift 1
            ;;
        --remove)
            [ -n "$VIASH_PAR_REMOVE" ] && ViashError Bad arguments for option \'--remove\': \'$VIASH_PAR_REMOVE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REMOVE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --remove. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        --remove=*)
            [ -n "$VIASH_PAR_REMOVE" ] && ViashError Bad arguments for option \'--remove=*\': \'$VIASH_PAR_REMOVE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REMOVE=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        --x)
            [ -n "$VIASH_PAR_REMOVE" ] && ViashError Bad arguments for option \'--x\': \'$VIASH_PAR_REMOVE\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_PAR_REMOVE="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to --x. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        ---engine)
            VIASH_ENGINE_ID="$2"
            shift 2
            ;;
        ---engine=*)
            VIASH_ENGINE_ID="$(ViashRemoveFlags "$1")"
            shift 1
            ;;
        ---setup)
            VIASH_MODE='setup'
            VIASH_SETUP_STRATEGY="$2"
            shift 2
            ;;
        ---setup=*)
            VIASH_MODE='setup'
            VIASH_SETUP_STRATEGY="$(ViashRemoveFlags "$1")"
            shift 1
            ;;
        ---dockerfile)
            VIASH_MODE='dockerfile'
            shift 1
            ;;
        ---docker_run_args)
            VIASH_DOCKER_RUN_ARGS+=("$2")
            shift 2
            ;;
        ---docker_run_args=*)
            VIASH_DOCKER_RUN_ARGS+=("$(ViashRemoveFlags "$1")")
            shift 1
            ;;
        ---docker_image_id)
            VIASH_MODE='docker_image_id'
            shift 1
            ;;
        ---debug)
            VIASH_MODE='debug'
            shift 1
            ;;
        ---cpus)
            [ -n "$VIASH_META_CPUS" ] && ViashError Bad arguments for option \'---cpus\': \'$VIASH_META_CPUS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_META_CPUS="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to ---cpus. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        ---cpus=*)
            [ -n "$VIASH_META_CPUS" ] && ViashError Bad arguments for option \'---cpus=*\': \'$VIASH_META_CPUS\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_META_CPUS=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        ---memory)
            [ -n "$VIASH_META_MEMORY" ] && ViashError Bad arguments for option \'---memory\': \'$VIASH_META_MEMORY\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_META_MEMORY="$2"
            [ $# -lt 2 ] && ViashError Not enough arguments passed to ---memory. Use "--help" to get more information on the parameters. && exit 1
            shift 2
            ;;
        ---memory=*)
            [ -n "$VIASH_META_MEMORY" ] && ViashError Bad arguments for option \'---memory=*\': \'$VIASH_META_MEMORY\' \& \'$2\' - you should provide exactly one argument for this option. && exit 1
            VIASH_META_MEMORY=$(ViashRemoveFlags "$1")
            shift 1
            ;;
        *)  # positional arg or unknown option
            # since the positional args will be eval'd, can we always quote, instead of using ViashQuote
            VIASH_POSITIONAL_ARGS="$VIASH_POSITIONAL_ARGS '$1'"
            [[ $1 == -* ]] && ViashWarning $1 looks like a parameter but is not a defined parameter and will instead be treated as a positional argument. Use "--help" to get more information on the parameters.
            shift # past argument
            ;;
    esac
done

# parse positional parameters
eval set -- $VIASH_POSITIONAL_ARGS


if   [ "$VIASH_ENGINE_ID" == "native" ]  ; then
  VIASH_ENGINE_TYPE='native'
elif   [ "$VIASH_ENGINE_ID" == "docker" ]  ; then
  VIASH_ENGINE_TYPE='docker'
else
  ViashError "Engine '$VIASH_ENGINE_ID' is not recognized. Options are: docker, native."
  exit 1
fi

if [[ "$VIASH_ENGINE_TYPE" == "docker" ]]; then
  # check if docker is installed properly
  ViashDockerInstallationCheck

  # determine docker image id
  if [[ "$VIASH_ENGINE_ID" == 'docker' ]]; then
    VIASH_DOCKER_IMAGE_ID='images.viash-hub.com/vsh/biobox/bcftools/bcftools_annotate:add_sgdemux'
  fi

  # print dockerfile
  if [ "$VIASH_MODE" == "dockerfile" ]; then
    ViashDockerfile "$VIASH_ENGINE_ID"
    exit 0

  elif [ "$VIASH_MODE" == "docker_image_id" ]; then
    echo "$VIASH_DOCKER_IMAGE_ID"
    exit 0
  
  # enter docker container
  elif [[ "$VIASH_MODE" == "debug" ]]; then
    VIASH_CMD="docker run --entrypoint=bash ${VIASH_DOCKER_RUN_ARGS[@]} -v '$(pwd)':/pwd --workdir /pwd -t $VIASH_DOCKER_IMAGE_ID"
    ViashNotice "+ $VIASH_CMD"
    eval $VIASH_CMD
    exit 

  # build docker image
  elif [ "$VIASH_MODE" == "setup" ]; then
    ViashDockerSetup "$VIASH_DOCKER_IMAGE_ID" "$VIASH_SETUP_STRATEGY"
    ViashDockerCheckCommands "$VIASH_DOCKER_IMAGE_ID" 'ps' 'bash'
    exit 0
  fi

  # check if docker image exists
  ViashDockerSetup "$VIASH_DOCKER_IMAGE_ID" ifneedbepullelsecachedbuild
  ViashDockerCheckCommands "$VIASH_DOCKER_IMAGE_ID" 'ps' 'bash'
fi

# setting computational defaults

# helper function for parsing memory strings
function ViashMemoryAsBytes {
  local memory=`echo "$1" | tr '[:upper:]' '[:lower:]' | tr -d '[:space:]'`
  local memory_regex='^([0-9]+)([kmgtp]i?b?|b)$'
  if [[ $memory =~ $memory_regex ]]; then
    local number=${memory/[^0-9]*/}
    local symbol=${memory/*[0-9]/}
    
    case $symbol in
      b)      memory_b=$number ;;
      kb|k)   memory_b=$(( $number * 1000 )) ;;
      mb|m)   memory_b=$(( $number * 1000 * 1000 )) ;;
      gb|g)   memory_b=$(( $number * 1000 * 1000 * 1000 )) ;;
      tb|t)   memory_b=$(( $number * 1000 * 1000 * 1000 * 1000 )) ;;
      pb|p)   memory_b=$(( $number * 1000 * 1000 * 1000 * 1000 * 1000 )) ;;
      kib|ki)   memory_b=$(( $number * 1024 )) ;;
      mib|mi)   memory_b=$(( $number * 1024 * 1024 )) ;;
      gib|gi)   memory_b=$(( $number * 1024 * 1024 * 1024 )) ;;
      tib|ti)   memory_b=$(( $number * 1024 * 1024 * 1024 * 1024 )) ;;
      pib|pi)   memory_b=$(( $number * 1024 * 1024 * 1024 * 1024 * 1024 )) ;;
    esac
    echo "$memory_b"
  fi
}
# compute memory in different units
if [ ! -z ${VIASH_META_MEMORY+x} ]; then
  VIASH_META_MEMORY_B=`ViashMemoryAsBytes $VIASH_META_MEMORY`
  # do not define other variables if memory_b is an empty string
  if [ ! -z "$VIASH_META_MEMORY_B" ]; then
    VIASH_META_MEMORY_KB=$(( ($VIASH_META_MEMORY_B+999) / 1000 ))
    VIASH_META_MEMORY_MB=$(( ($VIASH_META_MEMORY_KB+999) / 1000 ))
    VIASH_META_MEMORY_GB=$(( ($VIASH_META_MEMORY_MB+999) / 1000 ))
    VIASH_META_MEMORY_TB=$(( ($VIASH_META_MEMORY_GB+999) / 1000 ))
    VIASH_META_MEMORY_PB=$(( ($VIASH_META_MEMORY_TB+999) / 1000 ))
    VIASH_META_MEMORY_KIB=$(( ($VIASH_META_MEMORY_B+1023) / 1024 ))
    VIASH_META_MEMORY_MIB=$(( ($VIASH_META_MEMORY_KIB+1023) / 1024 ))
    VIASH_META_MEMORY_GIB=$(( ($VIASH_META_MEMORY_MIB+1023) / 1024 ))
    VIASH_META_MEMORY_TIB=$(( ($VIASH_META_MEMORY_GIB+1023) / 1024 ))
    VIASH_META_MEMORY_PIB=$(( ($VIASH_META_MEMORY_TIB+1023) / 1024 ))
  else
    # unset memory if string is empty
    unset $VIASH_META_MEMORY_B
  fi
fi
# unset nproc if string is empty
if [ -z "$VIASH_META_CPUS" ]; then
  unset $VIASH_META_CPUS
fi


# check whether required parameters exist
if [ -z ${VIASH_PAR_INPUT+x} ]; then
  ViashError '--input' is a required argument. Use "--help" to get more information on the parameters.
  exit 1
fi
if [ -z ${VIASH_PAR_OUTPUT+x} ]; then
  ViashError '--output' is a required argument. Use "--help" to get more information on the parameters.
  exit 1
fi
if [ -z ${VIASH_META_NAME+x} ]; then
  ViashError 'name' is a required argument. Use "--help" to get more information on the parameters.
  exit 1
fi
if [ -z ${VIASH_META_FUNCTIONALITY_NAME+x} ]; then
  ViashError 'functionality_name' is a required argument. Use "--help" to get more information on the parameters.
  exit 1
fi
if [ -z ${VIASH_META_RESOURCES_DIR+x} ]; then
  ViashError 'resources_dir' is a required argument. Use "--help" to get more information on the parameters.
  exit 1
fi
if [ -z ${VIASH_META_EXECUTABLE+x} ]; then
  ViashError 'executable' is a required argument. Use "--help" to get more information on the parameters.
  exit 1
fi
if [ -z ${VIASH_META_CONFIG+x} ]; then
  ViashError 'config' is a required argument. Use "--help" to get more information on the parameters.
  exit 1
fi
if [ -z ${VIASH_META_TEMP_DIR+x} ]; then
  ViashError 'temp_dir' is a required argument. Use "--help" to get more information on the parameters.
  exit 1
fi

# filling in defaults
if [ -z ${VIASH_PAR_FORCE+x} ]; then
  VIASH_PAR_FORCE="false"
fi
if [ -z ${VIASH_PAR_KEEP_SITES+x} ]; then
  VIASH_PAR_KEEP_SITES="false"
fi
if [ -z ${VIASH_PAR_NO_VERSION+x} ]; then
  VIASH_PAR_NO_VERSION="false"
fi
if [ -z ${VIASH_PAR_SINGLE_OVERLAPS+x} ]; then
  VIASH_PAR_SINGLE_OVERLAPS="false"
fi

# check whether required files exist
if [ ! -z "$VIASH_PAR_INPUT" ]; then
  IFS=';'
  set -f
  for file in $VIASH_PAR_INPUT; do
    unset IFS
    if [ ! -e "$file" ]; then
      ViashError "Input file '$file' does not exist."
      exit 1
    fi
  done
  set +f
fi
if [ ! -z "$VIASH_PAR_ANNOTATIONS" ] && [ ! -e "$VIASH_PAR_ANNOTATIONS" ]; then
  ViashError "Input file '$VIASH_PAR_ANNOTATIONS' does not exist."
  exit 1
fi
if [ ! -z "$VIASH_PAR_COLUMNS_FILE" ] && [ ! -e "$VIASH_PAR_COLUMNS_FILE" ]; then
  ViashError "Input file '$VIASH_PAR_COLUMNS_FILE' does not exist."
  exit 1
fi
if [ ! -z "$VIASH_PAR_HEADER_LINES" ] && [ ! -e "$VIASH_PAR_HEADER_LINES" ]; then
  ViashError "Input file '$VIASH_PAR_HEADER_LINES' does not exist."
  exit 1
fi
if [ ! -z "$VIASH_PAR_REGIONS_FILE" ] && [ ! -e "$VIASH_PAR_REGIONS_FILE" ]; then
  ViashError "Input file '$VIASH_PAR_REGIONS_FILE' does not exist."
  exit 1
fi
if [ ! -z "$VIASH_PAR_RENAME_ANNOTATIONS" ] && [ ! -e "$VIASH_PAR_RENAME_ANNOTATIONS" ]; then
  ViashError "Input file '$VIASH_PAR_RENAME_ANNOTATIONS' does not exist."
  exit 1
fi
if [ ! -z "$VIASH_PAR_RENAME_CHROMOSOMES" ] && [ ! -e "$VIASH_PAR_RENAME_CHROMOSOMES" ]; then
  ViashError "Input file '$VIASH_PAR_RENAME_CHROMOSOMES' does not exist."
  exit 1
fi
if [ ! -z "$VIASH_PAR_SAMPLES_FILE" ] && [ ! -e "$VIASH_PAR_SAMPLES_FILE" ]; then
  ViashError "Input file '$VIASH_PAR_SAMPLES_FILE' does not exist."
  exit 1
fi

# check whether parameters values are of the right type
if [[ -n "$VIASH_PAR_FORCE" ]]; then
  if ! [[ "$VIASH_PAR_FORCE" =~ ^(true|True|TRUE|false|False|FALSE|yes|Yes|YES|no|No|NO)$ ]]; then
    ViashError '--force' has to be a boolean_true. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_PAR_KEEP_SITES" ]]; then
  if ! [[ "$VIASH_PAR_KEEP_SITES" =~ ^(true|True|TRUE|false|False|FALSE|yes|Yes|YES|no|No|NO)$ ]]; then
    ViashError '--keep_sites' has to be a boolean_true. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_PAR_NO_VERSION" ]]; then
  if ! [[ "$VIASH_PAR_NO_VERSION" =~ ^(true|True|TRUE|false|False|FALSE|yes|Yes|YES|no|No|NO)$ ]]; then
    ViashError '--no_version' has to be a boolean_true. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_PAR_SINGLE_OVERLAPS" ]]; then
  if ! [[ "$VIASH_PAR_SINGLE_OVERLAPS" =~ ^(true|True|TRUE|false|False|FALSE|yes|Yes|YES|no|No|NO)$ ]]; then
    ViashError '--single_overlaps' has to be a boolean_true. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_CPUS" ]]; then
  if ! [[ "$VIASH_META_CPUS" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'cpus' has to be an integer. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_B" ]]; then
  if ! [[ "$VIASH_META_MEMORY_B" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_b' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_KB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_KB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_kb' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_MB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_MB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_mb' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_GB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_GB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_gb' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_TB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_TB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_tb' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_PB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_PB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_pb' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_KIB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_KIB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_kib' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_MIB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_MIB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_mib' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_GIB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_GIB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_gib' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_TIB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_TIB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_tib' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi
if [[ -n "$VIASH_META_MEMORY_PIB" ]]; then
  if ! [[ "$VIASH_META_MEMORY_PIB" =~ ^[-+]?[0-9]+$ ]]; then
    ViashError 'memory_pib' has to be a long. Use "--help" to get more information on the parameters.
    exit 1
  fi
fi

# check whether value is belongs to a set of choices
if [ ! -z "$VIASH_PAR_OUTPUT_TYPE" ]; then
  VIASH_PAR_OUTPUT_TYPE_CHOICES=("u;z;b;v")
  IFS=';'
  set -f
  if ! [[ ";${VIASH_PAR_OUTPUT_TYPE_CHOICES[*]};" =~ ";$VIASH_PAR_OUTPUT_TYPE;" ]]; then
    ViashError '--output_type' specified value of \'$VIASH_PAR_OUTPUT_TYPE\' is not in the list of allowed values. Use "--help" to get more information on the parameters.
    exit 1
  fi
  set +f
  unset IFS
fi

if [ ! -z "$VIASH_PAR_PAIR_LOGIC" ]; then
  VIASH_PAR_PAIR_LOGIC_CHOICES=("snps;indels;both;all;some;exact")
  IFS=';'
  set -f
  if ! [[ ";${VIASH_PAR_PAIR_LOGIC_CHOICES[*]};" =~ ";$VIASH_PAR_PAIR_LOGIC;" ]]; then
    ViashError '--pair_logic' specified value of \'$VIASH_PAR_PAIR_LOGIC\' is not in the list of allowed values. Use "--help" to get more information on the parameters.
    exit 1
  fi
  set +f
  unset IFS
fi

if [ ! -z "$VIASH_PAR_REGIONS_OVERLAP" ]; then
  VIASH_PAR_REGIONS_OVERLAP_CHOICES=("pos;record;variant;0;1;2")
  IFS=';'
  set -f
  if ! [[ ";${VIASH_PAR_REGIONS_OVERLAP_CHOICES[*]};" =~ ";$VIASH_PAR_REGIONS_OVERLAP;" ]]; then
    ViashError '--regions_overlap' specified value of \'$VIASH_PAR_REGIONS_OVERLAP\' is not in the list of allowed values. Use "--help" to get more information on the parameters.
    exit 1
  fi
  set +f
  unset IFS
fi

# create parent directories of output files, if so desired
if [ ! -z "$VIASH_PAR_OUTPUT" ] && [ ! -d "$(dirname "$VIASH_PAR_OUTPUT")" ]; then
  mkdir -p "$(dirname "$VIASH_PAR_OUTPUT")"
fi

if  [ "$VIASH_ENGINE_ID" == "native" ]  ; then
  if [ "$VIASH_MODE" == "run" ]; then
    VIASH_CMD="bash"
  else
    ViashError "Engine '$VIASH_ENGINE_ID' does not support mode '$VIASH_MODE'."
    exit 1
  fi
fi

if [[ "$VIASH_ENGINE_TYPE" == "docker" ]]; then
  # detect volumes from file arguments
  VIASH_CHOWN_VARS=()
if [ ! -z "$VIASH_PAR_INPUT" ]; then
  VIASH_TEST_INPUT=()
  IFS=';'
  for var in $VIASH_PAR_INPUT; do
    unset IFS
    VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$var")" )
    var=$(ViashDockerAutodetectMount "$var")
    VIASH_TEST_INPUT+=( "$var" )
  done
  VIASH_PAR_INPUT=$(IFS=';' ; echo "${VIASH_TEST_INPUT[*]}")
fi
if [ ! -z "$VIASH_PAR_OUTPUT" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_PAR_OUTPUT")" )
  VIASH_PAR_OUTPUT=$(ViashDockerAutodetectMount "$VIASH_PAR_OUTPUT")
  VIASH_CHOWN_VARS+=( "$VIASH_PAR_OUTPUT" )
fi
if [ ! -z "$VIASH_PAR_ANNOTATIONS" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_PAR_ANNOTATIONS")" )
  VIASH_PAR_ANNOTATIONS=$(ViashDockerAutodetectMount "$VIASH_PAR_ANNOTATIONS")
fi
if [ ! -z "$VIASH_PAR_COLUMNS_FILE" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_PAR_COLUMNS_FILE")" )
  VIASH_PAR_COLUMNS_FILE=$(ViashDockerAutodetectMount "$VIASH_PAR_COLUMNS_FILE")
fi
if [ ! -z "$VIASH_PAR_HEADER_LINES" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_PAR_HEADER_LINES")" )
  VIASH_PAR_HEADER_LINES=$(ViashDockerAutodetectMount "$VIASH_PAR_HEADER_LINES")
fi
if [ ! -z "$VIASH_PAR_REGIONS_FILE" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_PAR_REGIONS_FILE")" )
  VIASH_PAR_REGIONS_FILE=$(ViashDockerAutodetectMount "$VIASH_PAR_REGIONS_FILE")
fi
if [ ! -z "$VIASH_PAR_RENAME_ANNOTATIONS" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_PAR_RENAME_ANNOTATIONS")" )
  VIASH_PAR_RENAME_ANNOTATIONS=$(ViashDockerAutodetectMount "$VIASH_PAR_RENAME_ANNOTATIONS")
fi
if [ ! -z "$VIASH_PAR_RENAME_CHROMOSOMES" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_PAR_RENAME_CHROMOSOMES")" )
  VIASH_PAR_RENAME_CHROMOSOMES=$(ViashDockerAutodetectMount "$VIASH_PAR_RENAME_CHROMOSOMES")
fi
if [ ! -z "$VIASH_PAR_SAMPLES_FILE" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_PAR_SAMPLES_FILE")" )
  VIASH_PAR_SAMPLES_FILE=$(ViashDockerAutodetectMount "$VIASH_PAR_SAMPLES_FILE")
fi
if [ ! -z "$VIASH_META_RESOURCES_DIR" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_META_RESOURCES_DIR")" )
  VIASH_META_RESOURCES_DIR=$(ViashDockerAutodetectMount "$VIASH_META_RESOURCES_DIR")
fi
if [ ! -z "$VIASH_META_EXECUTABLE" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_META_EXECUTABLE")" )
  VIASH_META_EXECUTABLE=$(ViashDockerAutodetectMount "$VIASH_META_EXECUTABLE")
fi
if [ ! -z "$VIASH_META_CONFIG" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_META_CONFIG")" )
  VIASH_META_CONFIG=$(ViashDockerAutodetectMount "$VIASH_META_CONFIG")
fi
if [ ! -z "$VIASH_META_TEMP_DIR" ]; then
  VIASH_DIRECTORY_MOUNTS+=( "$(ViashDockerAutodetectMountArg "$VIASH_META_TEMP_DIR")" )
  VIASH_META_TEMP_DIR=$(ViashDockerAutodetectMount "$VIASH_META_TEMP_DIR")
fi
  
  # get unique mounts
  VIASH_UNIQUE_MOUNTS=($(for val in "${VIASH_DIRECTORY_MOUNTS[@]}"; do echo "$val"; done | sort -u))
fi

if [[ "$VIASH_ENGINE_TYPE" == "docker" ]]; then
  # change file ownership
  function ViashPerformChown {
    if (( ${#VIASH_CHOWN_VARS[@]} )); then
      set +e
      VIASH_CMD="docker run --entrypoint=bash --rm ${VIASH_UNIQUE_MOUNTS[@]} $VIASH_DOCKER_IMAGE_ID -c 'chown $(id -u):$(id -g) --silent --recursive ${VIASH_CHOWN_VARS[@]}'"
      ViashDebug "+ $VIASH_CMD"
      eval $VIASH_CMD
      set -e
    fi
  }
  trap ViashPerformChown EXIT
fi

if [[ "$VIASH_ENGINE_TYPE" == "docker" ]]; then
  # helper function for filling in extra docker args
  if [ ! -z "$VIASH_META_MEMORY_B" ]; then
    VIASH_DOCKER_RUN_ARGS+=("--memory=${VIASH_META_MEMORY_B}")
  fi
  if [ ! -z "$VIASH_META_CPUS" ]; then
    VIASH_DOCKER_RUN_ARGS+=("--cpus=${VIASH_META_CPUS}")
  fi
fi

if [[ "$VIASH_ENGINE_TYPE" == "docker" ]]; then
  VIASH_CMD="docker run --entrypoint=bash ${VIASH_DOCKER_RUN_ARGS[@]} ${VIASH_UNIQUE_MOUNTS[@]} $VIASH_DOCKER_IMAGE_ID"
fi


# set dependency paths


ViashDebug "Running command: $(echo $VIASH_CMD)"
cat << VIASHEOF | eval $VIASH_CMD
set -e
tempscript=\$(mktemp "$VIASH_META_TEMP_DIR/viash-run-bcftools_annotate-XXXXXX").sh
function clean_up {
  rm "\$tempscript"
}
function interrupt {
  echo -e "\nCTRL-C Pressed..."
  exit 1
}
trap clean_up EXIT
trap interrupt INT SIGINT
cat > "\$tempscript" << 'VIASHMAIN'
#!/bin/bash

## VIASH START
# The following code has been auto-generated by Viash.
$( if [ ! -z ${VIASH_PAR_INPUT+x} ]; then echo "${VIASH_PAR_INPUT}" | sed "s#'#'\"'\"'#g;s#.*#par_input='&'#" ; else echo "# par_input="; fi )
$( if [ ! -z ${VIASH_PAR_OUTPUT+x} ]; then echo "${VIASH_PAR_OUTPUT}" | sed "s#'#'\"'\"'#g;s#.*#par_output='&'#" ; else echo "# par_output="; fi )
$( if [ ! -z ${VIASH_PAR_ANNOTATIONS+x} ]; then echo "${VIASH_PAR_ANNOTATIONS}" | sed "s#'#'\"'\"'#g;s#.*#par_annotations='&'#" ; else echo "# par_annotations="; fi )
$( if [ ! -z ${VIASH_PAR_COLUMNS+x} ]; then echo "${VIASH_PAR_COLUMNS}" | sed "s#'#'\"'\"'#g;s#.*#par_columns='&'#" ; else echo "# par_columns="; fi )
$( if [ ! -z ${VIASH_PAR_COLUMNS_FILE+x} ]; then echo "${VIASH_PAR_COLUMNS_FILE}" | sed "s#'#'\"'\"'#g;s#.*#par_columns_file='&'#" ; else echo "# par_columns_file="; fi )
$( if [ ! -z ${VIASH_PAR_EXCLUDE+x} ]; then echo "${VIASH_PAR_EXCLUDE}" | sed "s#'#'\"'\"'#g;s#.*#par_exclude='&'#" ; else echo "# par_exclude="; fi )
$( if [ ! -z ${VIASH_PAR_FORCE+x} ]; then echo "${VIASH_PAR_FORCE}" | sed "s#'#'\"'\"'#g;s#.*#par_force='&'#" ; else echo "# par_force="; fi )
$( if [ ! -z ${VIASH_PAR_HEADER_LINE+x} ]; then echo "${VIASH_PAR_HEADER_LINE}" | sed "s#'#'\"'\"'#g;s#.*#par_header_line='&'#" ; else echo "# par_header_line="; fi )
$( if [ ! -z ${VIASH_PAR_HEADER_LINES+x} ]; then echo "${VIASH_PAR_HEADER_LINES}" | sed "s#'#'\"'\"'#g;s#.*#par_header_lines='&'#" ; else echo "# par_header_lines="; fi )
$( if [ ! -z ${VIASH_PAR_SET_ID+x} ]; then echo "${VIASH_PAR_SET_ID}" | sed "s#'#'\"'\"'#g;s#.*#par_set_id='&'#" ; else echo "# par_set_id="; fi )
$( if [ ! -z ${VIASH_PAR_INCLUDE+x} ]; then echo "${VIASH_PAR_INCLUDE}" | sed "s#'#'\"'\"'#g;s#.*#par_include='&'#" ; else echo "# par_include="; fi )
$( if [ ! -z ${VIASH_PAR_KEEP_SITES+x} ]; then echo "${VIASH_PAR_KEEP_SITES}" | sed "s#'#'\"'\"'#g;s#.*#par_keep_sites='&'#" ; else echo "# par_keep_sites="; fi )
$( if [ ! -z ${VIASH_PAR_MERGE_LOGIC+x} ]; then echo "${VIASH_PAR_MERGE_LOGIC}" | sed "s#'#'\"'\"'#g;s#.*#par_merge_logic='&'#" ; else echo "# par_merge_logic="; fi )
$( if [ ! -z ${VIASH_PAR_MARK_SITES+x} ]; then echo "${VIASH_PAR_MARK_SITES}" | sed "s#'#'\"'\"'#g;s#.*#par_mark_sites='&'#" ; else echo "# par_mark_sites="; fi )
$( if [ ! -z ${VIASH_PAR_MIN_OVERLAP+x} ]; then echo "${VIASH_PAR_MIN_OVERLAP}" | sed "s#'#'\"'\"'#g;s#.*#par_min_overlap='&'#" ; else echo "# par_min_overlap="; fi )
$( if [ ! -z ${VIASH_PAR_NO_VERSION+x} ]; then echo "${VIASH_PAR_NO_VERSION}" | sed "s#'#'\"'\"'#g;s#.*#par_no_version='&'#" ; else echo "# par_no_version="; fi )
$( if [ ! -z ${VIASH_PAR_OUTPUT_TYPE+x} ]; then echo "${VIASH_PAR_OUTPUT_TYPE}" | sed "s#'#'\"'\"'#g;s#.*#par_output_type='&'#" ; else echo "# par_output_type="; fi )
$( if [ ! -z ${VIASH_PAR_PAIR_LOGIC+x} ]; then echo "${VIASH_PAR_PAIR_LOGIC}" | sed "s#'#'\"'\"'#g;s#.*#par_pair_logic='&'#" ; else echo "# par_pair_logic="; fi )
$( if [ ! -z ${VIASH_PAR_REGIONS+x} ]; then echo "${VIASH_PAR_REGIONS}" | sed "s#'#'\"'\"'#g;s#.*#par_regions='&'#" ; else echo "# par_regions="; fi )
$( if [ ! -z ${VIASH_PAR_REGIONS_FILE+x} ]; then echo "${VIASH_PAR_REGIONS_FILE}" | sed "s#'#'\"'\"'#g;s#.*#par_regions_file='&'#" ; else echo "# par_regions_file="; fi )
$( if [ ! -z ${VIASH_PAR_REGIONS_OVERLAP+x} ]; then echo "${VIASH_PAR_REGIONS_OVERLAP}" | sed "s#'#'\"'\"'#g;s#.*#par_regions_overlap='&'#" ; else echo "# par_regions_overlap="; fi )
$( if [ ! -z ${VIASH_PAR_RENAME_ANNOTATIONS+x} ]; then echo "${VIASH_PAR_RENAME_ANNOTATIONS}" | sed "s#'#'\"'\"'#g;s#.*#par_rename_annotations='&'#" ; else echo "# par_rename_annotations="; fi )
$( if [ ! -z ${VIASH_PAR_RENAME_CHROMOSOMES+x} ]; then echo "${VIASH_PAR_RENAME_CHROMOSOMES}" | sed "s#'#'\"'\"'#g;s#.*#par_rename_chromosomes='&'#" ; else echo "# par_rename_chromosomes="; fi )
$( if [ ! -z ${VIASH_PAR_SAMPLES+x} ]; then echo "${VIASH_PAR_SAMPLES}" | sed "s#'#'\"'\"'#g;s#.*#par_samples='&'#" ; else echo "# par_samples="; fi )
$( if [ ! -z ${VIASH_PAR_SAMPLES_FILE+x} ]; then echo "${VIASH_PAR_SAMPLES_FILE}" | sed "s#'#'\"'\"'#g;s#.*#par_samples_file='&'#" ; else echo "# par_samples_file="; fi )
$( if [ ! -z ${VIASH_PAR_SINGLE_OVERLAPS+x} ]; then echo "${VIASH_PAR_SINGLE_OVERLAPS}" | sed "s#'#'\"'\"'#g;s#.*#par_single_overlaps='&'#" ; else echo "# par_single_overlaps="; fi )
$( if [ ! -z ${VIASH_PAR_REMOVE+x} ]; then echo "${VIASH_PAR_REMOVE}" | sed "s#'#'\"'\"'#g;s#.*#par_remove='&'#" ; else echo "# par_remove="; fi )
$( if [ ! -z ${VIASH_META_NAME+x} ]; then echo "${VIASH_META_NAME}" | sed "s#'#'\"'\"'#g;s#.*#meta_name='&'#" ; else echo "# meta_name="; fi )
$( if [ ! -z ${VIASH_META_FUNCTIONALITY_NAME+x} ]; then echo "${VIASH_META_FUNCTIONALITY_NAME}" | sed "s#'#'\"'\"'#g;s#.*#meta_functionality_name='&'#" ; else echo "# meta_functionality_name="; fi )
$( if [ ! -z ${VIASH_META_RESOURCES_DIR+x} ]; then echo "${VIASH_META_RESOURCES_DIR}" | sed "s#'#'\"'\"'#g;s#.*#meta_resources_dir='&'#" ; else echo "# meta_resources_dir="; fi )
$( if [ ! -z ${VIASH_META_EXECUTABLE+x} ]; then echo "${VIASH_META_EXECUTABLE}" | sed "s#'#'\"'\"'#g;s#.*#meta_executable='&'#" ; else echo "# meta_executable="; fi )
$( if [ ! -z ${VIASH_META_CONFIG+x} ]; then echo "${VIASH_META_CONFIG}" | sed "s#'#'\"'\"'#g;s#.*#meta_config='&'#" ; else echo "# meta_config="; fi )
$( if [ ! -z ${VIASH_META_TEMP_DIR+x} ]; then echo "${VIASH_META_TEMP_DIR}" | sed "s#'#'\"'\"'#g;s#.*#meta_temp_dir='&'#" ; else echo "# meta_temp_dir="; fi )
$( if [ ! -z ${VIASH_META_CPUS+x} ]; then echo "${VIASH_META_CPUS}" | sed "s#'#'\"'\"'#g;s#.*#meta_cpus='&'#" ; else echo "# meta_cpus="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_B+x} ]; then echo "${VIASH_META_MEMORY_B}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_b='&'#" ; else echo "# meta_memory_b="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_KB+x} ]; then echo "${VIASH_META_MEMORY_KB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_kb='&'#" ; else echo "# meta_memory_kb="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_MB+x} ]; then echo "${VIASH_META_MEMORY_MB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_mb='&'#" ; else echo "# meta_memory_mb="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_GB+x} ]; then echo "${VIASH_META_MEMORY_GB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_gb='&'#" ; else echo "# meta_memory_gb="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_TB+x} ]; then echo "${VIASH_META_MEMORY_TB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_tb='&'#" ; else echo "# meta_memory_tb="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_PB+x} ]; then echo "${VIASH_META_MEMORY_PB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_pb='&'#" ; else echo "# meta_memory_pb="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_KIB+x} ]; then echo "${VIASH_META_MEMORY_KIB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_kib='&'#" ; else echo "# meta_memory_kib="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_MIB+x} ]; then echo "${VIASH_META_MEMORY_MIB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_mib='&'#" ; else echo "# meta_memory_mib="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_GIB+x} ]; then echo "${VIASH_META_MEMORY_GIB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_gib='&'#" ; else echo "# meta_memory_gib="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_TIB+x} ]; then echo "${VIASH_META_MEMORY_TIB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_tib='&'#" ; else echo "# meta_memory_tib="; fi )
$( if [ ! -z ${VIASH_META_MEMORY_PIB+x} ]; then echo "${VIASH_META_MEMORY_PIB}" | sed "s#'#'\"'\"'#g;s#.*#meta_memory_pib='&'#" ; else echo "# meta_memory_pib="; fi )

## VIASH END

# Exit on error
set -eo pipefail

# Unset parameters
unset_if_false=(
    par_force
    par_keep_sites
    par_no_version
    par_single_overlaps
)

for par in \${unset_if_false[@]}; do
    test_val="\${!par}"
    [[ "\$test_val" == "false" ]] && unset \$par
done

# Execute bcftools annotate with the provided arguments
bcftools annotate \\
    \${par_annotations:+-a "\$par_annotations"} \\
    \${par_columns:+-c "\$par_columns"} \\
    \${par_columns_file:+-C "\$par_columns_file"} \\
    \${par_exclude:+-e "\$par_exclude"} \\
    \${par_force:+--force} \\
    \${par_header_line:+-H "\$par_header_line"} \\
    \${par_header_lines:+-h "\$par_header_lines"} \\
    \${par_set_id:+-I "\$par_set_id"} \\
    \${par_include:+-i "\$par_include"} \\
    \${par_keep_sites:+-k} \\
    \${par_merge_logic:+-l "\$par_merge_logic"} \\
    \${par_mark_sites:+-m "\$par_mark_sites"} \\
    \${par_min_overlap:+--min-overlap "\$par_min_overlap"} \\
    \${par_no_version:+--no-version} \\
    \${par_samples_file:+-S "\$par_samples_file"} \\
    \${par_output_type:+-O "\$par_output_type"} \\
    \${par_pair_logic:+--pair-logic "\$par_pair_logic"} \\
    \${par_regions:+-r "\$par_regions"} \\
    \${par_regions_file:+-R "\$par_regions_file"} \\
    \${par_regions_overlap:+--regions-overlap "\$par_regions_overlap"} \\
    \${par_rename_annotations:+--rename-annots "\$par_rename_annotations"} \\
    \${par_rename_chromosomes:+--rename-chrs "\$par_rename_chromosomes"} \\
    \${par_samples:+-s "\$par_samples"} \\
    \${par_single_overlaps:+--single-overlaps} \\
    \${par_threads:+--threads "\$par_threads"} \\
    \${par_remove:+-x "\$par_remove"} \\
    -o \$par_output \\
    \$par_input
    

    
VIASHMAIN
bash "\$tempscript" &
wait "\$!"

VIASHEOF


if [[ "$VIASH_ENGINE_TYPE" == "docker" ]]; then
  # strip viash automount from file paths
  
  if [ ! -z "$VIASH_PAR_INPUT" ]; then
    unset VIASH_TEST_INPUT
    IFS=';'
    for var in $VIASH_PAR_INPUT; do
      unset IFS
      if [ -z "$VIASH_TEST_INPUT" ]; then
      VIASH_TEST_INPUT="$(ViashDockerStripAutomount "$var")"
    else
      VIASH_TEST_INPUT="$VIASH_TEST_INPUT;""$(ViashDockerStripAutomount "$var")"
    fi
    done
    VIASH_PAR_INPUT="$VIASH_TEST_INPUT"
  fi
  if [ ! -z "$VIASH_PAR_OUTPUT" ]; then
    VIASH_PAR_OUTPUT=$(ViashDockerStripAutomount "$VIASH_PAR_OUTPUT")
  fi
  if [ ! -z "$VIASH_PAR_ANNOTATIONS" ]; then
    VIASH_PAR_ANNOTATIONS=$(ViashDockerStripAutomount "$VIASH_PAR_ANNOTATIONS")
  fi
  if [ ! -z "$VIASH_PAR_COLUMNS_FILE" ]; then
    VIASH_PAR_COLUMNS_FILE=$(ViashDockerStripAutomount "$VIASH_PAR_COLUMNS_FILE")
  fi
  if [ ! -z "$VIASH_PAR_HEADER_LINES" ]; then
    VIASH_PAR_HEADER_LINES=$(ViashDockerStripAutomount "$VIASH_PAR_HEADER_LINES")
  fi
  if [ ! -z "$VIASH_PAR_REGIONS_FILE" ]; then
    VIASH_PAR_REGIONS_FILE=$(ViashDockerStripAutomount "$VIASH_PAR_REGIONS_FILE")
  fi
  if [ ! -z "$VIASH_PAR_RENAME_ANNOTATIONS" ]; then
    VIASH_PAR_RENAME_ANNOTATIONS=$(ViashDockerStripAutomount "$VIASH_PAR_RENAME_ANNOTATIONS")
  fi
  if [ ! -z "$VIASH_PAR_RENAME_CHROMOSOMES" ]; then
    VIASH_PAR_RENAME_CHROMOSOMES=$(ViashDockerStripAutomount "$VIASH_PAR_RENAME_CHROMOSOMES")
  fi
  if [ ! -z "$VIASH_PAR_SAMPLES_FILE" ]; then
    VIASH_PAR_SAMPLES_FILE=$(ViashDockerStripAutomount "$VIASH_PAR_SAMPLES_FILE")
  fi
  if [ ! -z "$VIASH_META_RESOURCES_DIR" ]; then
    VIASH_META_RESOURCES_DIR=$(ViashDockerStripAutomount "$VIASH_META_RESOURCES_DIR")
  fi
  if [ ! -z "$VIASH_META_EXECUTABLE" ]; then
    VIASH_META_EXECUTABLE=$(ViashDockerStripAutomount "$VIASH_META_EXECUTABLE")
  fi
  if [ ! -z "$VIASH_META_CONFIG" ]; then
    VIASH_META_CONFIG=$(ViashDockerStripAutomount "$VIASH_META_CONFIG")
  fi
  if [ ! -z "$VIASH_META_TEMP_DIR" ]; then
    VIASH_META_TEMP_DIR=$(ViashDockerStripAutomount "$VIASH_META_TEMP_DIR")
  fi
fi


# check whether required files exist
if [ ! -z "$VIASH_PAR_OUTPUT" ] && [ ! -e "$VIASH_PAR_OUTPUT" ]; then
  ViashError "Output file '$VIASH_PAR_OUTPUT' does not exist."
  exit 1
fi


exit 0
