Warning

Iceberg reaches end-of-life on 30th November 2020. If you are running jobs on Iceberg then you need to take urgent action to ensure that your jobs/scripts will run on ShARC or Bessemer. If you have never used ShARC or Bessemer then now is the time to test your scripts. Not all software on Iceberg is available on ShARC/Bessemer.

GATK¶

The Genome Analysis Toolkit or GATK is a software package for analysis of high-throughput sequencing data, developed by the Data Science and Data Engineering group at the Broad Institute. The toolkit offers a wide variety of tools, with a primary focus on variant discovery and genotyping as well as strong emphasis on data quality assurance. Its robust architecture, powerful processing engine and high-performance computing features make it capable of taking on projects of any size.

Interactive Usage¶

After connecting to Iceberg (see Establishing a SSH connection), start an interactive sesssion with the qsh or qrsh command.

The latest version of GATK (currently 3.4-46) is made available with the command

module load apps/binapps/GATK

Alternatively, you can load a specific version with

module load apps/binapps/GATK/3.4-46
module load apps/binapps/GATK/2.6-5

Version 3.4-46 of GATK also changes the environment to use Java 1.7 since this is required by GATK 3.4-46. An environment variable called GATKHOME is created by the module command that contains the path to the requested version of GATK.

Thus, you can run the program with the command

java -jar $GATKHOME/GenomeAnalysisTK.jar -h

Which will give a large amount of help, beginning with the version information

---------------------------------------------------------------------------------
The Genome Analysis Toolkit (GATK) v3.4-46-gbc02625, Compiled 2015/07/09 17:38:12
Copyright (c) 2010 The Broad Institute
For support and documentation go to http://www.broadinstitute.org/gatk
---------------------------------------------------------------------------------

Documentation¶

The GATK manual is available online https://www.broadinstitute.org/gatk/guide/

Installation notes¶

The entire install is just a .jar file. Put it in the install directory and you’re done.

Modulefile¶

Version 3.4-46

The module file is on the system at /usr/local/modulefiles/apps/binapps/GATK/3.4-46

Its contents are

#%Module1.0#####################################################################
##
## GATK 3.4-46 modulefile
##

## Module file logging
source /usr/local/etc/module_logging.tcl
##

#This version of GATK needs Java 1.7
module load apps/java/1.7

proc ModulesHelp { } {
        puts stderr "Makes GATK 3.4-46 available"
}


set version 3.4-46
set GATK_DIR /usr/local/packages6/apps/binapps/GATK/$version

module-whatis   "Makes GATK 3.4-46 available"

prepend-path GATKHOME $GATK_DIR

Version 2.6-5

The module file is on the system at /usr/local/modulefiles/apps/binapps/GATK/2.6-5

Its contents are

#%Module1.0#####################################################################
##
## GATK 2.6-5 modulefile
##

## Module file logging
source /usr/local/etc/module_logging.tcl
##

proc ModulesHelp { } {
        puts stderr "Makes GATK 3.4-46 available"
}


set version 2.6.5
set GATK_DIR /usr/local/packages6/apps/binapps/GATK/$version

module-whatis   "Makes GATK 2.6-5 available"

prepend-path GATKHOME $GATK_DIR