|
1
|
1 <tool id="muTect" name="muTect" version="1.0.0">
|
|
|
2 <description>identify somatic point mutations</description>
|
|
|
3 <requirements>
|
|
3
|
4 <requirement type="set_environment">MUTECTPATH</requirement>
|
|
1
|
5 </requirements>
|
|
|
6 <command interpreter="java">
|
|
|
7 -Xmx2g -jar \$MUTECTPATH/muTect-1.1.1.jar
|
|
|
8 --analysis_type MuTect
|
|
|
9 --reference_sequence $reference
|
|
|
10 #if ($cosmic.has_cosmic):
|
|
|
11 --cosmic $cosmic.cosmicfile
|
|
|
12 #end if
|
|
|
13 --dbsnp $dbsnp
|
|
|
14 #if str($intervals) != ""
|
|
|
15 --intervals $intervals
|
|
|
16 #end if
|
|
|
17 #if str($intervals) == "" and str($intervalfile) != ""
|
|
|
18 --intervals $intervalfile
|
|
|
19 #end if
|
|
|
20 --input_file:normal $normal
|
|
|
21 --input_file:tumor $tumor
|
|
|
22 --out call_stats.txt
|
|
|
23 --coverage_file coverage.wig.txt
|
|
|
24 </command>
|
|
|
25
|
|
|
26 <inputs>
|
|
|
27 <param name="reference" type="data" format="bed, fasta" label="Select a reference genome">
|
|
|
28 </param>
|
|
|
29 <param name="dbsnp" type="data" format="vcf" label="dbsnp.vcf file" />
|
|
|
30 <param name="intervals" type="text" label="Intervals to process" help="In format 'chr1:1500-2500; chr2:2500-3500', separated by semicolons" optional="true"/>
|
|
|
31 <param name="intervalfile" type="data" format="txt" label="Intervals to process (.txt file)" help="'chr1:1500-2500', one entry per line" optional="true"/>
|
|
|
32 <param name="normal" type="bam" label="Normal BAM file"/>
|
|
|
33 <param name="tumor" type="bam" label="Tumor BAM file"/>
|
|
|
34 <conditional name="cosmic">
|
|
|
35 <param name="has_cosmic" type="boolean" label="Cosmic vcf file available for specimen" value="false" />
|
|
|
36 <when value="true">
|
|
|
37 <param name="cosmicfile" type="data" format="vcf" label="Cosmic vcf file"/>
|
|
|
38 </when>
|
|
|
39 </conditional>
|
|
|
40
|
|
|
41
|
|
|
42 </inputs>
|
|
|
43
|
|
|
44 <outputs>
|
|
|
45 <data format="txt" name="callstats" label="Detailed caller output" from_work_dir="call_stats.txt"/>
|
|
|
46 <data format="txt" name="coverage" label="Filename for coverage output" from_work_dir="coverage.wig.txt"/>
|
|
|
47 </data>
|
|
|
48 </outputs>
|
|
3
|
49
|
|
|
50
|
|
1
|
51 <help>
|
|
|
52
|
|
|
53 .. class:: infomark
|
|
|
54
|
|
|
55 **License**
|
|
|
56
|
|
|
57 TmuTect Revision: 44830
|
|
|
58 GATK Revision: 2.1-202-g2fe6a31
|
|
|
59 Copyright (c) 2012 The Broad Institute
|
|
|
60 Please view our documentation at http://www.broadinstitute.org/gsa/wiki
|
|
|
61 For support, please view our support site at http://getsatisfaction.com/gsa
|
|
|
62
|
|
|
63 **What it does**
|
|
|
64
|
|
|
65 MuTect is a method developed at the Broad Institute for the reliable
|
|
|
66 and accurate identification of somatic point mutations in next generation
|
|
|
67 sequencing data of cancer genomes.
|
|
|
68
|
|
|
69
|
|
|
70
|
|
|
71 **Input**
|
|
|
72
|
|
|
73 Reference Genome:
|
|
|
74
|
|
|
75 Fasta file of ref gnome
|
|
|
76
|
|
|
77 Normal Sample:
|
|
|
78
|
|
|
79 bam normal sample
|
|
|
80
|
|
|
81 Tumor Sample:
|
|
|
82
|
|
|
83 bam tumor sample
|
|
|
84
|
|
|
85
|
|
|
86 Intervals:
|
|
|
87
|
|
|
88 A list of genomic intervals over which to operate.
|
|
|
89
|
|
|
90
|
|
|
91
|
|
|
92 **Output**
|
|
|
93
|
|
|
94 Caller Output:
|
|
|
95
|
|
|
96 Detailed caller output
|
|
|
97
|
|
|
98 Coverage Output:
|
|
|
99
|
|
|
100 write out coverage in WIGGLE format to this file
|
|
|
101
|
|
|
102
|
|
|
103 </help>
|
|
|
104 </tool>
|
|
|
105
|
|
|
106
|
|
|
107
|
|
|
108
|
