To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at https://software.broadinstitute.org/firecloud/documentation/freecredits

HaplotypeCaller with WARNNING DepthPerSampleHC

Hello, sorry for my bad English.
While running GATK 4.0.1.1 Haplotypecaller with two or more bams together, this tool would send lots of lines of warning message about DepthPerSampleHC like follows:

10:32:58.674 WARN DepthPerSampleHC - Annotation will not be calculated, genotype is not called or alleleLikelihoodMap is null
10:33:02.097 WARN DepthPerSampleHC - Annotation will not be calculated, genotype is not called or alleleLikelihoodMap is null
10:33:02.102 WARN DepthPerSampleHC - Annotation will not be calculated, genotype is not called or alleleLikelihoodMap is null
10:33:02.114 WARN DepthPerSampleHC - Annotation will not be calculated, genotype is not called or alleleLikelihoodMap is null

And all those messages are gone while using "--annotation-group StandardAnnotation", which will remove "DP" in format. Maybe it due to the "./." genotype in some samples.
So, could I keep "DP" in vcf.gz without receiving too many WARN ?
Many thanks!

Here is my running command :

gatk HaplotypeCaller --java-options -Xmx40g 
-R GRCh37_latest_chr_genomic.fna 
--dbsnp All_20170710.vcf.gz 
-O pop.vcf.gz 
-I B1700.refine.bam 
-I NC.refine.bam 
--standard-min-confidence-threshold-for-calling 30 
1>pop.HaplotypeCaller.log 

Best Answer

Answers

Sign In or Register to comment.