This is a test version of Biostars. For the public version, visit https://www.biostars.org.
haplotypecaller can not run in single bam file

Hi all,

I an new to WGS analysis, i am trying to run haplotypecaller in one bam file. 

java -jar GenomeAnalysisTK.jar -T HaplotypeCaller -R GRCh37/hs37d5.fa -I LG_GWG.1.bam  -D reference/gatkbundle/dbsnp_138.b37.vcf -o LG_GWG.output.g.vcf -ERC GVCF -nct 4

can anyone help me to figure out this issue. thanks

the error that put out :

ERROR MESSAGE: Invalid command line: Argument emitRefConfidence has a bad value: Can only be used in single sample mode currently. Use the sample_name argument to run on a single sample out of a multi-sample BAM file
haplotypecaller gatk3 -erc wgs
• 5,880 views • link
written by deepaksrna • deepaksrna

It works without GVCFspecification, but to generate the gvcf file I need the GVCF argument

• 1 views • link
written by deepaksrna • deepaksrna

Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized.

• 1 views • link
written by GenoMax GenoMax

What does that mean? Can you provide the actual commandline you used?

• 1 views • link
written by GenoMax GenoMax

1 answer

try omitting the -ERC GVCF argument or provide the sample name to sample_name argument if it's a multi-sample BAM file

• 1 views • link
written by egeulgen egeulgen

Log in to answer this question.

More posts like this