|
2
|
1 <!--Sequenza - developed by Jocelyn Brayet <jocelyn.brayet@curie.fr>
|
|
|
2 Copyright (C) 2015 Institut Curie
|
|
|
3
|
|
|
4 This program is free software: you can redistribute it and/or modify
|
|
|
5 it under the terms of the GNU General Public License as published by
|
|
|
6 the Free Software Foundation, either version 3 of the License, or
|
|
|
7 (at your option) any later version.
|
|
|
8
|
|
|
9 This program is distributed in the hope that it will be useful,
|
|
|
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
|
12 GNU General Public License for more details.
|
|
|
13
|
|
|
14 You should have received a copy of the GNU General Public License
|
|
|
15 along with this program. If not, see <http://www.gnu.org/licenses/>.-->
|
|
|
16 <tool id="sequenza_tool" name="Sequenza" version="1.2">
|
|
|
17 <description>allele-specific copy number and mutation profiles</description>
|
|
|
18 <command interpreter="python">
|
|
|
19 sequenza_wrapper.py -name $sampleName -outGalaxy $HTMLFile -format $fileFormat.format -estimation $usePersonalEstimation.estimation -gcContent $gc_content_file
|
|
|
20 #if $fileFormat.format == "BAM":
|
|
|
21 #if $fileFormat.reference_source.reference_source_selector=="cached":
|
|
|
22 -selector $fileFormat.reference_source.reference_source_selector
|
|
|
23 -normal $fileFormat.reference_source.normal_file
|
|
|
24 -tumor $fileFormat.reference_source.tumor_file
|
|
|
25 -ref_file $fileFormat.reference_source.ref_file.fields.path
|
|
|
26 -samtools_options '
|
|
|
27 $fileFormat.skip_anomalous_read_pairs
|
|
|
28 $fileFormat.disable_probabilistic_realignment
|
|
|
29 -d "$fileFormat.max_reads_per_bam"
|
|
|
30 -q "$fileFormat.minimum_mapping_quality"
|
|
|
31 -Q "$fileFormat.minimum_base_quality"
|
|
|
32 '
|
|
|
33 #else:
|
|
|
34 -selector $fileFormat.reference_source.reference_source_selector
|
|
|
35 -normal $fileFormat.reference_source.normal_file
|
|
|
36 -tumor $fileFormat.reference_source.tumor_file
|
|
|
37 -ref_file $fileFormat.reference_source.ref_file
|
|
|
38 #end if
|
|
|
39 #else:
|
|
|
40 -normal $fileFormat.normal_file
|
|
|
41 -tumor $fileFormat.tumor_file
|
|
|
42 #end if
|
|
|
43 #if $usePersonalEstimation.estimation == "yes":
|
|
|
44 -cellularity $usePersonalEstimation.cellularity
|
|
|
45 -ploidy $usePersonalEstimation.ploidy
|
|
|
46 #end if
|
|
|
47 </command>
|
|
|
48 <inputs>
|
|
|
49 <param name="sampleName" type="text" value="sample" size="30" label="Sample name">
|
|
|
50 <sanitizer invalid_char="">
|
|
|
51 <valid initial="string.letters,string.digits"><add value="_"/></valid>
|
|
|
52 </sanitizer>
|
|
|
53 </param>
|
|
|
54 <conditional name="fileFormat">
|
|
|
55 <param name="format" type="select" label="File format" >
|
|
|
56 <option value="BAM" selected="true">BAM</option>
|
|
|
57 <option value="pileup" >Pileup</option>
|
|
|
58 <option value="pileup_gz" >Pileup.gz</option>
|
|
|
59 </param>
|
|
|
60 <when value="BAM">
|
|
|
61 <conditional name="reference_source">
|
|
|
62 <param name="reference_source_selector" type="select" label="Choose the source for the reference list">
|
|
|
63 <option value="cached">Locally cached</option>
|
|
|
64 <option value="history">History</option>
|
|
|
65 </param>
|
|
|
66 <when value="cached">
|
|
|
67 <param name="normal_file" type="data" format="bam" label="Normal BAM file" >
|
|
|
68 <validator type="unspecified_build" />
|
|
|
69 <validator type="dataset_metadata_in_data_table" table_name="sam_fa_indexes" metadata_name="dbkey" metadata_column="value" message="Sequences are not currently available for the specified build." />
|
|
|
70 </param>
|
|
|
71 <param name="tumor_file" type="data" format="bam" label="Tumor BAM file" >
|
|
|
72 <validator type="unspecified_build" />
|
|
|
73 <validator type="dataset_metadata_in_data_table" table_name="sam_fa_indexes" metadata_name="dbkey" metadata_column="value" message="Sequences are not currently available for the specified build." />
|
|
|
74 </param>
|
|
|
75 <param name="ref_file" type="select" label="Using reference genome">
|
|
|
76 <options from_data_table="sam_fa_indexes">
|
|
|
77 </options>
|
|
|
78 </param>
|
|
|
79 </when>
|
|
|
80 <when value="history">
|
|
|
81 <param name="normal_file" type="data" format="bam" label="Normal BAM file" >
|
|
|
82 <validator type="metadata" check="bam_index" message="Metadata missing, click the pencil icon in the history item and use the auto-detect feature to correct this issue."/>
|
|
|
83 </param>
|
|
|
84 <param name="tumor_file" type="data" format="bam" label="Tumor BAM file" >
|
|
|
85 <validator type="metadata" check="bam_index" message="Metadata missing, click the pencil icon in the history item and use the auto-detect feature to correct this issue."/>
|
|
|
86 </param>
|
|
|
87 <param name="ref_file" type="data" format="fasta" label="Using reference file" />
|
|
|
88 </when>
|
|
|
89 </conditional>
|
|
|
90 <param name="skip_anomalous_read_pairs" type="boolean" truevalue="-A" falsevalue="" checked="False" label="Do not skip anomalous read pairs in variant calling" />
|
|
|
91 <param name="disable_probabilistic_realignment" type="boolean" truevalue="-B" falsevalue="" checked="False" label="Disable probabilistic realignment for the computation of base alignment quality (BAQ)" />
|
|
|
92 <param name="max_reads_per_bam" type="integer" value="100000" label="Max reads per BAM" />
|
|
|
93 <param name="minimum_mapping_quality" type="integer" value="20" label="Minimum mapping quality for an alignment to be used" />
|
|
|
94 <param name="minimum_base_quality" type="integer" value="13" label="Minimum base quality for a base to be considered" />
|
|
|
95 </when>
|
|
|
96 <when value="pileup">
|
|
|
97 <param name="normal_file" type="data" format="pileup" label="Normal pileup file" />
|
|
|
98 <param name="tumor_file" type="data" format="pileup" label="Tumor pileup file" />
|
|
|
99 </when>
|
|
|
100 <when value="pileup_gz">
|
|
|
101 <param name="normal_file" type="data" format="pileup.gz" label="Normal pileup.gz file" />
|
|
|
102 <param name="tumor_file" type="data" format="pileup.gz" label="Tumor pileup.gz file" />
|
|
|
103 </when>
|
|
|
104 </conditional >
|
|
|
105 <param name="gc_content_file" type="data" format="txt" label="GC content file" />
|
|
|
106 <conditional name="usePersonalEstimation">
|
|
|
107 <param name="estimation" type="select" label="Do you want to use personal cellularity and ploidy? Otherwise Sequenza estimates cellularity and ploidy" >
|
|
|
108 <option value="no" selected="true">No</option>
|
|
|
109 <option value="yes" >Yes</option>
|
|
|
110 </param>
|
|
|
111 <when value="no" />
|
|
|
112 <when value="yes">
|
|
|
113 <param name="cellularity" type="integer" value="30" min="0" max="100" label="Cellularity used (%)" />
|
|
|
114 <param name="ploidy" type="integer" value="2" label="Ploidy used" />
|
|
|
115 </when>
|
|
|
116 </conditional >
|
|
|
117 </inputs>
|
|
|
118 <outputs>
|
|
|
119 <data format="html" name="HTMLFile" label="HTML output - Sequenza results" />
|
|
|
120 </outputs>
|
|
|
121 <tests>
|
|
|
122 <test>
|
|
|
123 <param name="sampleName" value="test_curie"/>
|
|
|
124 <param name="fileFormat" value="pileup"/>
|
|
|
125 <param name="normal_file" value="constit_2.pileup"/>
|
|
|
126 <param name="tumor_file" value="tumor_2.pileup"/>
|
|
|
127 <param name="gc_content_file" value="hg19_chro1.gc50Base.txt"/>
|
|
|
128 <param name="estimation" value="yes"/>
|
|
|
129 <param name="cellularity" value="30"/>
|
|
|
130 <param name="ploidy" value="2"/>
|
|
|
131 <output name="HTMLFile" file="test_sequenza_1.dat" ftype="html"/>
|
|
|
132 </test>
|
|
|
133 </tests>
|
|
|
134 <help>
|
|
|
135
|
|
|
136 **What it does**
|
|
|
137
|
|
|
138 Tools to analyze genomic sequencing data from paired normal-tumor samples, including cellularity and ploidy estimation; mutation and copy number (allele-specific and total copy number) detection, quantification and visualization.
|
|
|
139
|
|
|
140 </help>
|
|
|
141 <citations>
|
|
|
142 <citation type="bibtex">@article{Favero01012015,
|
|
|
143 author = {Favero, F. and Joshi, T. and Marquard, A. M. and Birkbak, N. J. and Krzystanek, M. and Li, Q. and Szallasi, Z. and Eklund, A. C.},
|
|
|
144 title = {Sequenza: allele-specific copy number and mutation profiles from tumor sequencing data},
|
|
|
145 volume = {26},
|
|
|
146 number = {1},
|
|
|
147 pages = {64-70},
|
|
|
148 year = {2015},
|
|
|
149 doi = {10.1093/annonc/mdu479},
|
|
|
150 URL = {http://annonc.oxfordjournals.org/content/26/1/64.abstract},
|
|
|
151 eprint = {http://annonc.oxfordjournals.org/content/26/1/64.full.pdf+html},
|
|
|
152 journal = {Annals of Oncology}
|
|
|
153 }</citation>
|
|
|
154 </citations>
|
|
|
155 </tool>
|
