欢迎关注"生信修炼手册"!
本篇主要介绍index
, view
, query
, sort
, reheader
这五个命令。
1. index
index
命令用于对VCF文件建立索引,要求输入的VCF文件必须是使用bgzip
压缩之后的文件,支持.csi
和.tbi
两种索引,默认情况下建立的索引是.csi
格式, 用法如下
bgzip view.vcf
bcftools index view.vcf.gz
运行成功后,会生成索引文件view.vcf.gz.csi
。如果需要建立.tbi
格式的索引,用法如下
bcftools index -t view.vcf.gz
tbi索引文件为view.vcf.gz.tbi
。
2. view
view
命令可以用于VCF和BCF格式的转换,用法如下
bcftools view view.vcf.gz -O u -o view.bcf
-O
参数指定输出文件的类型,b
代表压缩后的BCF文件,u<