odelaneau / shapeit5

Segmented HAPlotype Estimation and Imputation Tool
https://odelaneau.github.io/shapeit5/
MIT License
66 stars 11 forks source link

Feature request: VCF output #105

Open pettyalex opened 3 months ago

pettyalex commented 3 months ago

My group still works with many tools that cannot read in BCF and need .vcf.gz still.

Would it be possible to allow outputting VCF? It seems like right now, the only reason that we cannot output VCF is because of bugs in xcftools, specifically because https://github.com/odelaneau/xcftools/commit/394bbaff4f90ea5971e018647171c534d53f9571 is not released in 5.1.1 and there's a similar bug with the line

        } else if (hts_fname.size() > 5 && hts_fname.substr(hts_fname.size()-5) == "vcf.gz") {

because ".vcf.gz" is 6 characters, not 5.

Thoughts on VCF output?