Question: VCFtools filter - genotype quality score changes
0
gravatar for AW
3.1 years ago by
AW350
United Kingdom
AW350 wrote:

Hi,

I ran the following filtering command

vcftools --vcf test.vcf --out test.test.vcf --remove-indels --remove-filtered-all --max-alleles 2 --max-alleles 2 --recode --recode-INFO-all

However, when I compare the new test.test.vcf compared to the original test.vcf the GQ scores of some snps have changed.

eg

test.vcf

scaffold11  389 .   T   C   .   PASS    ADP=22;WT=0;HET=0;HOM=1;NC=0    GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR    1/1:123:22:22:0:22:100%:4.7526E-13:0:41:0:0:17:5

test.test.vcf

scaffold11  389 .   T   C   .   PASS    ADP=22;WT=0;HET=0;HOM=1;NC=0    GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR    1/1:99:22:22:0:22:100%:4.7526E-13:0:41:0:0:17:5

The filter did not remove any snps, therefore why has the GQ score been downgraded?

Parameters as interpreted:

    --vcf test.vcf
    --recode-INFO-all
    --max-alleles 2
    --out test.test.vcf
    --recode
    --remove-filtered-all
    --remove-indels

After filtering, kept 1 out of 1 Individuals
Outputting VCF file...
After filtering, kept 4717895 out of a possible 4717895 Sites
Run Time = 156.00 seconds

Thanks!

A

vcftools gq genotype quality • 1.3k views
ADD COMMENTlink modified 3.1 years ago by Brice Sarver3.1k • written 3.1 years ago by AW350
0
gravatar for Brice Sarver
3.1 years ago by
Brice Sarver3.1k
United States
Brice Sarver3.1k wrote:

Genotype qualities are calculated by setting the phred-scaled likelihood of the best genotype to zero and then calculating the difference from the next most likely. My guess is that your variant caller has a higher max GQ than is allowed by VCFtools (which might be capped at 99 by default) and VCFtools is placing it at the maximum. Not really a ton of difference between 99 and 123 - still quite confident in your call relative to the others.

ADD COMMENTlink modified 3.1 years ago • written 3.1 years ago by Brice Sarver3.1k
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: 1696 users visited in the last hour