Convert SHAPEIT haps file to vcf

Convert SHAPEIT haps file to vcf

0

I’ve phased some sequencing using SHAPEIT, producing two files:

125QiPSC.haps
125QiPSC.log

I would now like to subset the file for my SNPs of interest. To achieve this, I’m first attempting to convert it to a vcf. The SHAPEIT manual suggests this can be done by –convert

However, when I try running this:

shapeit -convert 
--input-haps "$DIR"/shapeit/125QiPSC 
--output-vcf "$DIR"/shapeit/125QiPSC.vcf

I get the error message

Converter mode : unrecognised option '--output-vcf'

Am I doing something wrong?


convert


vcf


conda


shapeit

• 19 views

Read more here: Source link