SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
varscan2: copynumber detection mrfox Bioinformatics 2 07-10-2013 12:07 PM
too many indels output by varscan2 mrfox Bioinformatics 2 09-25-2012 08:17 AM
How to plot datasets from VarScan2 in R. pablosantamarta Bioinformatics 1 03-22-2012 03:47 AM
Good variant filters christophpale Bioinformatics 0 11-21-2011 02:20 PM
SV emPCR Filters Anthony.287 454 Pyrosequencing 0 11-07-2011 06:41 AM

Reply
 
Thread Tools
Old 10-04-2012, 10:45 AM   #1
mrfox
Senior Member
 
Location: USA

Join Date: Aug 2010
Posts: 103
Default varscan2: filters

Hi all and hopefully Dan,

I am still not very sure about how to use the Varscan2 filters for SNVs and Indels obtained. Below are my questions, any advice on the clarification of using the filters will be appreciated.

1) For a single sample,

Step 1. Use pileup2snp and pileup2indel call SNVs and Indels.

Step2. Run varscan2 filter on SNVs using the Indels as an input. Then filter indels using java -jar varscan2 filter output.indel >filtered.indel.

Aternatively,
Step2. Another option is to filter indel first and use the filtered indels as input to filter SNVs.

Which one is better for step 2?

2) For a Normal/Tumor sample,
From the context of the Varscan2 instruction, the "somaticFilter" is designed only for filtering SNVs, not for indels. So there is actually no filter for Indels. The Indel file obtained from calling Varscan2 somatic will be the final reported Indel. Is that correct? Is somaticFilter applicable to Indels?

Of course, we can change any somatic indel with read2>0 into Germline.

Thank you all!
mrfox is offline   Reply With Quote
Old 10-09-2012, 08:40 AM   #2
dkoboldt
Member
 
Location: St. Louis

Join Date: Mar 2009
Posts: 62
Default

Hao,

Thank you for the question. First of all, I would use *unfiltered* indels to filter SNPs (as long as your variant calling criteria are reasonable), because your goal is to remove variants that are due to alignment issues, and even unfiltered indel calls are a guide to these. But use your judgment: if the unfiltered indels remove a significant portion of dbSNPs, then there's a problem.

The somaticFilter works on somatic indel as well as SNV files from VarScan somatic. I do recommend that you use it, but you're correct to think of applying your own filters to those (especially to catch germline variants). Any time we see "somatic" indels with evidence in the normal sample, or that match dbSNP variants, we remove them.
dkoboldt 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 11:22 AM.


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