Ok, that's my impression as well. Thanks Pierre!!
• 0 views
•
link
Hi everyone!
I was wondering if GATK (specifically version 4.5.0) outputs partial genotypes, where only one allele in a diploid individual is called (e.g., ./0, 0/., ./1, 1/.? I ask this because I have super low-coverage WGS data, and the expectation based on average sample read depths is that only one allele should be called.
Log in to answer this question.
https://gatk.broadinstitute.org/hc/en-us/articles/6012243429531-GenotypeGVCFs-and-the-death-of-the-dot-obsolete-as-of-GATK-4-6-0-0
Thanks for sharing! :-)
I was confused because I do have missing genotypes coded as
./.instead of0/0in my VCF. Yet, I didn't find any partial genotypes, so I was curious if the program could even output that for a diploid. Based on the article, it seems like no.