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.

Did you remember to?


1. Search using the upper-right search box, e.g. using the error message.
2. Try the latest version of tools.
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.

Did we ask for a bug report?


Then follow instructions in Article#1894.

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.

Doubt Regarding Combine Variants and Multisample SNP calling

anupamanupam Member Posts: 8

Hi sir,
I used GATK combine variants to call SNP from my 7 vcf files. Now I am using unified genotyper to call all samples simultaneously for SNP.
When I focused on my gene of interest (DLX6), in combine variants output, there are 5 "./." among 7 samples. but in case of multisample SNP calling, DLX6 gave 0/0 in those samples and its corresponding information. So, I want to know what is the meaning of "./." in combine variants results and which results should I take for analysis?

Best Answer

Answers

  • anupamanupam Member Posts: 8

    Hi sir,
    I have attached my gene of interest DLX6 results came from combine variants and Multisample SNP calling. please look the DLX6 variants information in 7 samples. It's a two line results, so won't take your much time.

    Thanks in advance.

    xlsx
    xlsx
    DLX6_results.xlsx
    9K
  • anupamanupam Member Posts: 8

    Hi sir,
    Thank you for your reply. It clears my doubt regarding homozygous calls that are not coming in combine variants. But, there is some inconsistency in two of my samples. Younger son and wife sample contain heterozygous calls and satisfy other criteria (from multisample SNP calling),but it won't coming in combine variants results. I have attached my file and highlight the inconsistent samples in yellow. so please go through it once and let me know what's the reason for that inconsistency.

    Thank you for your time.

    xlsx
    xlsx
    DLX6_results.xlsx
    10K
  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie Posts: 11,631 admin

    Ah, this is a good thing actually. This can happen if the site is difficult to genotype in a sample, and it illustrates the power of multisample calling to fix such problem. In single sample calling, there is not enough evidence to call the variant, so the program decides the sample is hom-ref at this site, and does not emit a variant call (so CombineVariants will assign ./.). However, in multisample calling, the program sees more than one sample with some evidence for the variant. It would not be enough evidence to call the variant individually per sample, but taken together it is enough to call the variant in the cohort, so the program emits a variant call.

    This shows you how using multisample calling can help you discover variants that you would not be able to see in single sample calling.

    Geraldine Van der Auwera, PhD

Sign In or Register to comment.