Holiday Notice:
The Frontline Support team will be slow to respond December 17-18 due to an institute-wide retreat and offline December 22- January 1, while the institute is closed. Thank you for your patience during these next few weeks. Happy Holidays!

Haplotypecaller: Problem with using ploidy argument.

nasimstnasimst Leuven, BelgiumMember

Dear GATK team,

I tried to fix the ploidy issue with sex chromosome in my pipeline, based on the information that I got from the previously asked questions.
However, when I run that, there is No difference in my output with -ploidy 1 and the default. I tried it also with nightly version but it did not work also.

My commands:

cmd="/usr/bin/time -a -v -o ${samplePrefix}.xy.timing java \
-Xmx10g -jar ${gatkPath}/GenomeAnalysisTK.jar \
-T HaplotypeCaller \
-R ${refGenome} \
-I ${sampleInputPrefix}.${fileExt} \
--emitRefConfidence GVCF \
--variant_index_type LINEAR \
--variant_index_parameter 128000 \
--dbsnp ${dbSNPall} \
--sample_ploidy 1
-L ${targetxy} \
-o ${samplePrefix}.${build}.xy.HC.g.vcf.gz \
-nct ${threads} \
-log ${samplePrefix}.xy.log"

Thanks!

Answers

Sign In or Register to comment.