Hi GATK Users,

Happy Thanksgiving!
Our staff will be observing the holiday and will be unavailable from 22nd to 25th November. This will cause a delay in reaching out to you and answering your questions immediately. Rest assured we will get back to it on Monday November 26th. We are grateful for your support and patience.
Have a great holiday everyone!!!

GATK Staff

missing calls from VCF

1. I combined gVCFs of ~100 individuals. However, when I open the file there are missing genotypes. i.e.- not all rows have equal number of genotypes (see attached). I am using GATK_3.4.6 and the command I am using to merge the gVCFs is:
$cmd = "java -Xmx4g -jar $GATKdir/GenomeAnalysisTK.jar -T GenotypeGVCFs -R $genome --dbsnp $ResourcesDir/dbsnp_138.b37.vcf $filesToJoin -o example.vcf";

Where $filesToJoin is a string with -V individual.gvcf

  1. The total number of individuals is 400, and I combined every 100 together. Can I do the calibration with several -input arguments or must I merge all merged files into one?


Sign In or Register to comment.