GATK 3.8 do not exclude annotations
1
0
Entering edit mode
2.2 years ago
Filago ▴ 100

Hello,

I am using GATK 3.8 and want to call variants in GVCF mode:

gatk3 -T HaplotypeCaller -R human_g1k_v37.fasta -I test.bam -ERC GVCF -variant_index_type LINEAR -variant_index_parameter 128000 -o abc.vcf -L region.sort.bed -A QualByDepth

However, no QD values show up and in the VCF header it says:

excludeAnnotation=[ChromosomeCounts, FisherStrand, StrandOddsRatio, QualByDepth]

It seems like the excludeAnnotation has priority over the -A arguments. Do you have any idea how to avoid excludeAnnotation?

Best,

Andreas

GATK3 • 695 views
ADD COMMENT
1
Entering edit mode
2.2 years ago
Filago ▴ 100

Ok figured it out myself! The -ERC GVCF option prohibits QD

ADD COMMENT
0
Entering edit mode

Please accept your own answer to resolve the post.

ADD REPLY

Login before adding your answer.

Traffic: 1154 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6