SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
ChIP-Seq: Enabling Data Analysis on High-Throughput Data in Large Data Depository Usi Newsbot! Literature Watch 1 04-18-2018 09:50 PM
Newbie bwa question (HiSeq data)Hi, Sorry for the newbie question! I am trying to c bsmith030465 Bioinformatics 6 02-26-2015 09:08 AM
question about Solid data analysis using bowtie? wangleibio RNA Sequencing 5 03-26-2013 06:32 AM
question about formating of data mathew Helicos / Direct Genomics 0 10-20-2012 09:54 PM

Reply
 
Thread Tools
Old 11-30-2016, 12:15 PM   #1
lre1234
Senior Member
 
Location: US

Join Date: Aug 2011
Posts: 106
Default Data Analysis question

Hi all,
I'm going through some exome-seq data (10 samples), and following a "standard" pipeline - BWA, HaplotypeCaller, GenotypeGVCF, etc.. nothing too fancy.

My question is on the VQSR step. From my understanding, this works best when you have lots of samples to look at, but since I only have 10 samples, would it be better to simply do some hard-filtering on the dataset? Alternatively, I was thinking of bringing in some additional samples (such as a bunch of 1K genome exomes) and add these into the mix. This would increase the number of samples. Does anyone have any thoughts on this? Also, how many samples would be truly need to do the variant filtering with the VQSR rather than hard-filtering? (I'm having a little bit of a hard time finding what the minimum number of samples for it to work reliably).

Thanks for any advice
lre1234 is offline   Reply With Quote
Old 12-01-2016, 06:32 AM   #2
vivek_
PhD Student
 
Location: Denmark

Join Date: Jul 2012
Posts: 164
Default

Adding 1000 genomes samples might be counterintuitive because you probably already use them as a training set to build the VQSR classifier. There shouldn't be an issue with using hard filtering but regarding sample size for VQSR, you might ask the developers on the GATK forum to get the right answer.
vivek_ is offline   Reply With Quote
Reply

Tags
exome analsys, haplotypecaller, vqsr

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 12:33 PM.


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