Retain non-variable sites after phasing with shapeit4
0
0
Entering edit mode
21 months ago
pabe ▴ 30

Hi,

I noticed that after phasing with shapeit4, all homozygous reference sites present in my input (unphased) vcf were removed in the output vcf. Is there a way to retain these non-variable sites? I've searched the documentation but can't seem to find an option.

If not, I'm thinking to extract all the rows with . in the ALT column from the input vcf, then force all 0 / 0 calls to 0 | 0 but it seems kind of convoluted and I'm not sure how to then combine that with the phased vcf with only variable sites.

Thanks!

phasing vcf shapeit awk shapeit4 • 298 views
ADD COMMENT

Login before adding your answer.

Traffic: 2138 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6