Hi everyone,
I'm trying to calculate heterozigosity for one individual genome using vcftools. For that, I need to have GT information in the FORMAT field for all sites. However, using the command "samtools mpileup -C50 -uDgf ref.fa output.sorted.bam | bcftools view -cg - > output.vcf" in samtools, only adds the GT information for variable sites relative to the reference, which is apparently a problem for vcftools. Anyone has any idea how to tell samtools to include the GT field for all sites? I thought the -g option in both mpilup and bcfview would take care of it, but apparently not.
Thanks!
I'm trying to calculate heterozigosity for one individual genome using vcftools. For that, I need to have GT information in the FORMAT field for all sites. However, using the command "samtools mpileup -C50 -uDgf ref.fa output.sorted.bam | bcftools view -cg - > output.vcf" in samtools, only adds the GT information for variable sites relative to the reference, which is apparently a problem for vcftools. Anyone has any idea how to tell samtools to include the GT field for all sites? I thought the -g option in both mpilup and bcfview would take care of it, but apparently not.
Thanks!