whatsHap didn't phase all variants
0
0
Entering edit mode
5 months ago
s2004341 • 0

Hi there, I am using whatsHap to phase my filtered variants. The job was ran well and the output file generates normally. But when I was checking the contents of file, I found only a small proportion of SNPs has been phased. Anybody know why whatsHap didn't phase all SNPs, or is there some mechanism make this limit?

enter image description here

This screenshot was produced through "diff \ -biw \ bcftools_QUAL78_sample_1.vcf \ Phased_Q78_Horti_sample_1.vcf" which represents the phases part ("/" replaced by "|")

While other SNPs were still unphased.

Many thanks, Zitian

whatsHap phasing • 301 views
ADD COMMENT
1
Entering edit mode

WhatsHap can only phase variants if they are covered by reads which overlap with a different variant. So my guess is that the unphased variants didn't;t have that.

ADD REPLY
0
Entering edit mode

Thank you, I'll check that!

ADD REPLY
1
Entering edit mode

Your screenshot is rather unclear (hint: don't post screenshots but just paste the text and format it nicely with the 101010 button), but perhaps it is worth mentioning that phasing doesn't make a lot of sense for all these homozygous variants. You need heterozygous variants. In addition, I see most (if not all?) positions have a PS field in their genotype, and hence, are actually in a phase block.

ADD REPLY

Login before adding your answer.

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