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!

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.

#### ☞ Did you remember to?

1. Search using the upper-right search box, e.g. using the error message.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

#### ☞ 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.

GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

# Can I use UnifiedGenotyper to call somatic variants in tumor/control tissues?

Member Posts: 0
edited March 2013

Hi,

I am interested in the following scenario:
1. sequence tumor and control samples to separate fastq files.
2. Execute read alignment for the 2 samples separately.
3. Execute UnifiedGenotyper with the 2 BAM files (control and tumor) with the following command:

java -jar GenomeAnalysisTK.jar \
-T UnifiedGenotyper \
-R ReferenceGenome.fasta \
-I contro.bam -I tumor.bam \
-D /usr/sap/GenomicsPlatform/dbSNP/dbsnp_137.b37.vcf \
-o output.vcf

Is this a proper usage for UnifiedGenotyper?
How does UnifiedGenotyper refer to the 2 BAM files it recieves as input?
What do I expect to see in the output.vcf file? are they somatic variants which describe the variation between control and tumor samples?

Thank you,
Stas

• Member Posts: 4

Hello Geraldine,

How does the SomaticIndelDetector fits to overall picture? According to the document (http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_indels_SomaticIndelDetector.html) it can be used to detect indels in tumor-normal paired sample. Is it related to the workaround you mentioned above?

Thanks,
Amit

Hi Amit,

MuTect is designed to detect somatic SNPs while the SomaticIndelDetector is for somatic Indels. Please note that the SomaticIndelDetector is no longer included in GATK since version 2.4; that tool is now in the hands of a different group here at the Broad Institute. Specifically, the Cancer Group, who also produce MuTect. You can contact them to ask how you can obtain their tools for cancer-related analysis.

Geraldine Van der Auwera, PhD

• Member Posts: 1

Hi,
I was planning to get the somatic mutations by processing the normal and tumor samples separately then with a script filter out the mutations that are found in the normal sample.

Do you think, that works ?
Thanks

Suleyman