This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Sequence Haplotypes From 1Kg Project Variants

Greetings,

I need to create haplotype sequences (fastas) from the 1000 genomes data. I was wondering if anyone knows of a tool before I start writing my own?

I want to be able to pass in hg19 coordinates & 1kg.vcf & hg19.fasta and reconstruct all the haplotypes for that region.

1000genomes haplotype vcf
• 3,441 views • link
updated by Peixe Peixe • written by Zev.Kronenberg Zev.Kronenberg

2 answers

VCFtools should be able to do this:

cat ref.fa | vcf-consensus file.vcf.gz > out.fa

• 89 views • link
written by Adam Adam

Do you know if this supports multiple individuals? I don't want it to collapse the individuals into a consensus .

• 0 views • link
written by Zev.Kronenberg Zev.Kronenberg

I haven't used it for this, but it does have a --sample argument that sounds like it will do what you want.

• 0 views • link
written by matted matted

Take a look at Plink options for working with haplotypes here. I have been recently using it quite a lot for some haplotype calculations, and sounds like it might work for what you want, if you specify the haplotypes. But I don't know if sequences as large as fasta are gonna be accepted. Anyway looking is free, and the plink tutorial is really good! ;)

Cheers!

• 81 views • link
written by Peixe Peixe

Log in to answer this question.

More posts like this