annotate tophat2_wrapper.xml @ 5:5c5517d2a9e9 draft

planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
author devteam
date Thu, 14 May 2015 09:39:04 -0400
parents 81f97e12e573
children 4eb3c3beb9c7
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
1 <tool id="tophat2" name="Tophat" version="0.9">
0
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>
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
6 <requirement type="package" version="2.2.5">bowtie2</requirement>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
7 <requirement type="package" version="2.0.14">tophat</requirement>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
8 </requirements>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
9
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
10 <command>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
11 ##
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
12 ## Set path to index, building the reference if necessary.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
13 ##
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
14
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
15 #set index_path = ''
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
16 #if $refGenomeSource.genomeSource == "history":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
17 bowtie2-build "$refGenomeSource.ownFile" genome ; ln -s "$refGenomeSource.ownFile" genome.fa ;
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
18 #set index_path = 'genome'
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
19 #else:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
20 #set index_path = $refGenomeSource.index.fields.path
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
21 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
22
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
23 ##
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
24 ## Run tophat.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
25 ##
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
26
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
27 tophat2
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
28
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
29 ## Change this to accommodate the number of threads you have available.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
30 --num-threads \${GALAXY_SLOTS:-4}
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
31
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
32 ## Set params.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
33 #if $params.settingsType == "full":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
34 --read-mismatches $params.read_mismatches
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
35 #if str($params.bowtie_n) == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
36 --bowtie-n
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
37 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
38
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
39 --read-edit-dist $params.read_edit_dist
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
40 --read-realign-edit-dist $params.read_realign_edit_dist
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
41 -a $params.anchor_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
42 -m $params.splice_mismatches
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
43 -i $params.min_intron_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
44 -I $params.max_intron_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
45 -g $params.max_multihits
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
46 --min-segment-intron $params.min_segment_intron
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
47 --max-segment-intron $params.max_segment_intron
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
48 --segment-mismatches $params.seg_mismatches
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
49 --segment-length $params.seg_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
50 --library-type $params.library_type
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
51
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
52 ## Indel search.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
53 #if $params.indel_search.allow_indel_search == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
54 ## --allow-indels
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
55 --max-insertion-length $params.indel_search.max_insertion_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
56 --max-deletion-length $params.indel_search.max_deletion_length
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
57 #else:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
58 --no-novel-indels
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
59 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
60
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
61 ## Supplying junctions parameters.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
62 #if $params.own_junctions.use_junctions == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
63 #if $params.own_junctions.gene_model_ann.use_annotations == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
64 -G $params.own_junctions.gene_model_ann.gene_annotation_model
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
65 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
66 #if $params.own_junctions.raw_juncs.use_juncs == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
67 -j $params.own_junctions.raw_juncs.raw_juncs
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
68 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
69 #if str($params.own_junctions.no_novel_juncs) == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
70 --no-novel-juncs
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
71 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
72 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
73
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
74 #if $params.coverage_search.use_search == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
75 --coverage-search
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
76 --min-coverage-intron $params.coverage_search.min_coverage_intron
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
77 --max-coverage-intron $params.coverage_search.max_coverage_intron
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
78 #else:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
79 --no-coverage-search
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
80 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
81
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
82 #if str($params.microexon_search) == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
83 --microexon-search
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
84 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
85
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
86 #if $params.fusion_search.do_search == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
87 --fusion-search
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
88 --fusion-anchor-length $params.fusion_search.anchor_len
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
89 --fusion-min-dist $params.fusion_search.min_dist
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
90 --fusion-read-mismatches $params.fusion_search.read_mismatches
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
91 --fusion-multireads $params.fusion_search.multireads
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
92 --fusion-multipairs $params.fusion_search.multipairs
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
93 --fusion-ignore-chromosomes "$params.fusion_search.ignore_chromosomes"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
94 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
95
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
96 #if $params.bowtie2_settings.b2_settings == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
97 #if $params.bowtie2_settings.preset.b2_preset == "Yes":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
98 --b2-$params.bowtie2_settings.preset.b2_preset_select
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
99 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
100 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
101
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
102 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
103
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
104 ## Read group information.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
105 #if $readGroup.specReadGroup == "yes"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
106 --rg-id "$readGroup.rgid"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
107 --rg-library "$readGroup.rglb"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
108 --rg-platform "$readGroup.rgpl"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
109 --rg-sample "$readGroup.rgsm"
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
110 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
111
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
112 ## Set index path, inputs and parameters specific to paired data.
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
113 #if $singlePaired.sPaired != "single"
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
114 -r $singlePaired.mate_inner_distance
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
115 --mate-std-dev=$singlePaired.mate_std_dev
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
116
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
117 #if str($singlePaired.report_discordant_pairs) == "No":
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
118 --no-discordant
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
119 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
120
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
121 #if $singlePaired.sPaired == "paired"
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
122 ${index_path} "$singlePaired.input1" "$singlePaired.input2"
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
123 #else
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
124 ${index_path} "$singlePaired.input.forward" "$singlePaired.input.reverse"
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
125 #end if
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
126 #else
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
127 ${index_path} "$singlePaired.input1"
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
128 #end if
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
129 </command>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
130
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
131 <inputs>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
132 <conditional name="singlePaired">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
133 <param name="sPaired" type="select" label="Is this single-end or paired-end data?">
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
134 <option value="single">Single-end</option>
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
135 <option value="paired">Paired-end (as individual datasets)</option>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
136 <option value="paired_collection">Paired-end (as collection)</option>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
137 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
138 <when value="single">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
139 <param format="fastqsanger" name="input1" type="data" label="RNA-Seq FASTQ file" help="Must have Sanger-scaled quality values with ASCII offset 33"/>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
140 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
141 <when value="paired">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
142 <param format="fastqsanger" name="input1" type="data" label="RNA-Seq FASTQ file, forward reads" help="Must have Sanger-scaled quality values with ASCII offset 33" />
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
143 <param format="fastqsanger" name="input2" type="data" label="RNA-Seq FASTQ file, reverse reads" help="Must have Sanger-scaled quality values with ASCII offset 33" />
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
144 <expand macro="paired_parameters" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
145 </when>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
146 <when value="paired_collection">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
147 <param format="fastqsanger" name="input" type="data_collection" collection_type="paired" label="RNA-Seq FASTQ paired reads" help="Must have Sanger-scaled quality values with ASCII offset 33" />
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
148 <expand macro="paired_parameters" />
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
149 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
150 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
151 <expand macro="refGenomeSourceConditional">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
152 <options from_data_table="tophat2_indexes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
153 <filter type="sort_by" column="2"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
154 <validator type="no_options" message="No genomes are available for the selected input dataset"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
155 </options>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
156 </expand>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
157 <conditional name="params">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
158 <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
159 <option value="preSet">Use Defaults</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
160 <option value="full">Full parameter list</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
161 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
162 <when value="preSet" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
163 <!-- Full/advanced params. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
164 <when value="full">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
165 <param name="read_realign_edit_dist" type="integer" value="1000" label="Max realign edit distance" help="--read-realign-edit-dist; 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." />
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
166
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
167 <param name="read_edit_dist" type="integer" value="2" label="Max edit distance" help="--read-edit-dist; Final read alignments having more than these many edit distance are discarded." />
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
168
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
169 <param name="library_type" type="select" label="Library Type" help="--library-type; 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.">
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
170 <option value="fr-unstranded">FR Unstranded</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
171 <option value="fr-firststrand">FR First Strand</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
172 <option value="fr-secondstrand">FR Second Strand</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
173 </param>
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
174 <param name="read_mismatches" type="integer" value="2" label="Final read mismatches" help="--read-mismatches; Final read alignments having more than these many mismatches are discarded." />
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
175 <param name="bowtie_n" type="select" label="Use bowtie -n mode" help="--bowtie-n; TopHat uses &quot;-v&quot; in Bowtie for initial read mapping (the default), but with this option, &quot;-n&quot; is used instead. Read segments are always mapped using &quot;-v&quot; option.">
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
176 <option selected="true" value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
177 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
178 </param>
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
179 <param name="anchor_length" type="integer" value="8" label="Anchor length (at least 3)" help="-a/--min-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 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 read with this many bases on each side. This must be at least 3 and the default is 8." />
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
180 <param name="splice_mismatches" type="integer" value="0" label="Maximum number of mismatches that can appear in the anchor region of spliced alignment" help="-m/--splice-mismatches; The default is 0."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
181 <param name="min_intron_length" type="integer" value="70" label="The minimum intron length" help="-i/--min-intron-length; TopHat will ignore donor/acceptor pairs closer than this many bases apart. The default is 70." />
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
182 <param name="max_intron_length" type="integer" value="500000" label="The maximum intron length" help="-I/--max-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." />
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
183
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
184 <expand macro="indel_searchConditional" />
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
185
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
186 <param name="max_multihits" type="integer" value="20" label="Maximum number of alignments to be allowed" help="-g/--max-multihits; Instructs TopHat to allow up to this many alignments to the reference for a given read, and choose the alignments based on their alignment scores if there are more than this number. The default is 20 for read mapping. Unless you use --report-secondary-alignments, TopHat will report the alignments with the best alignment score. If there are more alignments with the same score than this number, TopHat will randomly report only this many alignments. In case of using --report-secondary-alignments, TopHat will try to report alignments up to this option value, and TopHat may randomly output some of the alignments with the same score to meet this number."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
187 <param name="min_segment_intron" type="integer" value="50" label="Minimum intron length that may be found during split-segment (default) search" help="--min-segment-intron; The minimum intron length that may be found during split-segment search. The default is 50."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
188 <param name="max_segment_intron" type="integer" value="500000" label="Maximum intron length that may be found during split-segment (default) search" help="--max-segment-intron; The maximum intron length that may be found during split-segment search. The default is 500000."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
189 <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" help="--segment-mismatches; Read segments are mapped independently, allowing up to this many mismatches in each segment alignment. The default is 2."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
190 <param name="seg_length" type="integer" value="25" label="Minimum length of read segments" help="--segment-length; Each read is cut up into segments, each at least this long. These segments are mapped independently. The default is 25."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
191 <param name="output_unmapped" type="boolean" truevalue="true" falsevalue="false" checked="false" label="Output unmapped reads" help="If checked, a BAM with the unmapped reads will be added to the history" />
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
192 <!-- Options for supplying own junctions. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
193 <expand macro="own_junctionsConditional" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
194 <!-- Coverage search. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
195 <conditional name="coverage_search">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
196 <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
197 <option selected="true" value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
198 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
199 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
200 <when value="Yes">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
201 <param name="min_coverage_intron" type="integer" value="50" label="Minimum intron length that may be found during coverage search" help="--min-coverage-intron; The minimum intron length that may be found during coverage search. The default is 50."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
202 <param name="max_coverage_intron" type="integer" value="20000" label="Maximum intron length that may be found during coverage search" help="--max-coverage-intron; The maximum intron length that may be found during coverage search. The default is 20000."/>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
203 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
204 <when value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
205 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
206
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
207 <!-- Microexon search params -->
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
208 <param name="microexon_search" type="select" label="Use Microexon Search" help="--microexon-search; With this option, the pipeline will attempt to find alignments incident to microexons. Works only for reads 50bp or longer.">
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
209 <option value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
210 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
211 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
212
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
213 <!-- Fusion mapping. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
214 <conditional name="fusion_search">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
215 <param name="do_search" type="select" label="Do Fusion Search" help="Reads can be aligned to potential fusion transcripts if the --fusion-search option is specified. The fusion alignments are reported in SAM format using custom fields XF and XP (see the output format) and some additional information about fusions will be reported (see fusions.out). Once mapping is done, you can run tophat-fusion-post to filter out fusion transcripts (see the TopHat-Fusion website for more details).">
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
216 <option selected="true" value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
217 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
218 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
219 <when value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
220 <when value="Yes">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
221 <param name="anchor_len" type="integer" value="20" label="Anchor Length" help="--fusion-anchor-length; A 'supporting' read must map to both sides of a fusion by at least this many bases. The default is 20."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
222 <param name="min_dist" type="integer" value="10000000" label="Minimum Distance" help="--fusion-min-dist; For intra-chromosomal fusions, TopHat-Fusion tries to find fusions separated by at least this distance. The default is 10000000."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
223 <param name="read_mismatches" type="integer" value="2" label="Read Mismatches" help="--fusion-read-mismatches; Reads support fusions if they map across fusion with at most this many mismatches. The default is 2."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
224 <param name="multireads" type="integer" value="2" label="Multireads" help="--fusion-multireads; 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. The default is 2."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
225 <param name="multipairs" type="integer" value="2" label="Multipairs" help="--fusion-multipairs; Pairs that map to more than this many places will be ignored. The default is 2."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
226 <param name="ignore_chromosomes" type="text" value='' label="--fusion-ignore-chromosomes; Ignore some chromosomes such as chrM when detecting fusion break points"/>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
227 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
228 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
229
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
230 <!-- Bowtie2 settings. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
231 <conditional name="bowtie2_settings">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
232 <param name="b2_settings" type="select" label="Set Bowtie2 settings">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
233 <option selected="true" value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
234 <option value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
235 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
236 <when value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
237 <when value="Yes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
238 <conditional name="preset">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
239 <param name="b2_preset" type="select" label="Use Preset options">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
240 <option selected="true" value="Yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
241 <option value="No">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
242 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
243 <when value="Yes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
244 <param name="b2_preset_select" type="select" label="Preset option">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
245 <option value="very-fast">Very fast</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
246 <option value="fast">Fast</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
247 <option selected="true" value="sensitive">Sensitive</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
248 <option value="very-sensitive">Very sensitive</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
249 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
250 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
251 <!-- TODO: -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
252 <when value="No" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
253 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
254 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
255 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
256 </when> <!-- full -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
257 </conditional> <!-- params -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
258 <conditional name="readGroup">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
259 <param name="specReadGroup" type="select" label="Specify read group?">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
260 <option value="yes">Yes</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
261 <option value="no" selected="True">No</option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
262 </param>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
263 <when value="yes">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
264 <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
265 <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
266 <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
267 <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
268 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
269 <when value="no" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
270 </conditional> <!-- readGroup -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
271 </inputs>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
272
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
273 <stdio>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
274 <regex match="Exception|Error" source="both" level="fatal" description="Tool execution failed"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
275 <regex match=".*" source="both" level="log" description="tool progress"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
276 </stdio>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
277
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
278 <outputs>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
279 <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
280 <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
281 <filter>(params['settingsType'] == 'full' and params['fusion_search']['do_search'] == 'Yes')</filter>
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="insertions" label="${tool.name} on ${on_string}: insertions" from_work_dir="tophat_out/insertions.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="deletions" label="${tool.name} on ${on_string}: deletions" from_work_dir="tophat_out/deletions.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="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
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 <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
293 <expand macro="dbKeyActions" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
294 </data>
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
295 <data format="bam" name="unmapped" label="${tool.name} on ${on_string}: unmapped" from_work_dir="tophat_out/unmapped.bam">
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
296 <filter>(params['settingsType'] == 'full' and params['output_unmapped'])</filter>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
297 <expand macro="dbKeyActions" />
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
298 </data>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
299
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
300 </outputs>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
301
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
302 <macros>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
303 <import>tophat_macros.xml</import>
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
304 <xml name="paired_parameters">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
305 <param name="mate_inner_distance" type="integer" value="300" label="Mean Inner Distance between Mate Pairs" help="-r/--mate-inner-dist; This is the expected (mean) inner distance between mate pairs. For, example, for paired end runs with fragments selected at 300bp, where each end is 50bp, you should set -r to be 200. The default is 50bp."/>
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
306 <param name="mate_std_dev" type="integer" value="20" label="Std. Dev for Distance between Mate Pairs" help="--mate-std-dev; The standard deviation for the distribution on inner distances between mate pairs. The default is 20bp."/>
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
307 <!-- Discordant pairs. -->
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
308 <param name="report_discordant_pairs" type="select" label="Report discordant pair alignments?" help="--no-discordant">
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
309 <option value="No">No</option>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
310 <option selected="True" value="Yes">Yes</option>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
311 </param>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
312 </xml>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
313 <macro name="dbKeyActions">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
314 <actions>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
315 <conditional name="refGenomeSource.genomeSource">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
316 <when value="indexed">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
317 <action type="metadata" name="dbkey">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
318 <option type="from_data_table" name="tophat2_indexes" column="1" offset="0">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
319 <filter type="param_value" column="0" value="#" compare="startswith" keep="False"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
320 <filter type="param_value" ref="refGenomeSource.index" column="0"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
321 </option>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
322 </action>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
323 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
324 <when value="history">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
325 <action type="metadata" name="dbkey">
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
326 <option type="from_param" name="refGenomeSource.ownFile" param_attribute="dbkey" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
327 </action>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
328 </when>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
329 </conditional>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
330 </actions>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
331 </macro>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
332 </macros>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
333
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
334 <tests>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
335 <!-- Test base-space single-end reads with pre-built index and preset parameters -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
336 <test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
337 <!-- TopHat commands:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
338 tophat2 -o tmp_dir -p 1 tophat_in1 test-data/tophat_in2.fastqsanger
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
339 Rename the files in tmp_dir appropriately
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
340 -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
341 <param name="sPaired" value="single" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
342 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
343 <param name="genomeSource" value="indexed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
344 <param name="index" value="tophat_test" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
345 <param name="settingsType" value="preSet" />
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
346 <param name="specReadGroup" value="no" />
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
347 <output name="junctions" file="tophat2_out1j.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
348 <output name="accepted_hits" file="tophat_out1h.bam" compare="sim_size" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
349 </test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
350 <!-- Test using base-space test data: paired-end reads, index from history. -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
351 <test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
352 <!-- TopHat commands:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
353 bowtie2-build -f test-data/tophat_in1.fasta tophat_in1
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
354 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
355 Rename the files in tmp_dir appropriately
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
356 -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
357 <param name="sPaired" value="paired" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
358 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
359 <param name="input2" ftype="fastqsanger" value="tophat_in3.fastqsanger" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
360 <param name="genomeSource" value="history" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
361 <param name="ownFile" ftype="fasta" value="tophat_in1.fasta" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
362 <param name="mate_inner_distance" value="20" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
363 <param name="settingsType" value="preSet" />
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
364 <param name="specReadGroup" value="no" />
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
365 <output name="junctions" file="tophat2_out2j.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
366 <output name="accepted_hits" file="tophat_out2h.bam" compare="sim_size" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
367 </test>
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
368 <test>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
369 <!-- Same test as above but with a collection. -->
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
370 <param name="sPaired" value="paired_collection" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
371 <param name="input">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
372 <collection type="paired">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
373 <element name="forward" value="tophat_in2.fastqsanger" ftype="fastqsanger" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
374 <element name="reverse" value="tophat_in3.fastqsanger" ftype="fastqsanger" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
375 </collection>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
376 </param>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
377 <param name="genomeSource" value="history" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
378 <param name="ownFile" ftype="fasta" value="tophat_in1.fasta" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
379 <param name="mate_inner_distance" value="20" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
380 <param name="settingsType" value="preSet" />
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
381 <param name="specReadGroup" value="no" />
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
382 <output name="junctions" file="tophat2_out2j.bed" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
383 <output name="accepted_hits" file="tophat_out2h.bam" compare="sim_size" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
384 </test>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
385 <!-- Test base-space single-end reads with user-supplied reference fasta and full parameters -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
386 <test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
387 <!-- Tophat commands:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
388 bowtie2-build -f test-data/tophat_in1.fasta tophat_in1
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
389 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
390 Replace the + with double-dash
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
391 Rename the files in tmp_dir appropriately
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
392 -->
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
393 <conditional name="singlePaired">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
394 <param name="sPaired" value="single"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
395 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
396 </conditional>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
397 <param name="genomeSource" value="history"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
398 <param name="ownFile" value="tophat_in1.fasta"/>
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
399 <conditional name="params">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
400 <param name="settingsType" value="full"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
401 <param name="library_type" value="FR Unstranded"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
402 <param name="read_mismatches" value="2"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
403 <param name="bowtie_n" value="No"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
404 <param name="anchor_length" value="8"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
405 <param name="splice_mismatches" value="0"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
406 <param name="min_intron_length" value="70"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
407 <param name="max_intron_length" value="500000"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
408 <param name="max_multihits" value="40"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
409 <param name="min_segment_intron" value="50" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
410 <param name="max_segment_intron" value="500000" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
411 <param name="seg_mismatches" value="2"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
412 <param name="seg_length" value="25"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
413 <conditional name="indel_search">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
414 <param name="allow_indel_search" value="Yes"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
415 <param name="max_insertion_length" value="3"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
416 <param name="max_deletion_length" value="3"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
417 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
418 <conditional name="own_junctions">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
419 <param name="use_junctions" value="Yes" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
420 <conditional name="gene_model_ann">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
421 <param name="use_annotations" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
422 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
423 <conditional name="raw_juncs">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
424 <param name="use_juncs" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
425 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
426 <conditional name="no_novel_juncs">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
427 <param name="no_novel_juncs" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
428 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
429 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
430 <conditional name="coverage_search">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
431 <param name="use_search" value="Yes" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
432 <param name="min_coverage_intron" value="50" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
433 <param name="max_coverage_intron" value="20000" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
434 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
435 <param name="microexon_search" value="Yes" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
436 <conditional name="bowtie2_settings">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
437 <param name="b2_settings" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
438 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
439 <!-- Fusion search params -->
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
440 <conditional name="fusion_search">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
441 <param name="do_search" value="Yes" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
442 <param name="anchor_len" value="21" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
443 <param name="min_dist" value="10000021" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
444 <param name="read_mismatches" value="3" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
445 <param name="multireads" value="4" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
446 <param name="multipairs" value="5" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
447 <param name="ignore_chromosomes" value="chrM"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
448 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
449 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
450 <conditional name="readGroup">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
451 <param name="specReadGroup" value="no" />
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
452 </conditional>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
453 <output name="insertions" file="tophat_out3i.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
454 <output name="deletions" file="tophat_out3d.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
455 <output name="junctions" file="tophat2_out3j.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
456 <output name="accepted_hits" file="tophat_out3h.bam" compare="sim_size" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
457 </test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
458 <!-- Test base-space paired-end reads with user-supplied reference fasta and full parameters -->
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
459 <test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
460 <!-- TopHat commands:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
461 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
462 Replace the + with double-dash
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
463 Rename the files in tmp_dir appropriately
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
464 -->
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
465 <conditional name="singlePaired">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
466 <param name="sPaired" value="paired"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
467 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
468 <param name="input2" ftype="fastqsanger" value="tophat_in3.fastqsanger"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
469 <param name="mate_inner_distance" value="20"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
470 <param name="report_discordant_pairs" value="Yes" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
471 </conditional>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
472 <param name="genomeSource" value="indexed"/>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
473 <param name="index" value="tophat_test"/>
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
474 <conditional name="params">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
475 <param name="settingsType" value="full"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
476 <param name="library_type" value="FR Unstranded"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
477 <param name="read_mismatches" value="5"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
478 <!-- Error: the read mismatches (5) and the read gap length (2) should be less than or equal to the read edit dist (2) -->
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
479 <param name="read_edit_dist" value="5" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
480 <param name="bowtie_n" value="Yes"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
481 <param name="mate_std_dev" value="20"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
482 <param name="anchor_length" value="8"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
483 <param name="splice_mismatches" value="0"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
484 <param name="min_intron_length" value="70"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
485 <param name="max_intron_length" value="500000"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
486 <param name="max_multihits" value="40"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
487 <param name="min_segment_intron" value="50" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
488 <param name="max_segment_intron" value="500000" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
489 <param name="seg_mismatches" value="2"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
490 <param name="seg_length" value="25"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
491 <conditional name="indel_search">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
492 <param name="allow_indel_search" value="No"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
493 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
494 <conditional name="own_junctions">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
495 <param name="use_junctions" value="Yes" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
496 <conditional name="gene_model_ann">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
497 <param name="use_annotations" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
498 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
499 <conditional name="raw_juncs">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
500 <param name="use_juncs" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
501 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
502 <conditional name="no_novel_juncs">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
503 <param name="no_novel_juncs" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
504 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
505 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
506 <conditional name="coverage_search">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
507 <param name="use_search" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
508 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
509 <param name="microexon_search" value="Yes" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
510 <conditional name="bowtie2_settings">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
511 <param name="b2_settings" value="No" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
512 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
513 <!-- Fusion search params -->
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
514 <conditional name="fusion_search">
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
515 <param name="do_search" value="Yes" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
516 <param name="anchor_len" value="21" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
517 <param name="min_dist" value="10000021" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
518 <param name="read_mismatches" value="3" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
519 <param name="multireads" value="4" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
520 <param name="multipairs" value="5" />
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
521 <param name="ignore_chromosomes" value="chrM"/>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
522 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
523 </conditional>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
524 <conditional name="readGroup">
5
5c5517d2a9e9 planemo upload commit a52cc16ed8d0d60e99742b55fccbdedcbb64b82c
devteam
parents: 3
diff changeset
525 <param name="specReadGroup" value="no" />
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
526 </conditional>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
527 <output name="junctions" file="tophat2_out4j.bed" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
528 <output name="accepted_hits" file="tophat_out4h.bam" compare="sim_size" />
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
529 </test>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
530 </tests>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
531 <help>
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
532 **Tophat Overview**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
533
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
534 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
535 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
536
1
ae06af1118dc Update Tophat URLs in help text.
Nate Coraor <nate@bx.psu.edu>
parents: 0
diff changeset
537 .. _Tophat: http://ccb.jhu.edu/software/tophat/
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
538
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
539 ------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
540
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
541 **Know what you are doing**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
542
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
543 .. class:: warningmark
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
544
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
545 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
546
1
ae06af1118dc Update Tophat URLs in help text.
Nate Coraor <nate@bx.psu.edu>
parents: 0
diff changeset
547 .. __: http://ccb.jhu.edu/software/tophat/manual.shtml
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
548
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
549 ------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
550
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
551 **Input formats**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
552
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
553 Tophat accepts files in Sanger FASTQ format. Use the FASTQ Groomer to prepare your files.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
554
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
555 ------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
556
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
557 **Outputs**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
558
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
559 Tophat produces two output files:
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
560
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
561 - 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
562 - accepted_hits -- A list of read alignments in BAM_ format.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
563
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
564 .. _BED: http://genome.ucsc.edu/FAQ/FAQformat.html#format1
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
565 .. _BAM: http://samtools.sourceforge.net/
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
566
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
567 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
568
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
569 -------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
570
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
571 **Tophat settings**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
572
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
573 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
574
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
575 ------
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
576
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
577 **Tophat parameter list**
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
578
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
579 This is a list of implemented Tophat options::
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
580
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
581 -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
582 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
583 is required for paired end runs.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
584 --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
585 -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
586 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
587 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
588 -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
589 -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
590 -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
591 -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
592 alignments. The default is 40.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
593 -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
594 -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
595 -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
596 --no-coverage-search Disables the coverage based search for junctions.
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
597 --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
598 --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
599 --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
600 --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
601 --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
602 --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
603 --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
604 --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
605 </help>
2
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
606 <citations>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
607 <citation type="doi">10.1186/gb-2013-14-4-r36</citation>
da1f39fe14bc Uploaded
devteam
parents: 1
diff changeset
608 </citations>
0
ffa30bedbee3 Imported from capsule None
devteam
parents:
diff changeset
609 </tool>