Description

Subworkflow to compute genotype likelihoods from BAM files using bcftools/mpileup and bcftools/call variants by chromosomes. The resulting VCF files are then merged by chromosomes and annotated with bcftools/annotate to give an ID for each variants. During samples merging and region concatenation, all metadata will be preserved and stored in the key metas at each step.

Input

name
description
pattern

ch_bam

Channel with input data

ch_posfile

Channel with position to call variants by chromosomes

ch_fasta

Channel with reference genome data

meta_sample_merge_key

List of keys that define the sample.
All file sharing the same values for these keys will be merged together.
e.g. [ “id” ]

meta_sample_merge_value

Shared value that will be set to meta_sample_merge_key for all samples to be merged together.
e.g. “all_samples”

meta_region_gather_keys

List of keys to be used to gather the VCF files by region with vcf_gather_bcftools.
All file sharing the same values for these keys will be gathered together.
e.g. [ “id”, “panel_id” ]

sort_region_gather

Whether to sort the VCF files by region after concatenation with
vcf_gather_bcftools.

annotate

Whether to annotate the merged VCF file with bcftools/annotate to give
an ID for each variants

Output

name
description
pattern

multiqc_files

Channel containing stat files of bcftools mpileup and call

vcf_index

Channel with one VCF files by chromosomes