Hello!
I'm trying to use GATK to call SNPs on multiple bam files (at once) but I encounter some "problems". I tried HaplotypeCaller and UnifiedGenotyper and it works well. The vcf is generated. But not in the way I expected:
CHROM ... FORMAT 20
...
01_contig_678504 ... GT:AD:DP:GQ:PL 1/1:1,249:251:99:11142,712,0
...
I expected to get one genotype per individual, instead of what I only got one sample/genotype named "20". I suppose it's normal and I probably missed something (down below is how I used these tools).
Is it possible to get all different genotypes (for 96 individuals) by snp in one file?
Thanks!
java -jar GenomeAnalysisTK.jar \
-R refs.fa \
-T HaplotypeCaller \
-I bams.list \
-- ERC GVCF \
-o snps.g.vcf
java -jar GenomeAnalysisTK.jar \
-T UnifiedGenotyper \
-R refs.fa -I bams.list \
--genotyping_mode DISCOVERY \
-o snps.raw.vcf -stand_call_conf 25.0 -stand_emit_conf 10.0