vcf filtering based on % of heterozygosity
1
0
Entering edit mode
8.8 years ago
Matt • 0

Hi all

Is it possible to use vcftools to filter vcf files based on % of heterozygosity ?!

best

M

SNP vcf vcftools hapmap • 4.0k views
ADD COMMENT
0
Entering edit mode
8.8 years ago
DG 7.3k

You can use vcflib for arbitrary filtering based on values in the INFO field. Once you download and install it there are specific examples of filtering based on criteria from the info field here.

ADD COMMENT
0
Entering edit mode

I don't think you can't do that using vcffilter

'GT=1|1' >= 0.05
ADD REPLY
0
Entering edit mode

If you're looking to do something like that you'll need to use a tool that already calculated the number of heterozygotes in your multi-sample VCF, write a wrapper script around another tool, or use the vcftools Perl API. I had assumed you already had your % heterozygosity calculated using an annotation tool and had it in your VCF in the INFO field. Otherwise you need to calculate this before you can filter on it.

ADD REPLY

Login before adding your answer.

Traffic: 2173 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