annotate tophat2_wrapper.xml @ 0:ffa30bedbee3

Imported from capsule None
author devteam
date Mon, 27 Jan 2014 09:26:57 -0500
parents
children ae06af1118dc
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
1 <tool id="tophat2" name="Tophat2" version="0.6">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
2 <!-- Wrapper compatible with Tophat version 2.0.0+ -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
3 <description>Gapped-read mapper for RNA-seq data</description>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
4 <version_command>tophat2 --version</version_command>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
5 <requirements>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
6 <requirement type="package" version="0.1.18">samtools</requirement>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
7 <requirement type="package" version="2.1.0">bowtie2</requirement>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
8 <requirement type="package" version="2.0.9">tophat2</requirement>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
9 </requirements>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
10
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
11 <command>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
12 ##
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
13 ## Set path to index, building the reference if necessary.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
14 ##
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
15
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
16 #set index_path = ''
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
17 #if $refGenomeSource.genomeSource == "history":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
18 bowtie2-build "$refGenomeSource.ownFile" genome ; ln -s "$refGenomeSource.ownFile" genome.fa ;
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
19 #set index_path = 'genome'
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
20 #else:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
21 #set index_path = $refGenomeSource.index.fields.path
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
22 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
23
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
24 ##
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
25 ## Run tophat.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
26 ##
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
27
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
28 tophat2
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
29
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
30 ## Change this to accommodate the number of threads you have available.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
31 --num-threads \${GALAXY_SLOTS:-4}
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
32
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
33 ## Set params.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
34 #if $params.settingsType == "full":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
35 --read-mismatches $params.read_mismatches
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
36 #if str($params.bowtie_n) == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
37 --bowtie-n
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
38 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
39
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
40 --read-edit-dist $params.read_edit_dist
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
41 --read-realign-edit-dist $params.read_realign_edit_dist
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
42 -a $params.anchor_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
43 -m $params.splice_mismatches
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
44 -i $params.min_intron_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
45 -I $params.max_intron_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
46 -g $params.max_multihits
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
47 --min-segment-intron $params.min_segment_intron
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
48 --max-segment-intron $params.max_segment_intron
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
49 --segment-mismatches $params.seg_mismatches
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
50 --segment-length $params.seg_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
51 --library-type $params.library_type
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
52
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
53 ## Indel search.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
54 #if $params.indel_search.allow_indel_search == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
55 ## --allow-indels
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
56 --max-insertion-length $params.indel_search.max_insertion_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
57 --max-deletion-length $params.indel_search.max_deletion_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
58 #else:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
59 --no-novel-indels
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
60 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
61
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
62 ## Supplying junctions parameters.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
63 #if $params.own_junctions.use_junctions == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
64 #if $params.own_junctions.gene_model_ann.use_annotations == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
65 -G $params.own_junctions.gene_model_ann.gene_annotation_model
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
66 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
67 #if $params.own_junctions.raw_juncs.use_juncs == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
68 -j $params.own_junctions.raw_juncs.raw_juncs
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
69 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
70 #if str($params.own_junctions.no_novel_juncs) == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
71 --no-novel-juncs
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
72 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
73 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
74
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
75 #if $params.coverage_search.use_search == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
76 --coverage-search
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
77 --min-coverage-intron $params.coverage_search.min_coverage_intron
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
78 --max-coverage-intron $params.coverage_search.max_coverage_intron
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
79 #else:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
80 --no-coverage-search
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
81 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
82
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
83 #if str($params.microexon_search) == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
84 --microexon-search
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
85 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
86
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
87 #if $params.fusion_search.do_search == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
88 --fusion-search
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
89 --fusion-anchor-length $params.fusion_search.anchor_len
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
90 --fusion-min-dist $params.fusion_search.min_dist
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
91 --fusion-read-mismatches $params.fusion_search.read_mismatches
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
92 --fusion-multireads $params.fusion_search.multireads
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
93 --fusion-multipairs $params.fusion_search.multipairs
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
94 --fusion-ignore-chromosomes "$params.fusion_search.ignore_chromosomes"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
95 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
96
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
97 #if $params.bowtie2_settings.b2_settings == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
98 #if $params.bowtie2_settings.preset.b2_preset == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
99 --b2-$params.bowtie2_settings.preset.b2_preset_select
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
100 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
101 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
102
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
103 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
104
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
105 ## Read group information.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
106 #if $readGroup.specReadGroup == "yes"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
107 --rg-id "$readGroup.rgid"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
108 --rg-library "$readGroup.rglb"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
109 --rg-platform "$readGroup.rgpl"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
110 --rg-sample "$readGroup.rgsm"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
111 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
112
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
113 ## Set index path, inputs and parameters specific to paired data.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
114 #if $singlePaired.sPaired == "paired"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
115 -r $singlePaired.mate_inner_distance
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
116 --mate-std-dev=$singlePaired.mate_std_dev
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
117
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
118 #if str($singlePaired.report_discordant_pairs) == "No":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
119 --no-discordant
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
120 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
121
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
122 ${index_path} $singlePaired.input1 $singlePaired.input2
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
123 #else
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
124 ${index_path} $singlePaired.input1
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
125 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
126 </command>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
127
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
128 <inputs>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
129 <conditional name="singlePaired">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
130 <param name="sPaired" type="select" label="Is this library mate-paired?">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
131 <option value="single">Single-end</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
132 <option value="paired">Paired-end</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
133 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
134 <when value="single">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
135 <param format="fastqsanger" name="input1" type="data" label="RNA-Seq FASTQ file" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
136 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
137 <when value="paired">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
138 <param format="fastqsanger" name="input1" type="data" label="RNA-Seq FASTQ file, forward reads" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
139 <param format="fastqsanger" name="input2" type="data" label="RNA-Seq FASTQ file, reverse reads" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
140 <param name="mate_inner_distance" type="integer" value="300" label="Mean Inner Distance between Mate Pairs" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
141 <param name="mate_std_dev" type="integer" value="20" label="Std. Dev for Distance between Mate Pairs" help="The standard deviation for the distribution on inner distances between mate pairs."/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
142 <!-- Discordant pairs. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
143 <param name="report_discordant_pairs" type="select" label="Report discordant pair alignments?">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
144 <option value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
145 <option selected="True" value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
146 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
147 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
148 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
149 <expand macro="refGenomeSourceConditional">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
150 <options from_data_table="tophat2_indexes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
151 <filter type="sort_by" column="2"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
152 <validator type="no_options" message="No genomes are available for the selected input dataset"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
153 </options>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
154 </expand>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
155 <conditional name="params">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
156 <param name="settingsType" type="select" label="TopHat settings to use" help="You can use the default settings or set custom values for any of Tophat's parameters.">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
157 <option value="preSet">Use Defaults</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
158 <option value="full">Full parameter list</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
159 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
160 <when value="preSet" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
161 <!-- Full/advanced params. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
162 <when value="full">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
163 <param name="read_realign_edit_dist" type="integer" value="1000" label="Max realign edit distance" help="Some of the reads spanning multiple exons may be mapped incorrectly as a contiguous alignment to the genome even though the correct alignment should be a spliced one - this can happen in the presence of processed pseudogenes that are rarely (if at all) transcribed or expressed. This option can direct TopHat to re-align reads for which the edit distance of an alignment obtained in a previous mapping step is above or equal to this option value. If you set this option to 0, TopHat will map every read in all the mapping steps (transcriptome if you provided gene annotations, genome, and finally splice variants detected by TopHat), reporting the best possible alignment found in any of these mapping steps. This may greatly increase the mapping accuracy at the expense of an increase in running time. The default value for this option is set such that TopHat will not try to realign reads already mapped in earlier steps." />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
164
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
165 <param name="read_edit_dist" type="integer" value="2" label="Max edit distance" help="Final read alignments having more than these many edit distance are discarded." />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
166
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
167 <param name="library_type" type="select" label="Library Type" help="TopHat will treat the reads as strand specific. Every read alignment will have an XS attribute tag. Consider supplying library type options below to select the correct RNA-seq protocol.">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
168 <option value="fr-unstranded">FR Unstranded</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
169 <option value="fr-firststrand">FR First Strand</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
170 <option value="fr-secondstrand">FR Second Strand</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
171 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
172 <param name="read_mismatches" type="integer" value="2" label="Final read mismatches" help="Final read alignments having more than these many mismatches are discarded." />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
173 <param name="bowtie_n" type="select" label="Use bowtie -n mode">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
174 <option selected="true" value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
175 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
176 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
177 <param name="anchor_length" type="integer" value="8" label="Anchor length (at least 3)" help="Report junctions spanned by reads with at least this many bases on each side of the junction." />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
178 <param name="splice_mismatches" type="integer" value="0" label="Maximum number of mismatches that can appear in the anchor region of spliced alignment" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
179 <param name="min_intron_length" type="integer" value="70" label="The minimum intron length" help="TopHat will ignore donor/acceptor pairs closer than this many bases apart." />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
180 <param name="max_intron_length" type="integer" value="500000" label="The maximum intron length" help="When searching for junctions ab initio, TopHat will ignore donor/acceptor pairs farther than this many bases apart, except when such a pair is supported by a split segment alignment of a long read." />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
181 <expand macro="indel_searchConditional" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
182 alignments (number of reads divided by average depth of coverage)" help="0.0 to 1.0 (0 to turn off)" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
183 <param name="max_multihits" type="integer" value="20" label="Maximum number of alignments to be allowed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
184 <param name="min_segment_intron" type="integer" value="50" label="Minimum intron length that may be found during split-segment (default) search" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
185 <param name="max_segment_intron" type="integer" value="500000" label="Maximum intron length that may be found during split-segment (default) search" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
186 <param name="seg_mismatches" type="integer" min="0" max="3" value="2" label="Number of mismatches allowed in each segment alignment for reads mapped independently" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
187 <param name="seg_length" type="integer" value="25" label="Minimum length of read segments" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
188
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
189 <!-- Options for supplying own junctions. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
190 <expand macro="own_junctionsConditional" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
191 <!-- Coverage search. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
192 <conditional name="coverage_search">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
193 <param name="use_search" type="select" label="Use Coverage Search" help="Enables the coverage based search for junctions. Use when coverage search is disabled by default (such as for reads 75bp or longer), for maximum sensitivity.">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
194 <option selected="true" value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
195 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
196 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
197 <when value="Yes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
198 <param name="min_coverage_intron" type="integer" value="50" label="Minimum intron length that may be found during coverage search" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
199 <param name="max_coverage_intron" type="integer" value="20000" label="Maximum intron length that may be found during coverage search" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
200 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
201 <when value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
202 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
203
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
204 <!-- Microexon search params -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
205 <param name="microexon_search" type="select" label="Use Microexon Search" help="With this option, the pipeline will attempt to find alignments incident to microexons. Works only for reads 50bp or longer.">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
206 <option value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
207 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
208 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
209
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
210 <!-- Fusion mapping. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
211 <conditional name="fusion_search">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
212 <param name="do_search" type="select" label="Do Fusion Search">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
213 <option selected="true" value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
214 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
215 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
216 <when value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
217 <when value="Yes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
218 <param name="anchor_len" type="integer" value="20" label="Anchor Length" help="A 'supporting' read must map to both sides of a fusion by at least this many bases."/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
219 <param name="min_dist" type="integer" value="10000000" label="Minimum Distance" help="For intra-chromosomal fusions, TopHat-Fusion tries to find fusions separated by at least this distance."/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
220 <param name="read_mismatches" type="integer" value="2" label="Read Mismatches" help="Reads support fusions if they map across fusion with at most this many mismatches."/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
221 <param name="multireads" type="integer" value="2" label="Multireads" help="Reads that map to more than this many places will be ignored. It may be possible that a fusion is supported by reads (or pairs) that map to multiple places."/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
222 <param name="multipairs" type="integer" value="2" label="Multipairs" help="Pairs that map to more than this many places will be ignored."/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
223 <param name="ignore_chromosomes" type="text" value='' label="Ignore some chromosomes such as chrM when detecting fusion break points"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
224 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
225 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
226
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
227 <!-- Bowtie2 settings. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
228 <conditional name="bowtie2_settings">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
229 <param name="b2_settings" type="select" label="Set Bowtie2 settings">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
230 <option selected="true" value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
231 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
232 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
233 <when value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
234 <when value="Yes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
235 <conditional name="preset">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
236 <param name="b2_preset" type="select" label="Use Preset options">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
237 <option selected="true" value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
238 <option value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
239 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
240 <when value="Yes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
241 <param name="b2_preset_select" type="select" label="Preset option">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
242 <option value="very-fast">Very fast</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
243 <option value="fast">Fast</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
244 <option selected="true" value="sensitive">Sensitive</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
245 <option value="very-sensitive">Very sensitive</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
246 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
247 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
248 <!-- TODO: -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
249 <when value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
250 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
251 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
252 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
253 </when> <!-- full -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
254 </conditional> <!-- params -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
255 <conditional name="readGroup">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
256 <param name="specReadGroup" type="select" label="Specify read group?">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
257 <option value="yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
258 <option value="no" selected="True">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
259 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
260 <when value="yes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
261 <param name="rgid" type="text" size="25" label="Read group identifier (ID). Each @RG line must have a unique ID. The value of ID is used in the RG tags of alignment records. Must be unique among all read groups in header section." help="Required if RG specified. Read group IDs may be modified when merging SAM files in order to handle collisions." />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
262 <param name="rglb" type="text" size="25" label="Library name (LB)" help="Required if RG specified" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
263 <param name="rgpl" type="text" size="25" label="Platform/technology used to produce the reads (PL)" help="Required if RG specified. Valid values : CAPILLARY, LS454, ILLUMINA, SOLID, HELICOS, IONTORRENT and PACBIO" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
264 <param name="rgsm" type="text" size="25" label="Sample (SM)" help="Required if RG specified. Use pool name where a pool is being sequenced" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
265 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
266 <when value="no" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
267 </conditional> <!-- readGroup -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
268 </inputs>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
269
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
270 <stdio>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
271 <regex match="Exception|Error" source="both" level="fatal" description="Tool execution failed"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
272 <regex match=".*" source="both" level="log" description="tool progress"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
273 </stdio>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
274
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
275 <outputs>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
276 <data format="txt" name="align_summary" label="${tool.name} on ${on_string}: align_summary" from_work_dir="tophat_out/align_summary.txt"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
277 <data format="tabular" name="fusions" label="${tool.name} on ${on_string}: fusions" from_work_dir="tophat_out/fusions.out">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
278 <filter>(params['settingsType'] == 'full' and params['fusion_search']['do_search'] == 'Yes')</filter>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
279 </data>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
280 <data format="bed" name="insertions" label="${tool.name} on ${on_string}: insertions" from_work_dir="tophat_out/insertions.bed">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
281 <expand macro="dbKeyActions" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
282 </data>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
283 <data format="bed" name="deletions" label="${tool.name} on ${on_string}: deletions" from_work_dir="tophat_out/deletions.bed">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
284 <expand macro="dbKeyActions" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
285 </data>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
286 <data format="bed" name="junctions" label="${tool.name} on ${on_string}: splice junctions" from_work_dir="tophat_out/junctions.bed">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
287 <expand macro="dbKeyActions" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
288 </data>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
289 <data format="bam" name="accepted_hits" label="${tool.name} on ${on_string}: accepted_hits" from_work_dir="tophat_out/accepted_hits.bam">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
290 <expand macro="dbKeyActions" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
291 </data>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
292 </outputs>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
293
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
294 <macros>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
295 <import>tophat_macros.xml</import>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
296 <macro name="dbKeyActions">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
297 <actions>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
298 <conditional name="refGenomeSource.genomeSource">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
299 <when value="indexed">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
300 <action type="metadata" name="dbkey">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
301 <option type="from_data_table" name="tophat2_indexes" column="1" offset="0">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
302 <filter type="param_value" column="0" value="#" compare="startswith" keep="False"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
303 <filter type="param_value" ref="refGenomeSource.index" column="0"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
304 </option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
305 </action>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
306 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
307 <when value="history">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
308 <action type="metadata" name="dbkey">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
309 <option type="from_param" name="refGenomeSource.ownFile" param_attribute="dbkey" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
310 </action>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
311 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
312 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
313 </actions>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
314 </macro>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
315 </macros>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
316
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
317 <tests>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
318 <!-- Test base-space single-end reads with pre-built index and preset parameters -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
319 <test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
320 <!-- TopHat commands:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
321 tophat2 -o tmp_dir -p 1 tophat_in1 test-data/tophat_in2.fastqsanger
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
322 Rename the files in tmp_dir appropriately
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
323 -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
324 <param name="sPaired" value="single" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
325 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
326 <param name="genomeSource" value="indexed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
327 <param name="index" value="tophat_test" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
328 <param name="settingsType" value="preSet" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
329 <param name="specReadGroup" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
330 <output name="junctions" file="tophat2_out1j.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
331 <output name="accepted_hits" file="tophat_out1h.bam" compare="sim_size" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
332 </test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
333 <!-- Test using base-space test data: paired-end reads, index from history. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
334 <test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
335 <!-- TopHat commands:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
336 bowtie2-build -f test-data/tophat_in1.fasta tophat_in1
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
337 tophat2 -o tmp_dir -p 1 -r 20 tophat_in1 test-data/tophat_in2.fastqsanger test-data/tophat_in3.fastqsanger
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
338 Rename the files in tmp_dir appropriately
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
339 -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
340 <param name="sPaired" value="paired" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
341 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
342 <param name="input2" ftype="fastqsanger" value="tophat_in3.fastqsanger" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
343 <param name="genomeSource" value="history" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
344 <param name="ownFile" ftype="fasta" value="tophat_in1.fasta" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
345 <param name="mate_inner_distance" value="20" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
346 <param name="settingsType" value="preSet" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
347 <param name="specReadGroup" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
348 <output name="junctions" file="tophat2_out2j.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
349 <output name="accepted_hits" file="tophat_out2h.bam" compare="sim_size" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
350 </test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
351 <!-- Test base-space single-end reads with user-supplied reference fasta and full parameters -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
352 <test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
353 <!-- Tophat commands:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
354 bowtie2-build -f test-data/tophat_in1.fasta tophat_in1
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
355 tophat2 -o tmp_dir -p 1 -a 8 -m 0 -i 70 -I 500000 -g 40 +coverage-search +min-coverage-intron 50 +max-coverage-intro 20000 +segment-mismatches 2 +segment-length 25 +microexon-search tophat_in1 test-data/tophat_in2.fastqsanger
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
356 Replace the + with double-dash
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
357 Rename the files in tmp_dir appropriately
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
358 -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
359 <param name="sPaired" value="single"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
360 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
361 <param name="genomeSource" value="history"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
362 <param name="ownFile" value="tophat_in1.fasta"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
363 <param name="settingsType" value="full"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
364 <param name="library_type" value="FR Unstranded"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
365 <param name="read_mismatches" value="2"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
366 <param name="bowtie_n" value="No"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
367 <param name="anchor_length" value="8"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
368 <param name="splice_mismatches" value="0"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
369 <param name="min_intron_length" value="70"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
370 <param name="max_intron_length" value="500000"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
371 <param name="max_multihits" value="40"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
372 <param name="min_segment_intron" value="50" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
373 <param name="max_segment_intron" value="500000" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
374 <param name="seg_mismatches" value="2"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
375 <param name="seg_length" value="25"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
376 <param name="allow_indel_search" value="Yes"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
377 <param name="max_insertion_length" value="3"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
378 <param name="max_deletion_length" value="3"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
379 <param name="use_junctions" value="Yes" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
380 <param name="use_annotations" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
381 <param name="use_juncs" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
382 <param name="no_novel_juncs" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
383 <param name="use_search" value="Yes" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
384 <param name="min_coverage_intron" value="50" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
385 <param name="max_coverage_intron" value="20000" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
386 <param name="microexon_search" value="Yes" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
387 <param name="b2_settings" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
388 <!-- Fusion search params -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
389 <param name="do_search" value="Yes" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
390 <param name="anchor_len" value="21" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
391 <param name="min_dist" value="10000021" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
392 <param name="read_mismatches" value="3" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
393 <param name="multireads" value="4" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
394 <param name="multipairs" value="5" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
395 <param name="ignore_chromosomes" value="chrM"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
396 <param name="specReadGroup" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
397 <output name="insertions" file="tophat_out3i.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
398 <output name="deletions" file="tophat_out3d.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
399 <output name="junctions" file="tophat2_out3j.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
400 <output name="accepted_hits" file="tophat_out3h.bam" compare="sim_size" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
401 </test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
402 <!-- Test base-space paired-end reads with user-supplied reference fasta and full parameters -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
403 <test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
404 <!-- TopHat commands:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
405 tophat2 -o tmp_dir -r 20 -p 1 -a 8 -m 0 -i 70 -I 500000 -g 40 +coverage-search +min-coverage-intron 50 +max-coverage-intro 20000 +segment-mismatches 2 +segment-length 25 +microexon-search +report_discordant_pairs tophat_in1 test-data/tophat_in2.fastqsanger test-data/tophat_in3.fastqsanger
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
406 Replace the + with double-dash
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
407 Rename the files in tmp_dir appropriately
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
408 -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
409 <param name="sPaired" value="paired"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
410 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
411 <param name="input2" ftype="fastqsanger" value="tophat_in3.fastqsanger"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
412 <param name="genomeSource" value="indexed"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
413 <param name="index" value="tophat_test"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
414 <param name="mate_inner_distance" value="20"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
415 <param name="settingsType" value="full"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
416 <param name="library_type" value="FR Unstranded"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
417 <param name="read_mismatches" value="5"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
418 <param name="bowtie_n" value="Yes"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
419 <param name="mate_std_dev" value="20"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
420 <param name="anchor_length" value="8"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
421 <param name="splice_mismatches" value="0"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
422 <param name="min_intron_length" value="70"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
423 <param name="max_intron_length" value="500000"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
424 <param name="max_multihits" value="40"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
425 <param name="min_segment_intron" value="50" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
426 <param name="max_segment_intron" value="500000" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
427 <param name="seg_mismatches" value="2"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
428 <param name="seg_length" value="25"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
429 <param name="allow_indel_search" value="No"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
430 <param name="use_junctions" value="Yes" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
431 <param name="use_annotations" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
432 <param name="use_juncs" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
433 <param name="no_novel_juncs" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
434 <param name="report_discordant_pairs" value="Yes" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
435 <param name="use_search" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
436 <param name="microexon_search" value="Yes" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
437 <param name="b2_settings" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
438 <!-- Fusion search params -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
439 <param name="do_search" value="Yes" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
440 <param name="anchor_len" value="21" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
441 <param name="min_dist" value="10000021" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
442 <param name="read_mismatches" value="3" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
443 <param name="multireads" value="4" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
444 <param name="multipairs" value="5" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
445 <param name="ignore_chromosomes" value="chrM"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
446 <param name="specReadGroup" value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
447 <output name="junctions" file="tophat2_out4j.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
448 <output name="accepted_hits" file="tophat_out4h.bam" compare="sim_size" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
449 </test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
450 </tests>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
451
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
452 <help>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
453 **Tophat Overview**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
454
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
455 TopHat_ is a fast splice junction mapper for RNA-Seq reads. It aligns RNA-Seq reads to mammalian-sized genomes using the ultra high-throughput short read aligner Bowtie(2), and then analyzes the mapping results to identify splice junctions between exons. Please cite: Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, and Salzberg SL. TopHat2: accurate alignment
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
456 of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14:R36, 2013.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
457
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
458 .. _Tophat: http://tophat.cbcb.umd.edu/
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
459
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
460 ------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
461
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
462 **Know what you are doing**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
463
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
464 .. class:: warningmark
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
465
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
466 There is no such thing (yet) as an automated gearshift in splice junction identification. 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.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
467
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
468 .. __: http://tophat.cbcb.umd.edu/manual.html
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
469
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
470 ------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
471
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
472 **Input formats**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
473
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
474 Tophat accepts files in Sanger FASTQ format. Use the FASTQ Groomer to prepare your files.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
475
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
476 ------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
477
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
478 **Outputs**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
479
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
480 Tophat produces two output files:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
481
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
482 - junctions -- A UCSC BED_ track of junctions reported by TopHat. Each junction consists of two connected BED blocks, where each block is as long as the maximal overhang of any read spanning the junction. The score is the number of alignments spanning the junction.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
483 - accepted_hits -- A list of read alignments in BAM_ format.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
484
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
485 .. _BED: http://genome.ucsc.edu/FAQ/FAQformat.html#format1
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
486 .. _BAM: http://samtools.sourceforge.net/
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
487
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
488 Two other possible outputs, depending on the options you choose, are insertions and deletions, both of which are in BED format.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
489
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
490 -------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
491
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
492 **Tophat settings**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
493
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
494 All of the options have a default value. You can change any of them. Some of the options in Tophat have been implemented here.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
495
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
496 ------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
497
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
498 **Tophat parameter list**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
499
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
500 This is a list of implemented Tophat options::
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
501
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
502 -r This is the expected (mean) inner distance between mate pairs. For, example, for paired end runs with fragments
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
503 selected at 300bp, where each end is 50bp, you should set -r to be 200. There is no default, and this parameter
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
504 is required for paired end runs.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
505 --mate-std-dev INT The standard deviation for the distribution on inner distances between mate pairs. The default is 20bp.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
506 -a/--min-anchor-length INT The "anchor length". TopHat will report junctions spanned by reads with at least this many bases on each side of the junction. Note that individual spliced
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
507 alignments may span a junction with fewer than this many bases on one side. However, every junction involved in spliced alignments is supported by at least one
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
508 read with this many bases on each side. This must be at least 3 and the default is 8.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
509 -m/--splice-mismatches INT The maximum number of mismatches that may appear in the "anchor" region of a spliced alignment. The default is 0.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
510 -i/--min-intron-length INT The minimum intron length. TopHat will ignore donor/acceptor pairs closer than this many bases apart. The default is 70.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
511 -I/--max-intron-length INT The maximum intron length. When searching for junctions ab initio, TopHat will ignore donor/acceptor pairs farther than this many bases apart, except when such a pair is supported by a split segment alignment of a long read. The default is 500000.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
512 -g/--max-multihits INT Instructs TopHat to allow up to this many alignments to the reference for a given read, and suppresses all alignments for reads with more than this many
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
513 alignments. The default is 40.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
514 -G/--GTF [GTF 2.2 file] Supply TopHat with a list of gene model annotations. TopHat will use the exon records in this file to build a set of known splice junctions for each gene, and will attempt to align reads to these junctions even if they would not normally be covered by the initial mapping.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
515 -j/--raw-juncs [juncs file] Supply TopHat with a list of raw junctions. Junctions are specified one per line, in a tab-delimited format. Records look like: [chrom] [left] [right] [+/-], left and right are zero-based coordinates, and specify the last character of the left sequenced to be spliced to the first character of the right sequence, inclusive.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
516 -no-novel-juncs Only look for junctions indicated in the supplied GFF file. (ignored without -G)
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
517 --no-coverage-search Disables the coverage based search for junctions.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
518 --coverage-search Enables the coverage based search for junctions. Use when coverage search is disabled by default (such as for reads 75bp or longer), for maximum sensitivity.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
519 --microexon-search With this option, the pipeline will attempt to find alignments incident to microexons. Works only for reads 50bp or longer.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
520 --segment-mismatches Read segments are mapped independently, allowing up to this many mismatches in each segment alignment. The default is 2.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
521 --segment-length Each read is cut up into segments, each at least this long. These segments are mapped independently. The default is 25.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
522 --min-coverage-intron The minimum intron length that may be found during coverage search. The default is 50.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
523 --max-coverage-intron The maximum intron length that may be found during coverage search. The default is 20000.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
524 --min-segment-intron The minimum intron length that may be found during split-segment search. The default is 50.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
525 --max-segment-intron The maximum intron length that may be found during split-segment search. The default is 500000.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
526 </help>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
527 </tool>