Question: Creating a set of rare variants shared by a child and a parents with VCF and Tabix
gravatar for solene.conrad
7 weeks ago by
solene.conrad10 wrote:


I have two single WES from a child and his father. I would like to create a set a of rare variants inherited from the father. Afterwards, my goal is to compare if these rare variants still be found in a WES where I have pooled 2 fathers, to tests the sensibility of pooling.

I have been told that I could use Tabix but I don't know where to find the information on allele frequency in general population from GnomAD, to see what variants are rares? I could work on FastQ, cram and VCF formats.

Thank you!


tabix pooling exome • 145 views
ADD COMMENTlink modified 6 weeks ago by Nicolas Rosewick8.6k • written 7 weeks ago by solene.conrad10
gravatar for Nicolas Rosewick
6 weeks ago by
Belgium, Brussels
Nicolas Rosewick8.6k wrote:

Some hints :

First filter only variants with a gnomAD AF under a certain threshold to be considered as rare ( < 1% ; < 0.1% , up to you) in the father and son vcf. To do that you could first annotate your VCFs using VEP ( and then filter with filter_vep or bcftools

Then intersect father and son vcf using bcftools isec -n +2 father_rare.vcf son_rare.vcf

ADD COMMENTlink written 6 weeks ago by Nicolas Rosewick8.6k
Please log in to add an answer.


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