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.

SelectVariants Starts Traversal but Does not Progress, High CPU Usage

Hi,

I am using the GATK tool SelectVariants to only select variants that have passed FilterMutectCalls. Both FilterMutectCalls and Mutect2 were run in multi-sample mode, so the VCF being input to SelectVariants is multisample. The GATK version is 4.1.0.0, HTSJDK Version 2.18.2, and Picard Version: 2.18.25.

The issue is the SelectVariants tool begins the traversal of the VCF, but does not make any progress. The log file stalls out as shown:
09:52:57.790 INFO SelectVariants - Done initializing engine 09:52:57.967 INFO ProgressMeter - Starting traversal 09:52:57.967 INFO ProgressMeter - Current Locus Elapsed Minutes Variants Processed Variants/Minute

This is the last line logged at which point the CPU usage jumps precipitously. I have alloted 10G of memory for the job and 1 hour of run-time, but there is till no progress.

The exact command I am using is:
$gatk_launcher --java-options -Xmx${mem}g SelectVariants \ -R $reference \ --variant $input_fn \ --output $output_fn \ --exclude-filtered true &>> $log_file

I ran ValidateVariants on the VCF in question, and all records validated. The size of the VCF is 18.7M, although I have been able to run the same command succesfully on larger multi-sample VCFs with the same command.

I have attached a subset of the records from the VCF that failed.

For now, I have circumvented the issue with a grep, but I figured I would point it out nonetheless.

Thank you!

Answers

Sign In or Register to comment.