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