HaplotypeCaller output mode
1
0
Entering edit mode
5 months ago

Hi,

I used GATK HaplotypeCaller to call both variant and non-variant sites in a specific region with the code below:

gatk HaplotypeCaller --heterozygosity 0.004 -ERC BP_RESOLUTION --native-pair-hmm-threads 16 -L 12:20140000-23600000 -I $bamdir/${i}.bam -O $outdir/${i}.vcf -R $refseq

Then I tried the other two modes: EMIT_ALL_ACTIVE_SITES, EMIT_ALL_CONFIDENT_SITES based on the code above, separately. For both of tests, I got 237 more sites with and alternative allele in the ATL column compared to the output-mode EMIT_VARIANTS_ONLY.

My objective is to do the calling to ensure specific coordinates have NO SNPs, and identify the possibility of mapping failure. Bearing this objective, which output-mode is a better choice?

Thanks in advance!

GATK HaplotypeCaller • 360 views
ADD COMMENT
0
0
Entering edit mode

Yes, thank you for your reply. I understand I have to use this parameter when I do joint genotyping using GenotypeGVCFs. I didn't make my words clear in the question, my current step is doing variant calling and getting GVCF for each sample, separately, and then I will combine and joint genotyping for a single VCF. I am quite confused that at my current step, except from adding BP_RESOLUTION, do I have to include output-mode as ACITVE/CONFIDENT?

ADD REPLY

Login before adding your answer.

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