SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Merging SNPs and indels VCF from strelka for ANNOVAR annotation lbeltrame Bioinformatics 0 12-11-2012 12:37 AM
Merging vcf files kjaja Bioinformatics 0 01-12-2012 11:06 AM
merging VCF files kjaja Bioinformatics 2 12-20-2011 06:50 AM
merging 2 sam files papori Bioinformatics 0 07-29-2011 04:44 AM
failed when merging the example VCF files using VCFtools jianfeng.mao Bioinformatics 1 02-02-2011 05:44 PM

Reply
 
Thread Tools
Old 06-04-2013, 03:01 AM   #1
marco12345
Member
 
Location: barcelona

Join Date: Mar 2013
Posts: 18
Question vcf files merging

Hi all,

I was wondering which is the easiest way to merge multiple vcf files (outputs of Varscan). I have some vcf files listing all the called snps for an individual (1 snp per line). I would like to put together all this info in a n output that shows for every line a snp (original base, variant base/bases, and eventually number of individuals sharing the variant).
I tried with GATK CombineVariants and had the ouput I wanted, but in case of different variant at the same position, only one of the variant bases was reported.

I am working on a Mac's Terminal, OS 10.6.

If needed I have access also to a Linux OS.

Thanks in advance!
marco12345 is offline   Reply With Quote
Old 06-04-2013, 03:08 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,795
Default

http://vcftools.sourceforge.net/

http://vcftools.sourceforge.net/docs.html#merge
GenoMax is offline   Reply With Quote
Old 06-04-2013, 06:33 AM   #3
marco12345
Member
 
Location: barcelona

Join Date: Mar 2013
Posts: 18
Default

Sorry GenoMax, I forgot to mention that I tried vcftools/vcf-merge, but still the output does not show different variants for the same position.
marco12345 is offline   Reply With Quote
Old 06-13-2013, 05:30 AM   #4
ernfrid
Member
 
Location: Saint Louis

Join Date: Jul 2011
Posts: 26
Default

You could try joinx's vcf-merge command which I'm fairly confident reports multiple variant alleles at a position.

The code is available as a debian package from:http://gmt.genome.wustl.edu/joinx/1.7/ but you may want to compile from source code (the github page is linked off the gmt page).
ernfrid is offline   Reply With Quote
Reply

Tags
osx, varscan, vcf, vcftools

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:45 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO