Bug Bulletin: The recent 3.2 release fixes many issues. If you run into a problem, please try the latest version before posting a bug report, as your problem may already have been solved.

UnifiedGenotyper --genotype_likelihoods_model INDEL --output_mode EMIT_ALL_CONFIDENT_SITES

vanheel123vanheel123 Posts: 5Member
edited January 2013 in Ask the GATK team

Hi,

When I use UnifiedGenotyper with --genotype_likelihoods_model SNP --output_mode EMIT_ALL_CONFIDENT_SITES I get the reference SNP homozygote calls (or ./. if insufficient depth/quality etc). Great!

But when I use UnifiedGenotyper with --genotype_likelihoods_model INDEL --output_mode EMIT_ALL_CONFIDENT_SITES I only get non-reference calls, everything else (i.e. reference homozygotes, and anything uncallable) is ./.

I want to be able to select variants (SNPs and INDELs) on call rate across samples - as one would do for array genotype data. And avoid case-control bias due to differential missingness.

thanks,

david van heel

Post edited by Geraldine_VdAuwera on

Comments

Sign In or Register to comment.