Test-drive the GATK tools and Best Practices pipelines on Terra


Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.

Mutect2 - Error: A USER ERROR has occurred: Badly formed genome unclippedLoc: Parameters to GenomeLo

Im trying to create a PON for the new mutect2 I got an error like this:

A USER ERROR has occurred: Badly formed genome unclippedLoc: Parameters to GenomeLocParser are incorrect:The genome loc coordinates 1-16571 exceed the contig size (16569)

This is the command i ran:
java -Dsamjdk.use_async_io_read_samtools=false -Dsamjdk.use_async_io_write_samtools=true -Dsamjdk.use_async_io_write_tribble=false -Dsamjdk.compression_level=2 -jar /gatk/gatk-package-4.1.2.0-local.jar Mutect2 -R /gscmnt/gc2560/core/model_data/2887491634/build21f22873ebe0486c8e6f69c15435aa96/all_sequences.fa -I /gscmnt/gc2648/wgs/model_data/fd6be422ed9a4172b4590bdad4b08163/build1700a85838d745cbafc1ece038394120/results/final.cram -tumor H_IJ-HG00514-HG00514_1 -L /gscmnt/gc2708/info/medseq/ppg/saimukund/CWLs/Mutect/my.interval_list -O 1_for_pon.vcf.gz

This is a sample of how my interval list looks like:
chr1:1-248956422
chr2:1-242193529
chr3:1-198295559
chr4:1-190214555
chr5:1-181538259
chr6:1-170805979
chr7:1-159345973
chr8:1-145138636
chr9:1-138394717
chr10:1-133797422
Sign In or Register to comment.