annotate cuffdiff_wrapper.xml @ 1:ebb9a992508d

Remove duplicate cds_diff_output option
author Jim Johnson <jj@umn.edu>
date Mon, 04 Feb 2013 21:23:20 -0600
parents da7241f92ecf
children fdf01b3c1841
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
1 <tool id="cuffdiff_cummerbund" name="Cuffdiff" version="0.0.6">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
2 <!-- Wrapper supports Cuffdiff versions v1.3.0-v2.0 -->
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
3 <description>find significant changes in transcript expression, splicing, and promoter use</description>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
4 <requirements>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
5 <requirement type="package">cufflinks</requirement>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
6 </requirements>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
7 <command interpreter="python">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
8 #set sel_outputs = $output_sel.__str__.split(',')
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
9 cuffdiff_wrapper.py
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
10 --FDR=$fdr
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
11 --num-threads="4"
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
12 --min-alignment-count=$min_alignment_count
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
13
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
14 #if 'cuffdata' in $sel_outputs or not $output_sel:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
15 --cuffdatadir=$cuffdata.extra_files_path
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
16 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
17 #if 'cummeRbund_db' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
18 --cummeRbund_db=$cummeRbund_db
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
19 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
20
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
21 #if 'isoforms_fpkm_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
22 --isoforms_fpkm_tracking_output=$isoforms_fpkm_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
23 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
24 #if 'genes_fpkm_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
25 --genes_fpkm_tracking_output=$genes_fpkm_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
26 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
27 #if 'cds_fpkm_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
28 --cds_fpkm_tracking_output=$cds_fpkm_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
29 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
30 #if 'tss_groups_fpkm_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
31 --tss_groups_fpkm_tracking_output=$tss_groups_fpkm_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
32 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
33 #if 'isoforms_exp_diff' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
34 --isoforms_exp_output=$isoforms_exp_diff
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
35 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
36 #if 'genes_exp_diff' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
37 --genes_exp_output=$genes_exp_diff
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
38 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
39 #if 'tss_groups_exp_diff' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
40 --tss_groups_exp_output=$tss_groups_exp_diff
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
41 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
42 #if 'cds_exp_fpkm_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
43 --cds_exp_fpkm_tracking_output=$cds_exp_fpkm_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
44 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
45 #if 'splicing_diff' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
46 --splicing_diff_output=$splicing_diff
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
47 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
48 #if 'cds_diff' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
49 --cds_diff_output=$cds_diff
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
50 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
51 #if 'promoters_diff' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
52 --promoters_diff_output=$promoters_diff
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
53 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
54 #if 'cds_read_group_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
55 --cds_read_group_tracking=$cds_read_group_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
56 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
57 #if 'tss_groups_read_group_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
58 --tss_groups_read_group_tracking=$tss_groups_read_group_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
59 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
60 #if 'genes_read_group_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
61 --genes_read_group_tracking=$genes_read_group_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
62 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
63 #if 'isoforms_read_group_tracking' in $sel_outputs:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
64 --isoforms_read_group_tracking=$isoforms_read_group_tracking
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
65 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
66
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
67 ## Set advanced data parameters?
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
68 #if $additional.sAdditional == "Yes":
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
69 -m $additional.frag_mean_len
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
70 -s $additional.frag_len_std_dev
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
71 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
72
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
73 ## Normalization?
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
74 #if str($do_normalization) == "Yes":
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
75 -N
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
76 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
77
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
78 ## Multi-read correct?
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
79 #if str($multiread_correct) == "Yes":
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
80 -u
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
81 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
82
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
83 ## Bias correction?
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
84 #if $bias_correction.do_bias_correction == "Yes":
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
85 -b
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
86 #if $bias_correction.seq_source.index_source == "history":
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
87 --ref_file=$bias_correction.seq_source.ref_file
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
88 #else:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
89 --ref_file="None"
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
90 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
91 --dbkey=${gtf_input.metadata.dbkey}
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
92 --index_dir=${GALAXY_DATA_INDEX_DIR}
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
93 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
94
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
95 ## Inputs.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
96 --inputA=$gtf_input
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
97 #if $group_analysis.do_groups == "No":
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
98 --input1=$aligned_reads1
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
99 --input2=$aligned_reads2
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
100 #else:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
101 ## Replicates.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
102 --labels
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
103 #for $group in $group_analysis.groups
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
104 ${group.group}
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
105 #end for
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
106 --files
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
107 #for $group in $group_analysis.groups
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
108 #for $file in $group.files:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
109 ${file.file}
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
110 #end for
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
111 ,
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
112 #end for
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
113 #end if
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
114
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
115 </command>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
116 <inputs>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
117 <param format="gtf,gff3" name="gtf_input" type="data" label="Transcripts" help="A transcript GFF3 or GTF file produced by cufflinks, cuffcompare, or other source."/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
118 <conditional name="group_analysis">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
119 <param name="do_groups" type="select" label="Perform replicate analysis" help="Perform cuffdiff with replicates in each group.">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
120 <option value="No">No</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
121 <option value="Yes">Yes</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
122 </param>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
123 <when value="Yes">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
124 <repeat name="groups" title="Group">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
125 <param name="group" title="Group name" type="text" label="Group name (no spaces or commas)"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
126 <repeat name="files" title="Replicate">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
127 <param name="file" label="Add file" type="data" format="sam,bam"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
128 </repeat>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
129 </repeat>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
130 </when>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
131 <when value="No">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
132 <param format="sam,bam" name="aligned_reads1" type="data" label="SAM or BAM file of aligned RNA-Seq reads" help=""/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
133 <param format="sam,bam" name="aligned_reads2" type="data" label="SAM or BAM file of aligned RNA-Seq reads" help=""/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
134 </when>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
135 </conditional>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
136
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
137 <param name="fdr" type="float" value="0.05" label="False Discovery Rate" help="The allowed false discovery rate."/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
138
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
139 <param name="min_alignment_count" type="integer" value="10" label="Min Alignment Count" help="The minimum number of alignments in a locus for needed to conduct significance testing on changes in that locus observed between samples."/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
140
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
141 <param name="do_normalization" type="select" label="Perform quartile normalization" help="Removes top 25% of genes from FPKM denominator to improve accuracy of differential expression calls for low abundance transcripts.">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
142 <option value="No">No</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
143 <option value="Yes">Yes</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
144 </param>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
145
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
146 <param name="multiread_correct" type="select" label="Use multi-read correct" help="Tells Cufflinks to do an initial estimation procedure to more accurately weight reads mapping to multiple locations in the genome.">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
147 <option value="No" selected="true">No</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
148 <option value="Yes">Yes</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
149 </param>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
150
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
151 <conditional name="bias_correction">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
152 <param name="do_bias_correction" type="select" label="Perform Bias Correction" help="Bias detection and correction can significantly improve accuracy of transcript abundance estimates.">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
153 <option value="No">No</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
154 <option value="Yes">Yes</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
155 </param>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
156 <when value="Yes">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
157 <conditional name="seq_source">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
158 <param name="index_source" type="select" label="Reference sequence data">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
159 <option value="cached">Locally cached</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
160 <option value="history">History</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
161 </param>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
162 <when value="cached"></when>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
163 <when value="history">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
164 <param name="ref_file" type="data" format="fasta" label="Using reference file" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
165 </when>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
166 </conditional>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
167 </when>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
168 <when value="No"></when>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
169 </conditional>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
170
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
171 <conditional name="additional">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
172 <param name="sAdditional" type="select" label="Set Additional Parameters? (not recommended)">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
173 <option value="No">No</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
174 <option value="Yes">Yes</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
175 </param>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
176 <when value="No"></when>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
177 <when value="Yes">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
178 <param name="frag_mean_len" type="integer" value="200" label="Average Fragment Length"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
179 <param name="frag_len_std_dev" type="integer" value="80" label="Fragment Length Standard Deviation"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
180 </when>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
181 </conditional>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
182
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
183 <param name="output_sel" type="select" multiple="true" display="checkboxes" force_select="true" label="Select outputs for history datasets">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
184 <option value="cuffdata">cuffdata - html page with links to cuffdiff outputs</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
185 <option value="cummeRbund_db">cummeRbund database</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
186 <option value="run_info">run.info</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
187 <option value="read_groups_info">read_groups.info</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
188 <option value="splicing_diff">splicing.diff</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
189 <option value="promoters_diff">promoters.diff</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
190 <option value="genes_exp_diff">genes_exp.diff</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
191 <option value="genes_fpkm_tracking">genes.fpkm_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
192 <option value="genes_count_tracking">genes.count_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
193 <option value="genes_read_group_tracking">genes.read_group_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
194 <option value="isoforms_exp_diff">isoforms.exp_diff</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
195 <option value="isoforms_fpkm_tracking">isoforms.fpkm_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
196 <option value="isoforms_count_tracking">isoforms.count_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
197 <option value="isoforms_read_group_tracking">isoforms.read_group_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
198 <option value="cds_diff">cds.diff</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
199 <option value="cds_exp_diff">cds_exp.diff</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
200 <option value="cds_fpkm_tracking">cds.fpkm_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
201 <option value="cds_count_tracking">cds.count_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
202 <option value="cds_read_group_tracking">cds.read_group_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
203 <option value="tss_groups_exp_diff">tss_groups_exp.diff</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
204 <option value="tss_groups_fpkm_tracking">tss_groups.fpkm_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
205 <option value="tss_groups_count_tracking">tss_groups.count_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
206 <option value="tss_groups_read_group_tracking">tss_groups.read_group_tracking</option>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
207 </param>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
208
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
209 </inputs>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
210
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
211 <outputs>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
212 <data format="text" name="run_info" label="${tool.name} on ${on_string}: run.info">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
213 <filter>output_sel and 'run_info' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
214 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
215 <data format="tabular" name="read_groups_info" label="${tool.name} on ${on_string}: read_groups.info">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
216 <filter>output_sel and 'read_groups_info' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
217 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
218 <data format="tabular" name="splicing_diff" label="${tool.name} on ${on_string}: splicing differential expression testing">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
219 <filter>output_sel and 'splicing_diff' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
220 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
221 <data format="tabular" name="promoters_diff" label="${tool.name} on ${on_string}: promoters differential expression testing">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
222 <filter>output_sel and 'promoters_diff' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
223 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
224 <data format="tabular" name="cds_diff" label="${tool.name} on ${on_string}: CDS overloading diffential expression testing">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
225 <filter>output_sel and 'cds_diff' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
226 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
227 <data format="tabular" name="cds_exp_diff" label="${tool.name} on ${on_string}: CDS differential expression testing">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
228 <filter>output_sel and 'cds_exp_diff' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
229 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
230 <data format="tabular" name="cds_fpkm_tracking" label="${tool.name} on ${on_string}: CDS FPKM tracking">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
231 <filter>output_sel and 'cds_fpkm_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
232 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
233 <data format="tabular" name="cds_count_tracking" label="${tool.name} on ${on_string}: CDS counts">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
234 <filter>output_sel and 'cds_count_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
235 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
236 <data format="tabular" name="cds_read_group_tracking" label="${tool.name} on ${on_string}: CDS Read Group tracking">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
237 <filter>output_sel and 'cds_read_group_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
238 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
239 <data format="tabular" name="tss_groups_exp_diff" label="${tool.name} on ${on_string}: TSS groups differential expression testing">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
240 <filter>output_sel and 'tss_groups_exp_diff' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
241 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
242 <data format="tabular" name="tss_groups_fpkm_tracking" label="${tool.name} on ${on_string}: TSS groups FPKM tracking">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
243 <filter>output_sel and 'tss_groups_fpkm_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
244 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
245 <data format="tabular" name="tss_groups_count_tracking" label="${tool.name} on ${on_string}: TSS groups counts">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
246 <filter>output_sel and 'tss_groups_count_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
247 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
248 <data format="tabular" name="tss_groups_read_group_tracking" label="${tool.name} on ${on_string}: TSS groups Read Group tracking">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
249 <filter>output_sel and 'tss_groups_read_group_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
250 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
251 <data format="tabular" name="isoforms_exp_diff" label="${tool.name} on ${on_string}: transcript differential expression testing">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
252 <filter>output_sel and 'isoforms_exp_diff' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
253 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
254 <data format="tabular" name="isoforms_fpkm_tracking" label="${tool.name} on ${on_string}: transcript FPKM tracking">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
255 <filter>output_sel and 'isoforms_fpkm_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
256 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
257 <data format="tabular" name="isoforms_count_tracking" label="${tool.name} on ${on_string}: transcript counts">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
258 <filter>output_sel and 'isoforms_count_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
259 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
260 <data format="tabular" name="isoforms_read_group_tracking" label="${tool.name} on ${on_string}: transcript Read Group tracking">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
261 <filter>output_sel and 'isoforms_read_group_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
262 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
263 <data format="tabular" name="genes_exp_diff" label="${tool.name} on ${on_string}: gene differential expression testing">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
264 <filter>output_sel and 'genes_exp_diff' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
265 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
266 <data format="tabular" name="genes_fpkm_tracking" label="${tool.name} on ${on_string}: gene FPKM tracking">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
267 <filter>output_sel and 'genes_fpkm_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
268 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
269 <data format="tabular" name="genes_count_tracking" label="${tool.name} on ${on_string}: gene counts">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
270 <filter>output_sel and 'genes_count_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
271 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
272 <data format="tabular" name="genes_read_group_tracking" label="${tool.name} on ${on_string}: gene Read Group tracking">
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
273 <filter>output_sel and 'genes_read_group_tracking' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
274 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
275 <data format="cuffdata" name="cuffdata" label="${tool.name} on ${on_string}: cuffdata" >
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
276 <filter>not output_sel or output_sel and 'cuffdata' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
277 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
278 <data format="cuffdatadb" name="cummeRbund_db" label="${tool.name} on ${on_string}: cummeRbund sqlite Database" >
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
279 <filter>output_sel and 'cummeRbund_db' in output_sel</filter>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
280 </data>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
281 </outputs>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
282 <stdio>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
283 <exit_code range="1:" level="fatal" description="Cufflinks Err" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
284 </stdio>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
285
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
286
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
287 <tests>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
288 <test>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
289 <!--
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
290 cuffdiff cuffcompare_out5.gtf cuffdiff_in1.sam cuffdiff_in2.sam
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
291 -->
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
292 <param name="gtf_input" value="cuffcompare_out5.gtf" ftype="gtf" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
293 <param name="do_groups" value="No" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
294 <param name="aligned_reads1" value="cuffdiff_in1.sam" ftype="sam" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
295 <param name="aligned_reads2" value="cuffdiff_in2.sam" ftype="sam" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
296 <!-- Defaults. -->
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
297 <param name="fdr" value="0.05" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
298 <param name="min_alignment_count" value="0" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
299 <param name="do_bias_correction" value="No" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
300 <param name="do_normalization" value="No" />
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
301 <param name="multiread_correct" value="No"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
302 <param name="sAdditional" value="No"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
303 <!--
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
304 Line diffs are needed because cuffdiff does not produce deterministic output.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
305 TODO: can we find datasets that lead to deterministic behavior?
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
306 -->
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
307 <output name="splicing_diff" file="cuffdiff_out9.txt"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
308 <output name="promoters_diff" file="cuffdiff_out10.txt"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
309 <output name="cds_diff" file="cuffdiff_out11.txt"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
310 <output name="cds_exp_fpkm_tracking" file="cuffdiff_out4.txt"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
311 <output name="cds_fpkm_tracking" file="cuffdiff_out8.txt"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
312 <output name="tss_groups_exp" file="cuffdiff_out3.txt" lines_diff="200"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
313 <output name="tss_groups_fpkm_tracking" file="cuffdiff_out7.txt"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
314 <output name="genes_exp" file="cuffdiff_out2.txt" lines_diff="200"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
315 <output name="genes_fpkm_tracking" file="cuffdiff_out6.txt" lines_diff="200"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
316 <output name="isoforms_exp" file="cuffdiff_out1.txt" lines_diff="200"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
317 <output name="isoforms_fpkm_tracking" file="cuffdiff_out5.txt" lines_diff="200"/>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
318 </test>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
319 </tests>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
320
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
321 <help>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
322 **Cuffdiff Overview**
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
323
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
324 Cuffdiff is part of Cufflinks_. Cuffdiff find significant changes in transcript expression, splicing, and promoter use. Please cite: Trapnell C, Williams BA, Pertea G, Mortazavi AM, Kwan G, van Baren MJ, Salzberg SL, Wold B, Pachter L. Transcript assembly and abundance estimation from RNA-Seq reveals thousands of new transcripts and switching among isoforms. Nature Biotechnology doi:10.1038/nbt.1621
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
325
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
326 .. _Cufflinks: http://cufflinks.cbcb.umd.edu/
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
327
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
328 ------
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
329
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
330 **Know what you are doing**
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
331
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
332 .. class:: warningmark
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
333
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
334 There is no such thing (yet) as an automated gearshift in expression analysis. It is all like stick-shift driving in San Francisco. In other words, running this tool with default parameters will probably not give you meaningful results. A way to deal with this is to **understand** the parameters by carefully reading the `documentation`__ and experimenting. Fortunately, Galaxy makes experimenting easy.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
335
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
336 .. __: http://cufflinks.cbcb.umd.edu/manual.html#cuffdiff
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
337
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
338 ------
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
339
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
340 **Input format**
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
341
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
342 Cuffdiff takes Cufflinks or Cuffcompare GTF files as input along with two SAM files containing the fragment alignments for two or more samples.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
343
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
344 ------
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
345
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
346 **Outputs**
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
347
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
348 Cuffdiff produces many output files:
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
349
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
350 1. Transcript FPKM expression tracking.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
351 2. Gene FPKM expression tracking; tracks the summed FPKM of transcripts sharing each gene_id
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
352 3. Primary transcript FPKM tracking; tracks the summed FPKM of transcripts sharing each tss_id
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
353 4. Coding sequence FPKM tracking; tracks the summed FPKM of transcripts sharing each p_id, independent of tss_id
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
354 5. Transcript differential FPKM.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
355 6. Gene differential FPKM. Tests difference sin the summed FPKM of transcripts sharing each gene_id
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
356 7. Primary transcript differential FPKM. Tests difference sin the summed FPKM of transcripts sharing each tss_id
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
357 8. Coding sequence differential FPKM. Tests difference sin the summed FPKM of transcripts sharing each p_id independent of tss_id
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
358 9. Differential splicing tests: this tab delimited file lists, for each primary transcript, the amount of overloading detected among its isoforms, i.e. how much differential splicing exists between isoforms processed from a single primary transcript. Only primary transcripts from which two or more isoforms are spliced are listed in this file.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
359 10. Differential promoter tests: this tab delimited file lists, for each gene, the amount of overloading detected among its primary transcripts, i.e. how much differential promoter use exists between samples. Only genes producing two or more distinct primary transcripts (i.e. multi-promoter genes) are listed here.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
360 11. Differential CDS tests: this tab delimited file lists, for each gene, the amount of overloading detected among its coding sequences, i.e. how much differential CDS output exists between samples. Only genes producing two or more distinct CDS (i.e. multi-protein genes) are listed here.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
361
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
362 -------
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
363
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
364 **Settings**
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
365
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
366 All of the options have a default value. You can change any of them. Most of the options in Cuffdiff have been implemented here.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
367
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
368 ------
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
369
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
370 **Cuffdiff parameter list**
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
371
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
372 This is a list of implemented Cuffdiff options::
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
373
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
374 -m INT Average fragement length; default 200
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
375 -s INT Fragment legnth standard deviation; default 80
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
376 -c INT The minimum number of alignments in a locus for needed to conduct significance testing on changes in that locus observed between samples. If no testing is performed, changes in the locus are deemed not significant, and the locus' observed changes don't contribute to correction for multiple testing. The default is 1,000 fragment alignments (up to 2,000 paired reads).
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
377 --FDR FLOAT The allowed false discovery rate. The default is 0.05.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
378 --num-importance-samples INT Sets the number of importance samples generated for each locus during abundance estimation. Default: 1000
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
379 --max-mle-iterations INT Sets the number of iterations allowed during maximum likelihood estimation of abundances. Default: 5000
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
380 -N With this option, Cufflinks excludes the contribution of the top 25 percent most highly expressed genes from the number of mapped fragments used in the FPKM denominator. This can improve robustness of differential expression calls for less abundant genes and transcripts.
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
381
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
382 </help>
da7241f92ecf Uploaded
jjohnson
parents:
diff changeset
383 </tool>