在群体遗传学相关的系数计算这一块，vcftools是个很好的工具，操作简单，可以从变异检测的vcf出发，直接计算Pi系数，TajimaD系数以及Fst系数，还是比较方便。

安装

在老版本不支持vcf4.2的情况下，可以选择安装新版vcftools。
下载地址为：https://sourceforge.net/projects/vcftools/files/

以最新版本0.1.13版本为例，首先我们需要下载源码安装包：

$ wget https://sourceforge.net/projects/vcftools/files/vcftools_0.1.13.tar.gz/download -o vcftools_0.1.13.tar.gz
$ tar zxvf vcftools_0.1.13.tar.gz
$ cd vcftools_0.1.13
$ make

这个安装包里面不包含configure文件，也就意味着不需要指定安装目录，软件会自动安装在当前目录，我们按上述代码解压之后直接make就好。

报错信息

安装过程可能会碰到错误提示，如

vcf_file.cpp:(.text+0xe72): undefined reference to `gzbuffer'
collect2: ld returned 1 exit status
make[1]: *** [vcftools] Error 1
make[1]: Leaving directory `/home/Public/Packages/vcftools_0.1.10/cpp'
/bin/sh: 2: cd: can't cd to perl
make: *** [install] Error 2

在出现这个错误的时候，我去geogle了一下（geogle到的网页是https://www.biostars.org/p/61236/），发现是zlib库版本不高，需要安装高版本的zlib，或者说你安装了高版本的zlib库，但是不在默认的库搜索路径，如/usr/local/lib、/usr/lib64目录下。

你可以对cpp/Makefile文件进行修改，指定LIB的库路径，-L 后接的是你高版本zlib的安装目录。
# Included libraries (zlib) #LIB = -lz LIB = -lz -I /usr/include/ -L /home/Public/Packages/zlib-1.2.7/

在修改完这个代码后，再次make即可安装成功。
vcftools的参数含义可以看vcftools参考文献，或者见博客.

vcftools的使用

在计算Pi系数的时候，代码如下：

vcftools --vcf all.gatk.fliter.vcf --remove-filtered-geno-all --site-pi --out watermelon.pi

运行信息如下：

VCFtools - v0.1.13
(C) Adam Auton and Anthony Marcketta 2009Parameters as interpreted:--vcf all.gatk.fliter.vcf--out watermelon.pi--site-pi--remove-filtered-geno-allAfter filtering, kept 19 out of 19 Individuals
Outputting Per-Site Nucleotide Diversity Statistics...
After filtering, kept 8079873 out of a possible 8079873 Sites
Run Time = 195.00 seconds

在计算Fst系数时候，代码如下：

vcftools --vcf all.gatk.fliter.vcf --remove-filtered-geno-all --weir-fst-pop population_1.txt --weir-fst-pop population_2.txt --out watermelon.pi

vcftools安装碰到的坑相关推荐

centos yum安装_centos7上yum安装碰到的坑
估计是之前更改过这个虚机的配置,故现在正常步骤安装完yum之后,无法使用, 具体表现在执行 yum makecache的时候报错,现在no module name yum 结果发现,无论怎么修改 u ...
bbqsql安装使用踩坑总结
bbqsql是一个知名的盲注SQL开源工具,可以用来做一些简单传统的sql注入的渗透测试或者靶场模拟等,但是由于其版本较老,使用的还是python2,且很多库都是很旧的,接近9年多没有维护更新,因此如 ...
Pytorch可视化工具tensorboardX（安装不踩坑）
安装方法(如何不踩坑,已修改): 首先默认你已经安装好pytorch(如果木有,为什么点进来呢) 安装tensorboardX pip install tensorboardX 安装tensorflo ...
说说初用 Mock 工具测试碰到的坑
说说初用 Mock 工具测试碰到的坑我是一个在校实习生,作为一个程序猿,是个菜鸟中战斗机!对于测试,只写过一点点简单到不能再简单了的 Junit 单元测试的例子(因为当时这足以应付学校课程的内容与要 ...
cmd命令安装composer踩坑
cmd命令安装composer踩坑很多童鞋在依照composer官网官网教程用命令行安装composer时回踩到以下坑,却不知道如何解决在此我分享下自己的经验,写个不好请多海涵打开命令行执行 ...
centos7安装mysql日志空白_centos7安装Mysql爬坑记录 - G
centos7安装Mysql爬坑记录查看是否已安装使用下列命令查看是否已经安装过mysql/mariadb/PostgreSQL 如果未安装,不返回任何结果(ECS的centos镜像默认未安装任何 ...
GPCC安装以及踩坑经历
gpcc安装以及踩坑经历官方下载地址文档 https://network.pivotal.io/products/pivotal-gpdb#/releases/29190 安装开始之前 chown ...
学习使用Visual studio 时碰到的坑
学习使用Visual studio 时碰到的坑文章目录学习使用Visual studio 时碰到的坑异常1:OpenCV中出现"Microsoft C++ 异常: cv::Except ...
前端项目开发中碰到的坑、移动端兼容性问题
前端开发中碰到的坑.移动端兼容性问题 1.ios兼容input光标高度问题:input输入光标,在安卓手机上显示没有问题,但是在苹果手机上当点击输入的时候,光标的高度和父盒子的高度一样. **出现问 ...
vant vant-list碰到的坑
vant vant-list碰到的坑最近新写的项目都是碰见了这个问题 onload在加载时只触发一次,还是在刚进页面是触发一次,页面向下滚动时,onload并不加载 (本人所尝试的解决方案) 我在d ...

vcftools安装碰到的坑

安装

报错信息

vcftools的使用

vcftools安装碰到的坑相关推荐

最新文章

热门文章