1、群體區(qū)分
將每個(gè)群體中所包含的個(gè)體一行一個(gè)，生成popA，popB文件
2、vcftools進(jìn)行對snp文件進(jìn)行劃窗計(jì)算fst和pi
窗口大小設(shè)置為1 Mbp，步長為0.1 Mbp
fst
vcftools?--vcf populations.snp.vcf \
?????????--weir-fst-pop?popA?\
?????????--weir-fst-pop?popB?\
?????????--out?popA_vs_popB?\
?????????--fst-window-size?1000000?\
?????????--fst-window-step?100000
pi
vcftools?\
??--vcf?populations.snp.vcf?\
??--window-pi?1000000?\
??--window-pi-step?100000?\
??--out?popA?\
??--keep?popA
vcftools?\
??--vcf?populations.snp.vcf?\
??--window-pi?1000000?\
??--window-pi-step?100000?\
??--out?popB?
??--keep?popB
繪圖：參考
