find position of snps in cM, given recombination map from ldhelmet
0
0
Entering edit mode
13 months ago
Jautis ▴ 300

Hello, I have a recombination map calculated from LDhelmet which has estimates of recombination rates for windows of the genome (in recombination events per generation per bp) and a separate file of SNP locations which I would like to transform from base number to location in centiMorgans. Is there an easy way to do this computationally? The best strategy I've been able to devise is, for each site, taking the sum of previous windows and then the proportion of the current window. However, that seems very inefficient when given hundreds of thousands of SNPs. Any suggestions for a way to calculate this (or program out there) that would be more efficient?

##Recombination rate data
 start     end      rate/bp     length 
 101       400      0.004       300
 401       900      0.0001      500

## SNPs
base pair
 800

# So for this SNP, the position would be 300*0.004 + (400/500)*0.0001
recombination genome physical position • 351 views
ADD COMMENT

Login before adding your answer.

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