We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

Options for haploid

Hi,

I'm working with SNP calling in a bacterium - I don't have a set of known SNPs, so prior to recalibration, so generate a mask file from all data. My question is, what should I put for the following two options:

--snpHets
--indelHets

because the bacterium is haploid (and I specify --ploidy 1) it seems like these options should be set to zero, as there are no heterozygous loci, but I worry that if I set them to zero, that it won't work as expected. Any advice? I was just setting them to 0.001

Thanks,
Gavin

Tagged:

Best Answer

Answers

  • Gavin_SherlockGavin_Sherlock StanfordMember

    sorry - those are the options for the wrapper script. They go into GATK as -hets and -indelHeterozygosity.

    Thanks.

Sign In or Register to comment.