vcffilter 工具bug以及解决办法

1，使用说明：

usage: vcffilter [options] <vcf file>

options:

-f, --info-filter     specifies a filter to apply to the info fields of records, removes alleles which do not pass the filter

-g, --genotype-filter specifies a filter to apply to the genotype fields of records

-s, --filter-sites    filter entire records, not just alleles

-t, --tag-pass        tag vcf records as positively filtered with this tag, print all records

-F, --tag-fail        tag vcf records as negatively filtered with this tag, print all records

-A, --append-filter   append the existing filter tag, don't just replace it

-a, --allele-tag      apply -t on a per-allele basis.  adds or sets the corresponding INFO field tag

-v, --invert          inverts the filter, e.g. grep -v

-o, --or              use logical OR instead of AND to combine filters

-r, --region          specify a region on which to target the filtering, requires a BGZF compressed file which has been indexed with tabix.  any number of regions may be specified.

Filter the specified vcf file using the set of filters. Filters are specified in the form "<ID> <operator> <value>:

-f "DP > 10"  # for info fields  -g "GT = 1|1" # for genotype fields  -f "CpG"  # for 'flag' fields

Operators can be any of: =, !, <, >, |, &

Any number of filters may be specified.  They are combined via logical AND unless --or is specified on the command line.  Obtain logical negation through the use of parentheses, e.g. ! "( DP = 10 )"

For convenience, you can specify "QUAL" to refer to the quality of the site, even though it does not appear in the INFO fields.

2，举例:

以上是说明，我举个例子，假如你对C18.vcf文件筛选的要求是：质量要大于30，深度要大于10，类型为SNP，那么执行如下命令：

vcffilter -f "QUAL > 30 & DP > 10 & TYPE = snp" C18.vcf > C18.filtered.vcf

会得到经过筛选过的vcf文件。

by freemao

FAFU.

free_mao@qq.com