Question: Meaning Of -A In Vcftools Vcfutils.Pl Varfilter?
0
gravatar for Dan
6.1 years ago by
Dan500
Cambridge
Dan500 wrote:

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 • 5.6k views
ADD COMMENTlink modified 6.1 years ago by pd3340 • written 6.1 years ago by Dan500

Anyone know why -p is an option ;-)

ADD REPLYlink written 6.1 years ago by Dan500
1
gravatar for pd3
6.1 years ago by
pd3340
pd3340 wrote:

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 COMMENTlink written 6.1 years ago by pd3340
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2008 users visited in the last hour