Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

Output files can't be found

0

6 comments

  • Avatar
    David Benjamin

    L. A. Liggett The line "java.lang.IllegalArgumentException: Cannot construct fragment from more than two reads at. . ." in your stack trace strongly suggests that running the most recent GATK release, 4.1.5.0, will resolve the issue.  Please let us know if this does not fix it!

    0
    Comment actions Permalink
  • Avatar
    L. A. Liggett

    David Benjamin It looks like the workspace I'm using was updated to 4.1.4.1. And I'm using v2.7.0 of the Mutect2_PON, but I do not see any way of adjusting the GATK version. Is this something that needs to and can be adjusted?

    0
    Comment actions Permalink
  • Avatar
    David Benjamin

    The ideal thing to do is to run Mutect2 as a WDL in Terra, not as a notebook.  In that case you would just set the gatk_docker input to the 4.1.5.0 docker.  However, if you want a quick fix to get things working within the notebook, here's what to do:

    The notebook's "Set Runtime Values" section mentions that you need to configure your notebook runtime with the startup script gs://gatk-tutorials/scripts/gatk_4141.sh.  Download this script and edit it by replacing the line "GATK_VERSION=4.1.4.1" with "GATK_VERSION=4.1.5.0".  Upload the edited script to your own google bucket and set the notebook runtime to use the path of your edited script.

    0
    Comment actions Permalink
  • Avatar
    L. A. Liggett

    David Benjamin I am running the WDL, not an associated jupyter notebook. And I am setting the gatk_docker to `workspace.gatk_docker`. And the field isn't suggesting an input for the 4.1.5.0 docker; how to I refer to it?

    0
    Comment actions Permalink
  • Avatar
    David Benjamin

    The workflow is referring to a global workspace variable.  You need to set the workspace gatk_docker variable under the "Data" tab —> Other Data —> Workspace Data.  

    0
    Comment actions Permalink
  • Avatar
    L. A. Liggett

    David Benjamin This seems to have fixed the issue. The job has been running for 2.5 hours at this point so hopefully works. Thanks.

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk