No Somatic SNV Calling with GATK?
1
1
Entering edit mode
8.9 years ago
sichan ▴ 90

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: Difference Between Somatic And Germline Variant Calling?

Thank you

SNP • 3.2k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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

ADD REPLY
2
Entering edit mode
8.9 years ago

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

ADD COMMENT
1
Entering edit mode

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

ADD REPLY
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 REPLY
1
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1421 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6