Question: about the use of varscan, somaticsnipper and samtools
0
gravatar for Bogdan
3.1 years ago by
Bogdan720
Palo Alto, CA, USA
Bogdan720 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 • 968 views
ADD COMMENTlink modified 3.1 years ago • written 3.1 years ago by Bogdan720
1
gravatar for Chris Miller
3.1 years ago by
Chris Miller20k
Washington University in St. Louis, MO
Chris Miller20k 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.1 years ago by Chris Miller20k

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.0 years ago by Len Trigg1.2k

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.0 years ago by Chris Miller20k
0
gravatar for Bogdan
3.1 years ago by
Bogdan720
Palo Alto, CA, USA
Bogdan720 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.1 years ago by Bogdan720
  • 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.1 years ago by Chris Miller20k
0
gravatar for Bogdan
3.1 years ago by
Bogdan720
Palo Alto, CA, USA
Bogdan720 wrote:

Chris, thank you very much !

ADD COMMENTlink written 3.1 years ago by Bogdan720
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: 2254 users visited in the last hour