Question: Merging vcf files ?
0
gravatar for Picasa
3.9 years ago by
Picasa530
Picasa530 wrote:

I have 8 individuals from whom I run variant analysis (with GATK). So I have now 8 .vcf files.

These 8 individuals can be splitted into 2 groups and that's what I wanted to know.

Is there a way to merge the 8 .vcf files and to perform a kind of GWAS analysis to discriminate my samples ? Do you have any recommendation for a method ?

Thanks.

vcf • 2.2k views
ADD COMMENTlink modified 15 months ago by jaybee40 • written 3.9 years ago by Picasa530

If you can clarify more. Do you perform non somatic variant call or variant call having both somatic and germline? Check this thread

Vcftool is there do so such.

Just check this post carefully and see before performing the merge.

ADD REPLYlink modified 3.9 years ago • written 3.9 years ago by ivivek_ngs4.9k

Hi,

How is it important between non somatic and with somatic ??

ADD REPLYlink written 3.9 years ago by Picasa530

It depends on your biological reasoning about what you are actually interested in, if you are trying to find mutations that is driving the disease specfic to your patients (exclusive or common) against control data of a family member who does not have the disease then it might help you in defining the mutations that are germline/somatic/ de novo. Then you can always associate them with any standard annotation softwares like ANNOVAR, SNPeff or VEP.

ADD REPLYlink written 3.9 years ago by ivivek_ngs4.9k

vcftools has vcc-merge module. It might help you through with the merging multiple vcf files.

ADD REPLYlink modified 3.7 years ago • written 3.7 years ago by morovatunc460
0
gravatar for jaybee
15 months ago by
jaybee40
South Korea
jaybee40 wrote:

You can try /PATH/to/bcftools merge Home/data/*vcf.gz -Oz -o Merged.vcf.gz

[considering that your files A.vcf.gz , B.vcf.gz and so on are in the Home directory, under the folder data]

[If your bcftools is installed in the usr/bin directory then simply use:

/usr/bin/bcftools merge Home/data/*vcf.gz -Oz -o Merged.vcf.gz

This should work fine!

ADD COMMENTlink written 15 months ago by jaybee40
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: 1109 users visited in the last hour