Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Problems using GATK for variant calling pooja_singh RNA Sequencing 1 07-07-2012 11:23 AM
Statistics behind variant/consensus calling ZHIHUA.LI Bioinformatics 1 06-28-2012 10:46 AM
GATK variant calling for multiple exome sehrrot Bioinformatics 2 04-24-2012 05:09 AM
GATK variant calling on uniquely mapped reads? pravee1216 Bioinformatics 5 03-26-2012 09:19 AM
GATK excludes some samples for cohort variant calling liu_xt005 Bioinformatics 2 02-01-2012 12:58 PM

Thread Tools
Old 01-23-2014, 06:24 AM   #1
Location: NJ

Join Date: Jun 2010
Posts: 34
Default GATK variant calling statistics...

Hi all,
We are calling variants on partial exome data using GATK. We want to see what fraction of our samples have a particular variant. As constituted, GATK will give us the set of samples containing the particular variant. Any sample that doesn't will simply exclude it from the VCF. My question is whether or not there is a way to distinguish the samples that have a true wild-type phenotype (and thus no variant call) from those that aren't called due to poor quality score (and thus should be excluded from our percentages since the phenotype cannot be determined). Thank you very much in advance for your help.
milesgr is offline   Reply With Quote
Old 01-23-2014, 10:01 AM   #2
Location: Boston, MA

Join Date: Apr 2012
Posts: 43

There are two ways to deal with this. One is to use the EMIT_ALL_SITES option of UnifiedGenotyper and enable calculation of reference PLs. The other, which is better, is to call variants with HalpotypeCaller (GATK's newer variant caller) and use the -ERC option to generate a gVCF.

For more details, see the online GATK documentation and please post any follow-up questions on the GATK support forum.
vdauwera is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 04:14 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO