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.4 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.

VCF contrasts - workflow

I am just starting with GATK but even though I have looked and looked I can not find a simple walkthrough of having many VCFs and running a range of contrasts based on sample data. I guess this must be one of the most common used tools.
I have 200 VCFs that I want to contrast in many different ways depending on sample phenotype.
Is there anyone that can direct me to a good guide on how to do that?
Thanks

Tagged:

Answers

  • CarneiroCarneiro Charlestown, MAMember

    I am not sure I understand what you want to do, but it seems to me like you are talking about the Plink tools ? You can make queries in the GATK on vcfs using SelectVariants, VariantEval, GenotypeConcordance, CombineVariants and other tools.

  • Well, lets say I want to group my vcfs in high and low BMI, high and low LDL etc. and look for the SNVs that differentiate the groups.

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie

    That is definitely something you can do with VariantEval and GenotypeConcordance. Those tools have a lot of capabilities, you should look them up in the documentation.

    We don't currently provide a set workflow for this type of analysis because it there are so many possible variations depending on your study design.

Sign In or Register to comment.