Meaning Of -A In Vcftools Vcfutils.Pl Varfilter?
1
0
Entering edit mode
11.8 years ago
Dan ▴ 540

The command line says:

Usage:   vcfutils.pl varFilter [options] <in.vcf>

Options: -Q INT    minimum RMS mapping quality for SNPs [10]
         -d INT    minimum read depth [2]
         -D INT    maximum read depth [10000000]
         -a INT    minimum number of alternate bases [2]
         -w INT    SNP within INT bp around a gap to be filtered [3]
         -W INT    window size for filtering adjacent gaps [10]
         -1 FLOAT  min P-value for strand bias (given PV4) [0.0001]
         -2 FLOAT  min P-value for baseQ bias [1e-100]
         -3 FLOAT  min P-value for mapQ bias [0]
         -4 FLOAT  min P-value for end distance bias [0.0001]
         -e FLOAT  min P-value for HWE (plus F<0) [0.0001]
         -p        print filtered variants

Note: Some of the filters rely on annotations generated by SAMtools/BCFtools.

The description of -a seems counter intuitive, why would you want two alternative alleles as the default?

Is this better written "minimum number of alternate reads supporting each alternative allele [2]"?

That would at least make sense as an option to have.

P.S. anyone know which filters rely on annotations generated by SAMtools/BCFtools?

vcf vcftools • 9.9k views
ADD COMMENT
0
Entering edit mode

Anyone know why -p is an option ;-)

ADD REPLY
1
Entering edit mode
11.8 years ago
pd3 ▴ 350

The help message wording may be confusing, it should read, as you suggested, the "minimum number of alternate reads supporting the alternative allele".

The most useful annotations in v0.1.19 are RPB, MDV, QBD, QUAL and SnpGap from vcftools. The exact combination of annotations depends on the data, a bit of experimentation is needed.

ADD COMMENT

Login before adding your answer.

Traffic: 1816 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6