0
点赞
收藏
分享

微信扫一扫

PSMC分析有效群体大小(Ne)

1 软件准备

bwa samtools bcftools picard  psmc

2 数据准备

二倍体物种(该分析是结合杂合度来分析的)!!

Ref_genome.fa: 组装好的核基因参考序列(个体A);

A1_R1.fa.gz,A1_R2.fa.gz: 二代重测序clean_reads(个体B)【注:必须两个个体A与B】

3 PSMC分析

##1 index reference by bwa对参考建立索引

##2 mapping 将readsmapping到参考序列

##3 将SAM格式转换为BAM   

##4 对BAM格式进行排序

##5 移除重复reads

##6 建立索引

##7 生成psmc的输入文件

##8  运行PSMC分析

##9 作图

##10 结果展示

将会生成一个.eps文件,用AI或PS打开即可。


参考来源:


PSMC软件分析群体历史有效群体大小步骤  https://blog.csdn.net/zaprily/article/details/108764219

PSMC分析流程  http://www.wuchangsong.com/?p=555

mpileup命令简介 https://blog.csdn.net/u013553061/article/details/53293302

举报

相关推荐

0 条评论