Attention:
The frontline support team will be unavailable to answer questions on April 15th and 17th 2019. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!

GenotypeGVCFs not working for 4.beta.3

jjmmiijjmmii Member
edited August 2017 in Ask the GATK team

I am using GATK beta 4.3 (downloaded yesterday) on trio whole-exome data, following best practices on germline SNP/indel calling up to HaplotypeCaller. I was unable to follow the example command in GenotypeGVCFs tool documentation "Perform joint genotyping on a set of GVCFs listed in a text file, one per line":
gatk-launch GenotypeGVCFs -R ~/hg38/Homo_sapiens_assembly38.fasta -V input_gvcfs.list -O trio11.vcf
It exited with this error:

[August 11, 2017 2:27:38 PM HKT] org.broadinstitute.hellbender.tools.walkers.GenotypeGVCFs done. Elapsed time: 0.02 minutes.
Runtime.totalMemory()=1495793664
***********************************************************************

A USER ERROR has occurred: Cannot read input_gvcfs.list because no suitable codecs found

My input_gvcfs.list looks like this:

F11.g.vcf
M11.g.vcf
P11.g.vcf

I have also tried listing the VCFs one by one, following the example command "Perform joint genotyping on a set of GVCFs enumerated in the command line":
gatk-launch GenotypeGVCFs -R ~/hg38/Homo_sapiens_assembly38.fasta -V F11.g.vcf -V M11.g.vcf -V P11.g.vcf -O trio11.vcf
But it exited with this error:

A USER ERROR has occurred: Argument '[V, variant]' cannot be specified more than once.

Please help! Thank you very much.
-Jamie

Best Answer

Answers

Sign In or Register to comment.