Generate an index for a PON file. This is a prerequisite for using the "--pon" option in mutectcaller.

prepon requires that the Contig field be present in the header of the input .vcf.gz file in order to do memory allocation at the start of execution. This field should include both the chromosome name and length:

##contig=<ID=chr1,length=248956422> ##contig=<ID=chr2,length=242193529> ##contig=<ID=chr3,length=198295559> ...

If your input .vcf.gz file does not include these value, run this command to update the header before running prepon:

bcftools reheader --fai YOUR_REFERENCE_FILE.fa.fai INPUT_PON.vcf.gz > UPDATED_PON.vcf.gz

For further information visit the prepon help page.

Quick Start

# This command assumes all the inputs are in <INPUT_DIR> and all the outputs go to <OUTPUT_DIR>.
$ docker run --rm --gpus all --volume <INPUT_DIR>:/workdir --volume <OUTPUT_DIR>:/outputdir
    -w /workdir \
    nvcr.io/nvidia/clara/clara-parabricks:<VERSION-TAG> \
    pbrun prepon \
    --in-pon-file /workdir/${INPUT_PON_VCF}