Question: Merging vcf files ?
0
gravatar for Picasa
2.9 years ago by
Picasa390
Picasa390 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 • 1.6k views
ADD COMMENTlink modified 3 months ago by jaybee10 • written 2.9 years ago by Picasa390

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 2.9 years ago • written 2.9 years ago by ivivek_ngs4.8k

Hi,

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

ADD REPLYlink written 2.9 years ago by Picasa390

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 2.9 years ago by ivivek_ngs4.8k

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

ADD REPLYlink modified 2.7 years ago • written 2.7 years ago by morovatunc400
0
gravatar for jaybee
3 months ago by
jaybee10
South Korea
jaybee10 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 3 months ago by jaybee10
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: 1508 users visited in the last hour