Question: about the use of varscan, somaticsnipper and samtools
0
gravatar for Bogdan
3.5 years ago by
Bogdan850
Palo Alto, CA, USA
Bogdan850 wrote:

Dear all,

please could you advise : before the use of eitther varscan, somaticsnipper or samtools to call the somatic mutations, would we need to :

a, mark the duplicates b. realign the data around indels (as we would do typically with GATK, before running Mutect) ?

many thanks !

bogdan

sequencing snp • 1.1k views
ADD COMMENTlink modified 3.5 years ago • written 3.5 years ago by Bogdan850
1
gravatar for Chris Miller
3.5 years ago by
Chris Miller21k
Washington University in St. Louis, MO
Chris Miller21k wrote:

a) yes, mark duplicates in your bam files

b) no, local indel realignment is not necessary

c) I wouldn't rely solely on somatic sniper or varscan for detecting indels. They'll provide some useful information, but Pindel-C and GATK generally perform better on indel calling.

ADD COMMENTlink written 3.5 years ago by Chris Miller21k

Hi Chris, Can you elaborate on why you say that local indel realignment is not necessary? Is it that you see no benefit, or that some callers do not need the local realignment preprocessing in order to deal with non-canonical alignments? Our (RTG) somatic caller internally does something similar to local realignment to candidate haplotypes when calling complex regions such as indels or dense snps, so realignment is not needed as a separate preprocessing step, and I imagine that callers incorporating local assembly also do not need realignment preprocessing.

ADD REPLYlink written 3.5 years ago by Len Trigg1.3k

Definitely the latter. I can't comment on whether there is any benefit, but the specific tools he mentioned do not require it.

ADD REPLYlink written 3.5 years ago by Chris Miller21k
0
gravatar for Bogdan
3.5 years ago by
Bogdan850
Palo Alto, CA, USA
Bogdan850 wrote:

Thanks Chris ! Talking about somatic variant calling, which one of these 4 would you recommend : STRELKA, MUTECT, VARSCAN or SOMATICSNIPPER ? Or perhaps some other somatic variant caller ?

what about the SV detection or CNV ? Any recommendation is welcome ! many thanks !

ADD COMMENTlink written 3.5 years ago by Bogdan850
  • We use all four of those and combine them in various ways to up our sensitivity.
  • For SVs, I like Manta or Lumpy. For CNV, I'm partial to my own package - copyCat (https://github.com/chrisamiller/copycat)
ADD REPLYlink written 3.5 years ago by Chris Miller21k
0
gravatar for Bogdan
3.5 years ago by
Bogdan850
Palo Alto, CA, USA
Bogdan850 wrote:

Chris, thank you very much !

ADD COMMENTlink written 3.5 years ago by Bogdan850
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1096 users visited in the last hour