Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

Is VariantAnnotator compatible with the GVCF workflow?

$
0
0

Hi,

In the Best Practices for Germline SNP & Indel Discovery in Whole Genome and Exome Sequence, https://software.broadinstitute.org/gatk/best-practices/bp_3step.php?case=GermShortWGS&p=2,

the steps are listed as:

  1. PRE-PROCESSING
    Map to reference with BWA mem - Mark duplicates with Picard - Base quality score recalibration

  2. VARIANT DISCOVERY
    Generate GVCF per-sample with HaplotypeCaller - Perform joint genotyping - Filter variants

  3. CALLSET REFINEMENT
    (Optional) Refine genotypes - Annotate variants - Evaluate callset

However, after step 2, we have a multi-sample VCF file. But, in step 3, it is suggested to use VariantAnnotator to annotate the variants, but in the tool documentation it seems VariantAnnotator accepts VCF files (not GVCF files) and seems to assume they are single-sample VCF files.

Thanks,

Henry


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>