Bug Bulletin: we have identified a bug that affects indexing when producing gzipped VCFs. This will be fixed in the upcoming 3.2 release; in the meantime you need to reindex gzipped VCFs using Tabix.

Reference calls by HaplotypeCaller

jianljianl Posts: 2Member

Hi there,

I'm trying to run the HaplotypeCaller in GATK2.5. I hope to get the confident reference calls as UnifiedGenotyper outputs, so I set the option "-out_mode EMIT_ALL_CONFIDENT_SITES". However, I notice the total number of reference calls HaplotypeCaller makes is really small. I'm wondering could you provide any information on how HaplotyperCaller decides where to/not to output the reference calls please?

Tagged:

Best Answer

Answers

  • jianljianl Posts: 2Member

    To be more specific, I'm working with single sample. I'm wondering is there any other option that I didn't correctly set or do I have to feed in multiple samples to get the ref call at positions where at least one sample showing variant?

    Any help would be appreciated!

    Thanks, Jian

Sign In or Register to comment.