Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

GenotypeGVCFs wgs samples

$
0
0

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


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>