How to force GATK haplotypecaller to output all genomic position (even non variant positions)
1
1
Entering edit mode
18 months ago

Hi,

Is there a way to force GATK haplotypecaller to output all genomic positions in the vcf; not only variant ones ?

Example here :

Top : Raw bam file

Middle : -outbam using : gatk HaplotypeCaller -R $genome -O test.vcf -I in.bam -ERC GVCF --max-alternate-alleles 3 -bamout test.bam -L interval.bed --output-mode EMIT_ALL_ACTIVE_SITES

bottom : -outbam using gatk HaplotypeCaller -R $genome -O test.vcf -I in.bam -ERC GVCF --max-alternate-alleles 3 -bamout test.bam -L interval.bed --output-mode EMIT_ALL_CONFIDENT_SITES

You can see that the region without variant sites are not processed by haplotypecaller ; even with EMIT_ALL_ACTIVE_SITES and EMIT_ALL_CONFIDENT SITES.

Any idea how to tune haplotypecaller to output all sites with enough coverage ?

gatk version : 4.1.7.0

enter image description here

Thanks

GATK haplotypecaller • 1.7k views
ADD COMMENT
2
Entering edit mode
18 months ago
curious ▴ 620

gatk HaplotypeCaller -R {reference} -I {bam_path} -L {bed_path} -ERC BP_RESOLUTION -O {out_vcf}

Maybe look at BP_RESOLUTION

ADD COMMENT

Login before adding your answer.

Traffic: 2377 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