Hi GATK Users,

Happy Thanksgiving!
Our staff will be observing the holiday and will be unavailable from 22nd to 25th November. This will cause a delay in reaching out to you and answering your questions immediately. Rest assured we will get back to it on Monday November 26th. We are grateful for your support and patience.
Have a great holiday everyone!!!

GATK Staff

I have 50 exome samples belong to 25 families. Do I run GenotypeVCFs on familywise or 50 together?

NandaNanda CanadaMember

We have exome sequenced data for 50 samples in total for a cardiac disease. But they have been sequenced in different batches. Even some of the batches were 2 years old. We have relationship information available for these 50 samples. So these 50 samples have been grouped to 25 families, that is each family has 2 samples. Each family relationship can be any one of the following: siblings, sisters, brothers, father & son, and mother & daughter. **Currently, I have GVCFs available for 50 samples. **

As per the article "GATK Tutorial: Variant Callset Evaluation & Filtering", there are two requirements for Variant Quality Score Recalibration (VQSR)
1) GATK requires atleast 30 exome samples or more or 1 whole genome sample
2) Known variant databases

Case1: If I run GenotypeVCFs on each family wise, then I won't be able to filter using VQSR. I need to go for hard filtering. (because I have only 2 exome samples under each family)

Case2: If I run GenotypeVCFs on 50 samples together, then I can filter using VQSR.

Do I need to run "GenotypeVCFs (Joint Calling)" on each family individually or 50 samples together?
If I opt for case2, won't I miss family specific mutations?

Best Answer


  • NandaNanda CanadaMember
    edited April 2017

    Thanks, Sheila. I read the article you provided. Last Thursday, I started GenotypeGVCFs for 50 samples together. But I didn't mention the "--useNewAFCalculator". Then I will submit another job mentioning this parameter. What is the significance of the new QUAL calculated from "usenewAFcalculator" option?

    Next step, I am running VQSR for the raw VCF generated for my exome samples as per the GATK best practices. For annotations of variants,
    1) I read that "DP-Depth of Coverage" should not be used for exome datasets.
    2) Also, the "InbreedingCoeff" requires at least 10 samples to be computed. I have 50 samples (25 families), but there is another line mentioning that I should omit this annotation
    - if I have fewer samples or
    - if I have closely related samples (such as a family). In my case, I have 25 families with different relationships aforementioned.
    Therefore, I omitted following options DP and InbreedingCoeff from my command. Is my understanding correct?

  • SheilaSheila Broad InstituteMember, Broadie, Moderator admin
Sign In or Register to comment.