This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Phasing a partly-phased VCF

Hello, I am trying to phase a VCF file for a non-model organism that has no reference panel. I have a VCF file of genotypes for 30 samples from various populations of a species, which I phased using read-based phasing with Whatshap. This managed to phase about 75% of my heterozygous sites, but a downstream analyses requires all heterozygous sites to be phased. Is there a method that can phase just the unphased sites/genotypes in my VCF, without disrupting the existing phasing? From the documentation, it seems like Beagle would treat the entire VCF as unphased and undo my read-based phasing. Thank you!

heterozygous whatshap sites phasing

Unfortunately if you don't have a reference panel it's not really possible to phase only 30 samples - the results would be very unreliable with that few samples.

0 answers

No answers yet.

Log in to answer this question.