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

Hello,

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!

Solène

tabix pooling exome • 298 views
ADD COMMENTlink modified 11 months ago by Nicolas Rosewick9.2k • written 12 months ago by solene.conrad10
1
gravatar for Nicolas Rosewick
11 months ago by
Belgium, Brussels
Nicolas Rosewick9.2k 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 (https://www.ensembl.org/info/docs/tools/vep/index.html) 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 11 months ago by Nicolas Rosewick9.2k
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: 1342 users visited in the last hour