Frequently Asked Questions
WES / WGS
What reference genomes are used?
CBIcall supports:
- GRCh37 / b37: GATK-compatible reference genome
- GRCh38 / hg38: GATK-compatible reference genome
For native WES/WGS workflows, support depends on the selected software stack. This table describes CBIcall-shipped workflows, not every possible use of GATK upstream:
| Stack | CBIcall WES reference | CBIcall WGS support/reference | Main use |
|---|---|---|---|
gatk-3.5 | b37 | No CBIcall WGS workflow | Legacy Bash WES workflows and mtDNA prerequisites |
gatk-4.6 | b37 | b37 or hg38 | Current native WES/WGS workflows |
Which exome regions are used for WES?
Native WES interval resources differ between the legacy and current stacks:
| Stack | WES interval resource | Used for |
|---|---|---|
gatk-3.5 | Agilent SureSelect hg19 BED files (hg19.chr*.bed and flanked hg19.chr*.flank100bp.bed) | BQSR, recalibrated BAM generation, UnifiedGenotyper variant calling, and coverage summaries |
gatk-4.6 | GATK bundle / Broad b37 exome interval list (b37_Broad.human.exome.b37.interval_list) | BQSR, HaplotypeCaller, GenotypeGVCFs, and coverage summaries for WES mode |
WGS mode does not use WES capture intervals.
mtDNA (MToolBox)
What reference genome is used?
mtDNA workflows use RSRS (rsrs), the Reconstructed Sapiens Reference Sequence.
VCF vs. prioritized variants allele notation
In rare cases, the allele reported in prioritized_variants.txt may differ from the ALT allele reported in the VCF.
The Variant_Allele column is generated during annotation and prioritization and does not always follow VCF semantics, where ALT is defined relative to the mapping reference, such as RSRS.
What does GT=1 mean in results?
In variant reports, the Genotype (GT) field shows the observed allele using VCF allele indices:
0= reference allele1= first alternate allele2,3, ... = additional alternate alleles in multiallelic records
For CBIcall mtDNA reports, GT=1 means an ALT allele was detected in that sample.
Biological interpretation should use:
HF: fraction of reads supporting the ALT alleleDP: total read depth at the variant position
For mtDNA, GT tells you which allele was detected, not how much of it was detected. Use HF and DP to interpret heteroplasmy or homoplasmy.
General
What is the difference between native and external workflows?
Native CBIcall workflows are maintained in the cbicall-core collection and launched through a supported workflow backend: Bash, Snakemake, Nextflow, or Cromwell. They use the CBIcall project layout, provenance files, run reports, and usually the CBIcall resource bundle. Here, CBIcall is the orchestration framework and cbicall-core identifies the bundled workflow collection.
External workflows are third-party workflows registered in CBIcall. Today this means selected nf-core workflows launched through Nextflow. CBIcall validates the YAML contract, pins the registered workflow, and writes provenance and run reports, while nf-core keeps its own output layout, profiles, containers, and reference-resource assumptions.
See Workflows, External nf-core, and Resource Validation.
How do I set up cbicall on an HPC system?
On most HPC systems, Docker is not available. CBIcall is designed to run with Apptainer, formerly Singularity, which is the recommended approach for HPC environments.
Apptainer can execute Docker images directly, requires no root privileges, and integrates cleanly with batch schedulers.
In this setup:
- the container image is read-only
- the parameters YAML and input data remain in a host project directory
- the CBIcall resource bundle, when needed, is downloaded outside the container and bind-mounted at runtime
Recommended workflow:
- Pull the CBIcall container image using Apptainer.
- Download the CBIcall resource bundle for bundled
cbicall-coreworkflows. - Bind the project and resource directories when running the image.
Normal runs do not require a writable copy of the packaged workflow directory. For external nf-core workflows, running CBIcall with the site's native Nextflow and Apptainer modules is usually simpler than nesting execution inside the CBIcall image.
Do you have Slurm examples?
Yes. Complete templates are available for both installation models:
How do I cite CBIcall?
Please cite:
Rueda M, Fernandez-Orth D, Gut IG. CBIcall: a configuration-driven framework for variant calling in large sequencing cohorts. Bioinformatics Advances. 2026. doi:10.1093/bioadv/vbag232.