Convert VCF + LDhelmet Recombination Rates to EIGENSTRAT format
0
0
Entering edit mode
17 months ago
Jautis ▴ 320

Hello, I have a vcf file of called variants (GATK) and genotypes per individual and a text file produced from LDhelmet which contains recombination rates for windows across the genome (4 columns: chromosome, start, end, and recombination rate). I'm looking to convert this to the eigenstrat format.

I feel like I could write something ugly to (i) convert bp to cM based on the linkage information and (ii) reformat the vcf file (maybe using Mega2?), but I also think this has got to be a problem someone else has already solved much more efficiently. Thanks, and I'd appreciate any suggestions about how to best approach this problem!

vcf linkage SNP genome • 450 views
ADD COMMENT

Login before adding your answer.

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