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