We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

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.