Call somatic variants with the accelerated MuSE variant caller.
The following runs MuSE in call mode:
$ pbrun muse \
--mode call \
--ref Ref/Homo_sapiens_assembly38.fasta \
--in-tumor-bam tumor.bam \
--in-normal-bam normal.bam \
--out-vcf out.vcf
pbrun muse will append .MuSE.txt to the end of the file name provided in
--out-file.
The following runs MuSE in sump mode:
$ pbrun muse \
--mode sump \
--in-callfile muse_call_results.txt \
--out-vcf muse_sump_results.vcf \
--datatype G \
--in-dbsnp input.dbsnp.vcf.gz
The command below is the original CPU counterpart of the Parabricks commands above. The output from this command will be identical to the output from the above command. See the Output Comparison page for comparing the results.
./MuSE call -f Ref/Homo_sapiens_assembly38.fasta -O out.vcf tumor.bam normal.bam
Call variants with accelerated Muse variant caller
Input/Output file options
- --in-tumor-bam IN_TUMOR_BAM
- --in-normal-bam IN_NORMAL_BAM
- --ref REF
- --in-dbsnp IN_DBSNP
- --in-callfile IN_CALLFILE
- --out-file OUT_FILE
- --out-vcf OUT_VCF
Path to the input tumor bam file. (default: None)
Path to the input normal bam file. (default: None)
Path to the input reference file. (default: None)
Path to the dbSNP vcf file that should be bgzip compressed, tabix indexed and based on the same reference genome used in call mode. (default: None)
Path to the output file of muse call command. (default: None)
Output stats file, .MuSE.txt will be automatically appended to this filename. (default: None)
Output vcf file. (default: None)
Options specific to this tool
- --datatype DATATYPE
- --mode MODE
- --num-threads NUM_THREADS
sample-specific error model to use. It can be either E (for whole exomes) or G (for whole genomes) (default: None)
Mode can be call or sump. User should first run muse with call mode and then run with sump mode to generate vcf. Please see MuSE somatic caller documentation for details. (default: None)
Option is required.
Number of threads to run. Use 32 threads for maximum performance if your system has enough CPU threads. (default: 12)
Common options:
- --logfile LOGFILE
- --tmp-dir TMP_DIR
- --with-petagene-dir WITH_PETAGENE_DIR
- --keep-tmp
- --license-file LICENSE_FILE
- --no-seccomp-override
- --version
Path to the log file. If not specified, messages will only be written to the standard error output. (default: None)
Full path to the directory where temporary files will be stored.
Full path to the PetaGene installation directory. By default, this should have been installed at /opt/petagene. Use of this option also requires that the PetaLink library has been preloaded by setting the LD_PRELOAD environment variable. Optionally set the PETASUITE_REFPATH and PGCLOUD_CREDPATH environment variables that are used for data and credentials (default: None)
Do not delete the directory storing temporary files after completion.
Path to license file license.bin if not in the installation directory.
Do not override seccomp options for docker (default: None).
View compatible software versions.