SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Viewing multiple VCF files in Artemis coldturkey Bioinformatics 0 02-01-2012 01:45 AM
Merging multiple BAM files unibegenomics Bioinformatics 1 08-25-2011 02:03 AM
GATK UnifiedGenotyper calling way too many SNPs in vcf swbarnes2 Bioinformatics 0 08-17-2011 01:33 PM
GATK calling Merged Bam files jayce_ocean Bioinformatics 3 03-16-2011 12:15 AM
calling SNPs asankaf General 2 02-04-2009 06:45 PM

Reply
 
Thread Tools
Old 04-18-2011, 08:49 AM   #1
newbietonextgen
Member
 
Location: USA

Join Date: Nov 2010
Posts: 56
Default Calling multiple BAM files for SNPs and vcf

Hi all,

I am in the process of calling SNPs using GATK. I do know that the unifiedgenotyper will call from multiple BAM files. The question is about the resulting VCF file.
1) Does the VCF file separate out the SNps? In other words, will i know which BAM sample file contributed to what SNP in the VCF file? In such as way i could just parse back individual sample SNPs back, if needed.
2) Also it it better to call SNPs from individual files or pool all the BAM files from the same treatment to get better depth?

The BAMs have unique Read Groups. Any help will a appreciated.

Thanks

Last edited by newbietonextgen; 04-18-2011 at 09:51 AM.
newbietonextgen is offline   Reply With Quote
Old 04-18-2011, 06:56 PM   #2
jstjohn
Member
 
Location: San Francisco, CA

Join Date: Jun 2010
Posts: 35
Default

1. I think GATK uses the Read Group tag to partition the data into separate organisms for variant calling. If you want them treated as the same sample then you give them all the same read group tag, while if you want them treated differently then give them different read group tags. I am actually not 100% sure that "read group" is the right tag for this, so someone please correct me if thats not the case.

There is a program in the Picard package to manipulate the read group along with other tags in bam files.

2. Read through this page which describes the different suggested protocals in order of how good they are thought to perform vs computational time involved. I think the best (and slowest) option is to pool everything. The page gives a good overview of the suggested set of steps to follow to call variants given different computational commitments:
http://www.broadinstitute.org/gsa/wi...th_the_GATK_v2
jstjohn is offline   Reply With Quote
Old 04-19-2011, 06:41 AM   #3
newbietonextgen
Member
 
Location: USA

Join Date: Nov 2010
Posts: 56
Default

Yes, read groups may help in the calling process but when you look a vcf file, can you actually distinguish which BAM file gave rise to particular SNPs?
newbietonextgen is offline   Reply With Quote
Old 04-19-2011, 11:29 AM   #4
d17
Member
 
Location: United States

Join Date: Sep 2008
Posts: 27
Default

Quote:
Originally Posted by newbietonextgen View Post
can you actually distinguish which BAM file gave rise to particular SNPs?
If you are doing multi-sample calling with each sample in a separate BAM, it doesn't always make sense to say that one BAM file gave rise to a particular SNP. For example you could have three samples each of which has weak evidence for a particular SNP: not enough to call the SNP if you are analyzing the samples separately but when you pool them together the combined evidence is sufficient to call a SNP.
d17 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 08:26 AM.


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