Question: Pindel Somatic Calls?
1
gravatar for richard.deborja
5.4 years ago by
Canada
richard.deborja110 wrote:

Does Pindel have a somatic mode for indel calling? I figure I would ask the community before I start making my own post-processing script.

pindel • 4.0k views
ADD COMMENTlink modified 5.2 years ago by Danielk570 • written 5.4 years ago by richard.deborja110
2

My custom GATK walker (with build instructions) is now available at github: https://github.com/dakl/gatk

Let me know if you have any issues with it. 

ADD REPLYlink written 5.2 years ago by Danielk570
"https://github.com/dakl/gatk" was compiled using commond :mvn verify -P\!queue, and external-example-1.0-SNAPSHOT.jar was generated.

But,  "java -jar ./public/external-example/target/external-example-1.0-SNAPSHOT.jar -T SomaticPindelFilter " reported errors:

##### ERROR ------------------------------------------------------------------------------------------
##### ERROR A USER ERROR has occurred (version exported): 
##### ERROR
##### ERROR This means that one or more arguments or inputs in your command are incorrect.
##### ERROR The error message below tells you what is the problem.
##### ERROR
##### ERROR If the problem is an invalid argument, please check the online documentation guide
##### ERROR (or rerun your command with --help) to view allowable command-line arguments for this tool.
##### ERROR
##### ERROR Visit our website and forum for extensive documentation and answers to 
##### ERROR commonly asked questions http://www.broadinstitute.org/gatk
##### ERROR
##### ERROR Please do NOT post this error to the GATK forum unless you have really tried to fix it yourself.
##### ERROR
##### ERROR MESSAGE: Invalid command line: Malformed walker argument: Could not find walker with name: SomaticPindelFilter
##### ERROR ------------------------------------------------------------------------------------------

ADD REPLYlink modified 4.7 years ago • written 4.7 years ago by xingyongma0

How did you compile? 

I get this, which is the expected output: 

$ git clone https://github.com/dakl/gatk.git

$ cd gatk

$ mvn verify -P\!queue

$ java -jar ./public/external-example/target/external-example-1.0-SNAPSHOT.jar -T SomaticPindelFilter

##### ERROR ------------------------------------------------------------------------------------------

##### ERROR A USER ERROR has occurred (version 3.2-36-gfe24d6a): 

##### ERROR

##### ERROR This means that one or more arguments or inputs in your command are incorrect.

##### ERROR The error message below tells you what is the problem.

##### ERROR

##### ERROR If the problem is an invalid argument, please check the online documentation guide

##### ERROR (or rerun your command with --help) to view allowable command-line arguments for this tool.

##### ERROR

##### ERROR Visit our website and forum for extensive documentation and answers to 

##### ERROR commonly asked questions http://www.broadinstitute.org/gatk

##### ERROR

##### ERROR Please do NOT post this error to the GATK forum unless you have really tried to fix it yourself.

##### ERROR

##### ERROR MESSAGE: Argument with name '--normalid' (-NID) is missing.

##### ERROR Argument with name '--tumorid' (-TID) is missing.

##### ERROR Argument with name '--variant' (-V) is missing.

##### ERROR ------------------------------------------------------------------------------------------
ADD REPLYlink written 4.7 years ago by Danielk570

I download the ZIP file , and compile as following as you:

$ cd gatk

$ mvn verify -P\!queue

$ java -jar ./public/external-example/target/external-example-1.0-SNAPSHOT.jar -T SomaticPindelFilter

 

[INFO] GATK Root .......................................... SUCCESS [  1.023 s]
[INFO] GATK Aggregator .................................... SUCCESS [  0.983 s]
[INFO] GATK GSALib ........................................ SUCCESS [  0.696 s]
[INFO] GATK Utils ......................................... SUCCESS [  1.364 s]
[INFO] GATK Engine ........................................ SUCCESS [  0.465 s]
[INFO] GATK Tools Public .................................. SUCCESS [ 12.864 s]
[INFO] External Example ................................... SUCCESS [  9.591 s]
[INFO] GATK Aggregator Public ............................. SUCCESS [  0.018 s]
[INFO] GATK Tools Protected ............................... SUCCESS [ 10.741 s]
[INFO] GATK Package Distribution .......................... SUCCESS [01:25 min]
[INFO] GATK Aggregator Protected .......................... SUCCESS [  0.014 s]

 

 

 

ADD REPLYlink written 4.7 years ago by xingyongma0

It worked:

git clone https://github.com/dakl/gatk.git


​instead

 

It is very useful

 

ADD REPLYlink modified 4.7 years ago • written 4.7 years ago by xingyongma0

Glad it worked - I will look into why the zip is not working. 

 

ADD REPLYlink written 4.7 years ago by Danielk570

I get the same (functional) results when downloading the master.zip, so I won't look into this any further. 

ADD REPLYlink written 4.7 years ago by Danielk570

I wonder if the SomaticPindelFilter tool is currently available in the repo at https://github.com/dakl/gatk.

When I enter the command: java -jar GenomeAnalysisTK-Klevebring.jar -T SomaticPindelFilter

I get an error message that contains the following line: ##### ERROR MESSAGE: Invalid command line: Malformed walker argument: Could not find walker with name: SomaticPindelFilter

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by limax0

I am also interested. Could you please release SomaticPindelFilter to the world?

ADD REPLYlink written 5.2 years ago by tayebwajb90
4
gravatar for Danielk
5.4 years ago by
Danielk570
Karolinska Institutet, Stockholm, Sweden
Danielk570 wrote:

Edit June 8th, 2019 ( by Kevin Blighe ): this answer (below) is old and now outdated. Latest versions of Pindel can perform somatic indel calling

-------------------------------------------

There's no way to do that in pindel, downstream filtering is required.

I use a couple of filters:

  1. Reasonable coverage in T and N (btwn 25-1000) and frequency in normal < 15% (settable on cli).
  2. Collect all variant that pass (1).
  3. Perform Fisher's exact test comparing ref and alt counts in T and N.
  4. Adjust for multiple testing using BH/FDR method.
  5. Keep variants with adjusted p < 0.05 (settable on cli).

I've implemented this in a custom GATK walker called SomaticPindelFilter. If you're interested, I could release this to the world.

ADD COMMENTlink modified 10 weeks ago by Kevin Blighe46k • written 5.4 years ago by Danielk570

I'm interested in your filter. How can I get it?

ADD REPLYlink written 5.3 years ago by volkansevim10
1
gravatar for mdm-two
5.4 years ago by
mdm-two230
mdm-two230 wrote:

The setup is described in the online user manual:

http://gmt.genome.wustl.edu/pindel/current/user-manual.html

There is a also a screencast from Kai Ye:

https://www.youtube.com/watch?v=SQKHXfPSJfE&feature=player_embedded

ADD COMMENTlink written 5.4 years ago by mdm-two230
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: 1657 users visited in the last hour