利用群体遗传数据估计基因组上重组率

本文介绍了如何利用FastEPRR软件来估计基因组上的重组率。首先,需要将数据转换为VCF格式,可以通过TASSEL或Mega2完成。接着,使用Beagle4.1进行未定位基因型的插补。然后,使用vcftools按染色体分割VCF文件。最后,在R环境中安装并运行FastEPRR,通过三个步骤计算单个染色体的重组率。若要计算整个基因组,需对所有染色体重复此过程,可以在Linux下编写脚本来自动化此任务。

开发板推荐:天空星STM32F407VET6开发板

超高性价比 STM32主控 | 超高主频 | 一板兼容百芯 | 比赛神器 | 沉金彩色丝印

计算所用软件:FastEPRR,需要熟悉和掌握的软件,PLINK,TASSEL, vcftools, Beagle4.1, R

准备工作:
1. 准备输入文件,phased genotype file with VCF format
2. 如何实现 phased genotype file? 假如我们手中的是unphased的 vcf 文件或者hapmap 或者plink 软件格式。
首先准备VCF 格式,可以采用Mega2软件或者TASSEL 进行转化,我的文件是plink 的格式,用TASSEL 转化(这个时候一定要设置最小等位基因频率进行过滤,否则产生的VCF文件有可能报错),如下:

./run_pipeline.pl -fork1 -plink -ped file.ped -map file.map -filterAlign -filterAlignMinFreq 0.05 -export file_vcf 

开发板推荐:天空星STM32F407VET6开发板

超高性价比 STM32主控 | 超高主频 | 一板兼容百芯 | 比赛神器 | 沉金彩色丝印

评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值