This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Phasing without a reference panel

Hi all,

I'm trying to make my own cohort from genotypes like HapMap project. I would like to know if there are any phasing tools which contain options for phasing genotypes when there's no reference panel?

Thank you.

phasing genotype hapmap population
• 5,496 views • link
updated by Wietje Wietje • written by chadapohn.chaosrikul • chadapohn.chaosrikul

1 answer

Beagle 4.0 offers the option to impute and phase without a reference panel

You should be able to run it with: java -jar ./beagle gt=[your.file.vcf] phase-its=5 nthreads=1 out=[your.output]

You can speed up the process with more threads and improve phasing accuracy by augmenting the phase-its number.

Documentation can be found here: https://faculty.washington.edu/browning/beagle/beagle.29Sep14.pdf

• 0 views • link
written by Wietje Wietje

Thank you for your help! I'd like to know if I used SHAPEIT without a reference (by not specify --input-ref option) Does it use the default reference e.g 1000GP or use MCMC algorithm to estimate haplotypes instead?

• 0 views • link
written by chadapohn.chaosrikul • chadapohn.chaosrikul

As far as I know, you have to invoke the reference panel option with --input-ref to make use of such a panel - otherwise SHAPEIT will phase without this kind of background information. If you want to use the 1000GP data for your phasing step, you'd have to download it as described on the SHAPEIT page (also check these instructions for further parameters and input/output options:

https://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html#prephasing

• 0 views • link
written by Wietje Wietje

I've already tried both Beagle and SHAPEIT. Thank you for your suggestions.

• 0 views • link
written by chadapohn.chaosrikul • chadapohn.chaosrikul

Log in to answer this question.

More posts like this