process HAPLINK_VARIANTS {
    tag "${prefix}"
    label 'process_medium'
    label 'haplink'

    container 'ghcr.io/ksumngs/haplink.jl:1.0.0'

    input:
    tuple val(prefix), path(bam), path(bai)
    path reference

    output:
    tuple val(prefix), path("*.vcf")

    publishDir "${params.outdir}/${task.process}", mode: 'copy'

    when:
    task.ext.when == null || task.ext.when

    script:
    """
    export JULIA_NUM_THREADS=${task.cpus}
    haplink variants \\
            "${reference}" \\
            "${bam}" \\
        > "${prefix}.vcf"
    """
}