The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

#### Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

#### ☞ Get notifications!

You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

#### ☞ Formatting tip!

Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks (  ) each to make a code block as demonstrated here.

##### Jump to another community
Picard 2.9.0 is now available. Download and read release notes here.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

# PlotGenotypingResults

Cambridge, MAMember, Administrator, Broadie Posts: 11,669 admin
edited September 2012

### 1. Introduction

The PlotGenotypingResults utility produces plots that help to visualize the
genotyping data for individual sites.

The plots show a histogram of read depth, decorated with information about
aberrant read pairs and/or breakpoint spanning reads and showing the called
genotype for each sample. This utility requires information from auxilliary
data files produced during genotyping which are written to the genotyping run
directory.

### 2. Inputs / Arguments

• -site &lt;site-list&gt; : The ID of a genotyped site (required). This should
match the third column in the output vcf file. This argument can be supplied
multiple times. It can also be a comma-separated list of sites or a file (with
extension .list) containing a list of sites, one site per line.

• -runDirectory &lt;dir&gt; : The run directory from a genotyping run which
contains the auxilliary output files. This is the typical way to specify the
location of the auxilliary output files.

• -partitionMapFile &lt;file&gt; : Full path to the partition map file that maps
site IDs to the partition used during genotyping (when genotyping is run in
parallel). The default value is partition.genotypes.map.dat in the specified
run directory. If you did not use parallel genotyping, you need to supply
-auxFilePrefix instead.

• -auxFilePrefix &lt;string&gt; : The path prefix to the auxilliary output files
for the listed sites. This is the part of the path without the suffixes such
as .genotypes.vcf or .genotypes.gts.dat. You can use either the auxilliary
files for a single partition (if all sites to plot come from that partition)
or the merged auxilliary files for the entire run, but the latter requires
more I/O and memory and is recommended only for small runs. You should supply
either the -runDirectory (and optionally the -partitionMapFile) argument or
the -auxFilePrefix argument.

• -genderMapFile &lt;map-file&gt; : Path to a map file specifying the gender of
each sample (optional). If supplied, the gender of each sample is shown for
sites on chrX and chrY by overlaying the sample with a colored dot (blue for
boys, red for girls). The map file should contain two tab-delimited columns,
SAMPLE and GENDER (M/F).

### 3. Outputs

• -O &lt;output-file&gt; : The destination output file (PDF) (required).

### 4. Running

java -Xmx2g -cp SVToolkit.jar:GenomeAnalysisTK.jar \
org.broadinstitute.sv.apps.PlotGenotypingResults \
-site DEL_P0001_1 \
-site my_site_file.list \
-O my_output_file.pdf \
-runDirectory genotyping_run_dir
`

Geraldine Van der Auwera, PhD

Post edited by Geraldine_VdAuwera on
Sign In or Register to comment.