Closed mlin closed 4 years ago
It's possible for the genotype revision to leave allele in the pVCF without any (hard) calls. htslib has a convenient subroutine we can use to drop them before writing out the final row: https://github.com/samtools/htslib/blob/123bd90394e6bc49ebc869b5cc43a31076d9523c/htslib/vcfutils.h#L47
It's possible for the genotype revision to leave allele in the pVCF without any (hard) calls. htslib has a convenient subroutine we can use to drop them before writing out the final row: https://github.com/samtools/htslib/blob/123bd90394e6bc49ebc869b5cc43a31076d9523c/htslib/vcfutils.h#L47