SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Using split/merge .bams for Cufflinks DerSeb Bioinformatics 2 03-15-2012 11:42 AM
Somatic SNP wanguan2000 Bioinformatics 2 03-13-2012 08:59 PM
Cufflinks reporting differnt FPKMS for the same gene nwalwort Bioinformatics 6 02-09-2012 08:47 AM
Problems with SPP read.bams.tags EdinG Bioinformatics 0 10-12-2011 03:27 AM
TCGA dbgap: where's the bams? Richard Finney Bioinformatics 0 06-13-2011 09:18 AM

Reply
 
Thread Tools
Old 02-06-2012, 01:31 PM   #1
Mamun
Junior Member
 
Location: UK

Join Date: Jul 2011
Posts: 2
Default Somatic calling comparing differnt platform BAMs

Hi guys,

I am trying to make some somatic calling using ( varscan / mpileup ) a tumor
sample from whole exome sequencing data and Normal sample from whole
genome data.

Size of the BAM files

Whole exome tumor ( 6.7 gb )
Whole genome tumor ( 18 gb )

I have ~6 million calls from this procedure.

I was wondering is it appropriate to perform somatic calling with two samples from
different platform.

thanks in advance,
Mamun
Mamun is offline   Reply With Quote
Old 02-06-2012, 04:23 PM   #2
adaptivegenome
Super Moderator
 
Location: US

Join Date: Nov 2009
Posts: 437
Default

If this is the only data you have, then you can certainly try. However, it is probably a bad idea for many reasons. Tools are often optimized to call genotypes and with the big differences in coverage (exome is usually ~100X; genome is usually ~30X) this could skew results especially when you are assessing changes that arise in the tumor.
adaptivegenome is offline   Reply With Quote
Old 02-07-2012, 07:24 AM   #3
sdvie
Member
 
Location: Spain

Join Date: Jul 2010
Posts: 68
Default

just one thought:

The exome sample must have been submitted to some enrichment in the experimental part of the workflow. Usually you would have a .bed file containing the enriched regions. So one thing to do, would be to select from the whole genome data only the regions that were enriched for the exome. Then, the coverage might still be different between the two bam files.... you might consider adjusting the amount of reads from the exome experiment by downsampling the number of reads in the exome bam file, as discussed here.

Last edited by sdvie; 02-07-2012 at 07:42 AM.
sdvie 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 04:35 PM.


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