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

GenotypeGVCFs wgs samples

I have about 130 wgs vcf files (created by GATK HaplotypeCaller). When I try to use GenotypeGVCFs, it goes out of memory. Can I split the sample vcf files by chromosome, call using GenotypeGVCFs, and then combine the chromowise results?

For example, can I do something like (for chr5, where ‘mysamplelist.list’ contains filepaths for sample vcf files):

java -Xmx32g -jar GenomeAnalysisTK.jar —L 5 -R myref.fa -T GenotypeGVCFs —variant mysamplelist.list -o out_chr5.vcf

and then combine the ‘out_chr{1..22}.vcf’ files using CombineGVCFs?

Or how else should I go about doing this?

thanks

Tagged:

Answers

Sign In or Register to comment.