Notice:
If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!

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.
Attention:
We will be out of the office on October 14, 2019, due to the U.S. holiday. We will return to monitoring the forum on October 15.

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.