Question: Distinguish alleles of a SNP in trio VCF files
1
gravatar for charlesberkn
4.7 years ago by
charlesberkn20
United States
charlesberkn20 wrote:

Hi All,

I have three VCF files, one for father, one for mother, and one for son. I want to distinguish each allele of a SNP in son's VCF file from father or mother, just like a Mendelian analysis of inheritance.  Are there any tools or proper way to achieve this? Thanks.

snp • 1.7k views
ADD COMMENTlink modified 4.7 years ago • written 4.7 years ago by charlesberkn20

Thank you all guys, I will try your suggestions. 

ADD REPLYlink written 4.7 years ago by charlesberkn20

Hi All,

I have tried three ways suggested above, but based on my personal experience, they are not easy to use. I just want to know if there are any other tools available. Thanks.

ADD REPLYlink modified 6 months ago by RamRS27k • written 4.7 years ago by charlesberkn20

what would be a tool easy to use for your personal experience ? what's missing in the other tools ?

ADD REPLYlink written 4.7 years ago by Pierre Lindenbaum128k

For example, I used CombineVariants to combine vcf files, but it cannot handle different records in the CHROM field as I posted here.

ADD REPLYlink modified 6 months ago by RamRS27k • written 4.7 years ago by charlesberkn20
0
gravatar for reza.jabal
4.7 years ago by
reza.jabal370
New York, USA
reza.jabal370 wrote:

Hi,

1) Convert your VCF files to annovar files (see here)

2) Use Microsoft excel VLOOKUP function to match segregation consistent variants.

ADD COMMENTlink written 4.7 years ago by reza.jabal370
0
gravatar for cpad0112
4.7 years ago by
cpad011213k
India
cpad011213k wrote:

If they are in single VCF file, I would have suggested, vcf-contrast or GATK select variant walker.
 

ADD COMMENTlink written 4.7 years ago by cpad011213k
0
gravatar for Pierre Lindenbaum
4.7 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum128k wrote:

merge with gatk CombineVariants https://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_tools_walkers_variantutils_CombineVariants.php

you might use my tool https://github.com/lindenb/jvarkit/wiki/FixVcfMissingGenotypes to find wether the uncalled sites were REF/REF or ./.

extract the mendelian incompatibilities with PhaseByTransmission

ADD COMMENTlink written 4.7 years ago by Pierre Lindenbaum128k

Pierre have you got a piece of code that query intersections in ANNOVAR files?

ADD REPLYlink written 4.6 years ago by reza.jabal370
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: 1280 users visited in the last hour