The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.
Register now for the upcoming GATK Best Practices workshop, Nov 7-8 at the Broad in Cambridge, MA. Open to all comers! More info and signup at


When I use EMIT_ALL_CONFIDENT_SITES for SNPs, I get an expected very large list of genotypes regardless if the genotypes vary from the reference. When I use the same command line but I switch the model to Indels, I only get a VCF of variant sites. Is the EMIT_ALL_CONFIDENT_SITES option not compatible with Indel discovery?

I'm grateful for any clarification.

Best Answer


  • dmittelmandmittelman Posts: 3Member

    Eric thanks for the super fast response. Let me just ask, when calling indels, how do I know which regions are callable but simply not variant from the reference?

  • ebanksebanks Broad InstitutePosts: 698Member, Administrator, Broadie, Moderator, Dev admin

    That's actually (unwittingly) a loaded question around us. I personally believe that this can be inferred from the depth (DP) of a region; although there are enough folks here who disagree with me that it is a problem we do plan on tackling in the future (as discussed in other similar threads, and with no promise as to an ETA).

    Eric Banks, PhD -- Director, Data Sciences and Data Engineering, Broad Institute of Harvard and MIT

  • dmittelmandmittelman Posts: 3Member

    Eric, thanks. This actually makes perfect sense, it is definitely not straightforward.

Sign In or Register to comment.