Skip to content

Mismatches between VCF and ROH files when running with multiple samples #739

@fellen31

Description

@fellen31

Description of the bug

When running with multiple samples there's a risk of mixups of VCF and ROH files between samples in:

RHOCALL_VIZ(BCFTOOLS_VIEW_UNCOMPRESS.out.vcf, BCFTOOLS_ROH.out.roh)

Should be solved with a join+multiMap e.g:

BCFTOOLS_VIEW_UNCOMPRESS.out.vcf
        .join(BCFTOOLS_ROH.out.roh)
        .multiMap { meta, vcf, roh ->
            vcf: [meta, vcf]
            roh: [meta, roh]
        }
        .set { ch_rhocall_viz_input }

    RHOCALL_VIZ(
        ch_rhocall_viz_input.vcf,
        ch_rhocall_viz_input.roh,
    )

Command used and terminal output

Relevant files

No response

System information

No response

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions