bcftools学习笔记

软件简介

bcftools 是一个用于操作和解决 VCF/BCF 文件的软件工具集。它是 samtools 工具集的一部分，用于比照对后的 BAM 文件进行 SNP 和 Indel 的检测。bcftools 能够用于 SNP 和 Indel 的过滤、正文、统计和可视化等操作。

装置形式

bcftools 能够通过以下两种形式进行装置：

1. 通过源代码装置

从 bcftools 的官方网站 http://www.htslib.org/download/ 下载最新的源代码包。解压后进入解压后的目录，运行以下命令进行编译和装置：

./configuremakesudo make install

2. 应用包管理器装置

bcftools 能够通过包管理器进行装置，如应用 Ubuntu 零碎能够通过以下命令进行装置：

sudo apt-get install bcftools

应用办法

1. 查看 VCF/BCF 文件信息

bcftools view 命令能够用于查看 VCF/BCF 文件的信息，如查看文件头信息：

bcftools view -h file.vcf

查看文件中的 SNP 和 Indel 信息：

bcftools view -v snps file.vcfbcftools view -v indels file.vcf

2. 过滤 SNP 和 Indel

bcftools filter 命令能够用于过滤 SNP 和 Indel，如过滤掉 quality 值低于 20 的 SNP：

bcftools filter -i '%QUAL>20' file.vcf > filtered.vcf

3. 统计 SNP 和 Indel

bcftools stats 命令能够用于统计 SNP 和 Indel 的信息，如统计文件中 SNP 和 Indel 的数量：

bcftools stats file.vcf > stats.txt

4. 正文 SNP 和 Indel

bcftools annotate 命令能够用于正文 SNP 和 Indel，如正文 dbSNP 中的 rsID：

bcftools annotate -a dbsnp.vcf -c ID file.vcf > annotated.vcf

5. 可视化 SNP 和 Indel

bcftools view 命令能够用于可视化 SNP 和 Indel，如将文件中的 SNP 和 Indel 可视化为 HTML 文件：

bcftools view -H -v snps file.vcf | vcfutils.pl varFilter -d 5 > filtered.vcfbcftools view -H -v indels file.vcf | vcfutils.pl varFilter -d 5 >> filtered.vcfbcftools stats filtered.vcf > stats.txtplot-vcfstats -p plots/ stats.txt

外围性能

  • bcftools view

用于查看VCF文件的内容，包含样本信息、位点信息、正文信息等

  • bcftools query

用于提取VCF文件中的信息，能够依据须要抉择提取特定的字段或样本信息

  • bcftools stats

用于生成VCF文件的统计信息，包含变异类型、变异频率、基因型品质等

  • bcftools merge

用于合并多个VCF文件

  • bcftools norm

用于对VCF文件进行规范化解决，包含拆分多等位位点、合并相邻位点等

  • bcftools annotate

用于给VCF文件增加正文信息

参数解释

bcftools 命令中罕用参数：

  • -f: 用于指定参考基因组文件。
  • -i: 用于进行过滤，只保留符合条件的位点或样本信息。
  • -c: 用于对VCF文件进行压缩，进步文件的存储和传输效率。
  • -g: 用于将VCF文件转换为BGZF格局，进步文件的存储和传输效率。
  • -O: 用于指定输入文件格式，可选的格局有VCF、BCF、BAM等。
  • -s: 用于指定须要解决的样本信息。
  • -R: 用于指定须要解决的区域。

实战演练

首先利用bcftools软件将vcf格局生成gz格局和inde索引格局，因为bcftools的输出文件有这两项，须要用到“-Oz”和“index”命令，具体如下：

bcftools view ExAC.vcf -Oz -o ExAC.vcf.gzbcftools index ExAC.vcf.gz

拆分、筛选、合并

应用bcftools能够不便地对VCF.gz文件进行拆分、筛选和合并，以下是具体的命令：

1. 拆分VCF.gz文件

应用bcftools的split命令能够将一个VCF.gz文件拆分为多个小的VCF.gz文件，每个文件蕴含指定数量的记录。例如，将一个VCF.gz文件依照每个文件蕴含1000个记录的规范拆分，能够应用以下命令：

bcftools split -O z -o out.vcf.gz -l 1000 file.vcf.gz

其中，-O z选项指定输入文件的格局为.gz压缩的VCF格局，-o选项指定输入文件名的前缀，%02d示意输入文件名中的序号，-l选项指定每个输入文件蕴含的记录数量。

2. 筛选VCF.gz文件

应用bcftools view命令能够依据须要筛选VCF.gz文件中的记录。例如，筛选出染色体为chr1且地位在1000000到2000000之间的记录，能够应用以下命令：

bcftools view -O z -o output.vcf.gz file.vcf.gz chr1:1000000-2000000

其中，-O z选项指定输入文件的格局为.gz压缩的VCF格局，-o选项指定输入文件名为output.vcf.gz，chr1:1000000-2000000示意须要筛选的区间。

3. 合并VCF.gz文件

应用bcftools的merge命令能够将多个VCF.gz文件合并为一个文件，例如：

bcftools merge -O z -o merged.vcf.gz file1.vcf.gz file2.vcf.gz file3.vcf.gz

其中，-O z选项指定输入文件的格局为.gz压缩的VCF格局，-o选项指定输入文件名为merged.vcf.gz，前面跟着须要合并的多个VCF.gz文件名。

提取指定样本的变异信息

筹备样本ID文件，这里命名为sample.txt，一个样本一行，而后应用下列参数批改：

bcftools view -S sample.txt ordata.vcf.gz -Ov > sampleout.vcf

批改vcf文件的样本名称

应用bcftools的 rehead指令，sample.file 是指定要替换的样本名（蕴含两列，别离是新名称和旧名称）

bcftools reheader -s sample.file old.vcf -o new.sample.vcf

批改vcf.gz染色体编号形式

nohup bcftools annotate --rename-chrs chr_id.txt xx.vcf.gz -Oz -o xxx.vcf.gz --threads 23 &

统计vcf中样品和位点数量

办法：bcftools +counts xxx.vcf

$ bcftools +counts Chr.recode.vcfNumber of samples: 2232Number of SNPs:    698763Number of INDELs:  44021Number of MNPs:    0Number of others:  0Number of sites:   742784

提取vcf文件中的GT信息

bcftools annotate -x 'fmt' -O z -o out.vcf.gz in.vcf.gz

从VCF中提取区域的数据

每次须要批改的是染色体和地位信息，用regions参数传递，命令如下：

bcftools filter xx.vcf.gz --regions chr1A:5000-50900 > chr1A_5000-5090.vcf
参考资料：
https://www.jianshu.com/p/99895c7338b2https://www.jianshu.com/p/266c55c87978http://samtools.github.io/bcftools/bcftools.htmlhttps://www.biostars.org/p/95013/

本文由mdnice多平台公布