madprime / cgivar2gvcf

Lossy conversion of Complete Genomics var file to VCF
MIT License
5 stars 3 forks source link

Catch uncaught homozygous no-calls #7

Closed madprime closed 6 years ago

madprime commented 6 years ago

The algorithm is imperfect and falls back into making "no-calls"; this results in some complex var file regions being processed in the var_only=True condition but becoming a homozygous no-call, creating erroneous output that violates the "variant only" expectation.