Biostar Beta. Not for public use.
No Somatic SNV Calling with GATK?
1
Entering edit mode
2.2 years ago
sichan • 70
Canada

Hi,

I see that GATK should not be used to call somatic variants from tumor/normal pairs:
http://gatkforums.broadinstitute.org/discussion/2405/can-i-use-unifiedgenotyper-to-call-somatic-variants-in-tumor-control-tissues

Is this still true today? I ask because someone who is no longer in our group generated 'somatic calls' by

  • calling SNPs on the normal
  • calling SNPs on the matched tumor
  • subtracting the normal SNPs from the tumor SNPs

I'm trying to determine how useful those calls are (probably not much).

This post from ~2 years ago also suggests that GATK should not be used for somatic calls:
https://www.biostars.org/p/80826/

Thank you.

SNP • 2.0k views
ADD COMMENTlink
0
Entering edit mode

Thanks everyone! Yes, I forgot to mention that I had already started running MuTect, which is also from Broad.

ADD REPLYlink
0
Entering edit mode

Try to refrain from posting thanks as an answer. An upvote or accepted answer will do just fine!

ADD REPLYlink
2
Entering edit mode
17 months ago
National Institutes of Health, Bethesda…

You should definitely abandon those calls and use a tool designed for somatic variant calling. There are numerous tools to choose from.

ADD COMMENTlink
1
Entering edit mode

You can use MuTect, which is based on GATK and is available for download from the GATK website.

ADD REPLYlink
1
Entering edit mode

Yes, MuTect is one of many and a popular one. On the downside, it is pretty slow unless one takes the time to run it in parallel.

ADD REPLYlink
1
Entering edit mode

I could't agree more. This naive subtractive approach is in no way comparable to say VarScan2.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1