SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
merging fastq files shilo Illumina/Solexa 8 07-06-2016 02:15 PM
Merging bam files memento Bioinformatics 1 02-17-2012 04:39 PM
merging VCF files kjaja Bioinformatics 2 12-20-2011 07:50 AM
merging 2 sam files papori Bioinformatics 0 07-29-2011 05:44 AM
failed when merging the example VCF files using VCFtools jianfeng.mao Bioinformatics 1 02-02-2011 06:44 PM

Reply
 
Thread Tools
Old 01-12-2012, 12:06 PM   #1
kjaja
Member
 
Location: NY

Join Date: Aug 2011
Posts: 55
Default Merging vcf files

I am noticing a problem when merging two VCF file if the SNP is multiallelic, below is an example of how a multiallelic SNP should look like in a VCF file for the merged samples. You can see that the AF= Allele frequency below takes two values (0.333, 0.667) which corresponds to the frequency of Alternative allele G and the T .

1110696 rs6040355 A G,T 67 PASS NS=2;DP=10;AF=0.333,0.667;AA=T;DB GT:GQP:HQ 1|2:21:6:23,27 2|1:2:0:18,2 2/2:35:4

The problem that I am having is when merging two samples, I do not get AF values for each alternative allele. I tried to do this using VCF tools, below is my command.

perl vcf-merge file1.vcf.gz file2.vcf.gz >merged_files.vcf

Does anyone know how to solve this problem or if there are other tools to merge .vcf files that takes care of this problem?

Thanks,
kjaja is offline   Reply With Quote
Reply

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 10:32 PM.


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