googlegenomics / gcp-variant-transforms

GCP Variant Transforms
Apache License 2.0
134 stars 55 forks source link

The MISSING_GENOTYPE_VALUE value of -1 can conflict with genotype values that are in a VCF (scalar vs. list) #703

Open adaykin opened 2 years ago

adaykin commented 2 years ago

When running a bq_to_vcf command on data that has list values for genotype, e.g. [-1, -1] the densify_variants._densify_variants method will create a VariantCall object that defaults to using the MISSING_GENOTYPE_VALUE, which is set to -1. In vcf_parser.VariantCall.lt the genotypes are compared, so comparing a scalar value to a list value results in an exception being thrown.