Losing ~60% of variants from sequence data after phasing with SHAPEIT
0
0
Entering edit mode
3.2 years ago
kianalee ▴ 10

I have a sequenced locus on chromosome 9 that I am attempting to phase with SHAPEIT. My original vcf file has1617 variants. After phasing with SHAPEIT I am left with only 513 variants. I made sure to check the quality of my vcf with checkVCF and found no duplicates or reference mismatches.

I use the following code to run SHAPEIT:

shapeit4.2 --input my_data.vcf.gz --map chr9.b37.gmap.gz --region 9 --reference 1000g_phase3_nomulti_allelic/ALL.chr9.phase3_shapeit2_mvncall_integrated_v5b.20130502.genotypes.vcf.gz --thread 8 --log shapeit_chr9.log --output my_data_phased_SHAPEIT.vcf --mcmc-iterations 10b,1p,1b,1p,1b,1p,1b,1p,10m --pbwt-depth 8 &

Is there a way to minimize the loss of variants?

phasing SHAPEIT • 745 views
ADD COMMENT
0
Entering edit mode

It's likely that it's because these variants aren't in your reference panel - can you check for how many of your 1617 variants are in the reference panel?

ADD REPLY

Login before adding your answer.

Traffic: 2183 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