annotate tophat_wrapper.xml @ 1:af089ca8b4ee draft

planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
author devteam
date Tue, 13 Oct 2015 12:33:25 -0400
parents 51c6602b46b9
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
1 <tool id="tophat" name="TopHat for Illumina" version="1.5.0">
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
2 <!-- Wrapper compatible with TopHat versions 1.3.0 to 1.4.1 -->
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
3 <description>Find splice junctions using RNA-seq data</description>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
4 <version_command>tophat --version</version_command>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
5 <requirements>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
6 <requirement type="package" version="0.1.18">samtools</requirement>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
7 <requirement type="package" version="0.12.7">bowtie</requirement>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
8 <requirement type="package" version="1.4.0">tophat</requirement>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
9 </requirements>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
10 <command interpreter="python">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
11 tophat_wrapper.py
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
12 ## Change this to accommodate the number of threads you have available.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
13 --num-threads="\${GALAXY_SLOTS:-4}"
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
14 ## Provide outputs.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
15 --junctions-output=$junctions
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
16 --hits-output=$accepted_hits
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
17
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
18 ## Handle reference file.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
19 #if $refGenomeSource.genomeSource == "history":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
20 --own-file=$refGenomeSource.ownFile
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
21 #else:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
22 --indexes-path="${refGenomeSource.index.fields.path}"
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
23 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
24
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
25 ## Are reads single-end or paired?
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
26 --single-paired=$singlePaired.sPaired
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
27
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
28 ## First input file always required.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
29 --input1=$input1
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
30
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
31 ## Set params based on whether reads are single-end or paired.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
32 #if $singlePaired.sPaired == "single":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
33 --settings=$singlePaired.sParams.sSettingsType
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
34 #if $singlePaired.sParams.sSettingsType == "full":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
35 -a $singlePaired.sParams.anchor_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
36 -m $singlePaired.sParams.splice_mismatches
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
37 -i $singlePaired.sParams.min_intron_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
38 -I $singlePaired.sParams.max_intron_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
39 -g $singlePaired.sParams.max_multihits
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
40 --min-segment-intron $singlePaired.sParams.min_segment_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
41 --max-segment-intron $singlePaired.sParams.max_segment_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
42 --initial-read-mismatches=$singlePaired.sParams.initial_read_mismatches
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
43 --seg-mismatches=$singlePaired.sParams.seg_mismatches
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
44 --seg-length=$singlePaired.sParams.seg_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
45 --library-type=$singlePaired.sParams.library_type
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
46
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
47 ## Indel search.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
48 #if $singlePaired.sParams.indel_search.allow_indel_search == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
49 ## --allow-indels
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
50 --max-insertion-length $singlePaired.sParams.indel_search.max_insertion_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
51 --max-deletion-length $singlePaired.sParams.indel_search.max_deletion_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
52 #else:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
53 --no-novel-indels
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
54 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
55
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
56 ## Supplying junctions parameters.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
57 #if $singlePaired.sParams.own_junctions.use_junctions == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
58 #if $singlePaired.sParams.own_junctions.gene_model_ann.use_annotations == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
59 -G $singlePaired.sParams.own_junctions.gene_model_ann.gene_annotation_model
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
60 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
61 #if $singlePaired.sParams.own_junctions.raw_juncs.use_juncs == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
62 -j $singlePaired.sParams.own_junctions.raw_juncs.raw_juncs
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
63 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
64 ## TODO: No idea why a string cast is necessary, but it is:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
65 #if str($singlePaired.sParams.own_junctions.no_novel_juncs) == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
66 --no-novel-juncs
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
67 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
68 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
69
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
70 #if $singlePaired.sParams.closure_search.use_search == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
71 --closure-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
72 --min-closure-exon $singlePaired.sParams.closure_search.min_closure_exon
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
73 --min-closure-intron $singlePaired.sParams.closure_search.min_closure_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
74 --max-closure-intron $singlePaired.sParams.closure_search.max_closure_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
75 #else:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
76 --no-closure-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
77 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
78 #if $singlePaired.sParams.coverage_search.use_search == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
79 --coverage-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
80 --min-coverage-intron $singlePaired.sParams.coverage_search.min_coverage_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
81 --max-coverage-intron $singlePaired.sParams.coverage_search.max_coverage_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
82 #else:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
83 --no-coverage-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
84 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
85 ## TODO: No idea why the type conversion is necessary, but it seems to be.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
86 #if str($singlePaired.sParams.microexon_search) == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
87 --microexon-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
88 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
89 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
90 #else:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
91 --input2=$singlePaired.input2
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
92 -r $singlePaired.mate_inner_distance
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
93 --settings=$singlePaired.pParams.pSettingsType
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
94 #if $singlePaired.pParams.pSettingsType == "full":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
95 --mate-std-dev=$singlePaired.pParams.mate_std_dev
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
96 -a $singlePaired.pParams.anchor_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
97 -m $singlePaired.pParams.splice_mismatches
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
98 -i $singlePaired.pParams.min_intron_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
99 -I $singlePaired.pParams.max_intron_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
100 -g $singlePaired.pParams.max_multihits
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
101 --min-segment-intron $singlePaired.pParams.min_segment_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
102 --max-segment-intron $singlePaired.pParams.max_segment_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
103 --initial-read-mismatches=$singlePaired.pParams.initial_read_mismatches
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
104 --seg-mismatches=$singlePaired.pParams.seg_mismatches
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
105 --seg-length=$singlePaired.pParams.seg_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
106 --library-type=$singlePaired.pParams.library_type
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
107
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
108 ## Indel search.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
109 #if $singlePaired.pParams.indel_search.allow_indel_search == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
110 ## --allow-indels
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
111 --max-insertion-length $singlePaired.pParams.indel_search.max_insertion_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
112 --max-deletion-length $singlePaired.pParams.indel_search.max_deletion_length
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
113 #else:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
114 --no-novel-indels
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
115 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
116
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
117 ## Supplying junctions parameters.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
118 #if $singlePaired.pParams.own_junctions.use_junctions == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
119 #if $singlePaired.pParams.own_junctions.gene_model_ann.use_annotations == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
120 -G $singlePaired.pParams.own_junctions.gene_model_ann.gene_annotation_model
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
121 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
122 #if $singlePaired.pParams.own_junctions.raw_juncs.use_juncs == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
123 -j $singlePaired.pParams.own_junctions.raw_juncs.raw_juncs
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
124 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
125 ## TODO: No idea why type cast is necessary, but it is:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
126 #if str($singlePaired.pParams.own_junctions.no_novel_juncs) == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
127 --no-novel-juncs
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
128 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
129 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
130
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
131 #if $singlePaired.pParams.closure_search.use_search == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
132 --closure-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
133 --min-closure-exon $singlePaired.pParams.closure_search.min_closure_exon
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
134 --min-closure-intron $singlePaired.pParams.closure_search.min_closure_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
135 --max-closure-intron $singlePaired.pParams.closure_search.max_closure_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
136 #else:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
137 --no-closure-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
138 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
139 #if $singlePaired.pParams.coverage_search.use_search == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
140 --coverage-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
141 --min-coverage-intron $singlePaired.pParams.coverage_search.min_coverage_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
142 --max-coverage-intron $singlePaired.pParams.coverage_search.max_coverage_intron
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
143 #else:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
144 --no-coverage-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
145 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
146 ## TODO: No idea why the type conversion is necessary, but it seems to be.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
147 #if str ($singlePaired.pParams.microexon_search) == "Yes":
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
148 --microexon-search
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
149 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
150 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
151 #end if
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
152 </command>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
153 <inputs>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
154 <param format="fastqsanger" name="input1" type="data" label="RNA-Seq FASTQ file" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
155 <expand macro="refGenomeSourceConditional">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
156 <options from_data_table="tophat_indexes">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
157 <filter type="sort_by" column="2"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
158 <validator type="no_options" message="No genomes are available for the selected input dataset"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
159 </options>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
160 </expand>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
161 <conditional name="singlePaired">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
162 <param name="sPaired" type="select" label="Is this library mate-paired?">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
163 <option value="single">Single-end</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
164 <option value="paired">Paired-end</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
165 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
166 <when value="single">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
167 <conditional name="sParams">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
168 <param name="sSettingsType" type="select" label="TopHat settings to use" help="Use the Full parameter list to change default settings.">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
169 <option value="preSet">Default settings</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
170 <option value="full">Full parameter list</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
171 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
172 <when value="preSet" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
173 <!-- Full/advanced params. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
174 <when value="full">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
175 <param name="library_type" type="select" label="Library Type" help="TopHat will treat the reads as strand specific. Every read alignment will have an XS attribute tag. Consider supplying library type options below to select the correct RNA-seq protocol.">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
176 <option value="fr-unstranded">FR Unstranded</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
177 <option value="fr-firststrand">FR First Strand</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
178 <option value="fr-secondstrand">FR Second Strand</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
179 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
180 <param name="anchor_length" type="integer" value="8" label="Anchor length (at least 3)" help="Report junctions spanned by reads with at least this many bases on each side of the junction." />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
181 <param name="splice_mismatches" type="integer" value="0" label="Maximum number of mismatches that can appear in the anchor region of spliced alignment" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
182 <param name="min_intron_length" type="integer" value="70" label="The minimum intron length" help="TopHat will ignore donor/acceptor pairs closer than this many bases apart." />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
183 <param name="max_intron_length" type="integer" value="500000" label="The maximum intron length" help="When searching for junctions ab initio, TopHat will ignore donor/acceptor pairs farther than this many bases apart, except when such a pair is supported by a split segment alignment of a long read." />
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
184 <expand macro="indel_searchConditional" />
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
185 <param name="max_multihits" type="integer" value="20" label="Maximum number of alignments to be allowed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
186 <param name="min_segment_intron" type="integer" value="50" label="Minimum intron length that may be found during split-segment (default) search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
187 <param name="max_segment_intron" type="integer" value="500000" label="Maximum intron length that may be found during split-segment (default) search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
188 <param name="initial_read_mismatches" type="integer" min="0" value="2" label="Number of mismatches allowed in the initial read mapping" />
51c6602b46b9 Imported from capsule None
devteam
parents:
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" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
190 <param name="seg_length" type="integer" value="25" label="Minimum length of read segments" />
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
191
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
192 <!-- Options for supplying own junctions. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
193 <conditional name="own_junctions">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
194 <param name="use_junctions" type="select" label="Use Own Junctions">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
195 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
196 <option value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
197 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
198 <when value="Yes">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
199 <conditional name="gene_model_ann">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
200 <param name="use_annotations" type="select" label="Use Gene Annotation Model">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
201 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
202 <option value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
203 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
204 <when value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
205 <when value="Yes">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
206 <param format="gtf,gff3" name="gene_annotation_model" type="data" label="Gene Model Annotations" help="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."/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
207 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
208 </conditional>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
209 <conditional name="raw_juncs">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
210 <param name="use_juncs" type="select" label="Use Raw Junctions">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
211 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
212 <option value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
213 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
214 <when value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
215 <when value="Yes">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
216 <param format="interval" name="raw_juncs" type="data" label="Raw Junctions" help="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."/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
217 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
218 </conditional>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
219 <param name="no_novel_juncs" type="select" label="Only look for supplied junctions">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
220 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
221 <option value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
222 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
223 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
224 <when value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
225 </conditional> <!-- /own_junctions -->
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
226
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
227 <!-- Closure search. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
228 <conditional name="closure_search">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
229 <param name="use_search" type="select" label="Use Closure Search">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
230 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
231 <option value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
232 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
233 <when value="Yes">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
234 <param name="min_closure_exon" type="integer" value="50" label="During closure search for paired end reads, exonic hops in the potential splice graph must be at least this long. The default is 50." />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
235 <param name="min_closure_intron" type="integer" value="50" label="Minimum intron length that may be found during closure search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
236 <param name="max_closure_intron" type="integer" value="5000" label="Maximum intron length that may be found during closure search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
237 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
238 <when value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
239 </conditional>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
240 <!-- Coverage search. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
241 <conditional name="coverage_search">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
242 <param name="use_search" type="select" label="Use Coverage Search">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
243 <option selected="true" value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
244 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
245 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
246 <when value="Yes">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
247 <param name="min_coverage_intron" type="integer" value="50" label="Minimum intron length that may be found during coverage search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
248 <param name="max_coverage_intron" type="integer" value="20000" label="Maximum intron length that may be found during coverage search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
249 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
250 <when value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
251 </conditional>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
252 <param name="microexon_search" type="select" label="Use Microexon Search" help="With this option, the pipeline will attempt to find alignments incident to microexons. Works only for reads 50bp or longer.">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
253 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
254 <option value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
255 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
256 </when> <!-- full -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
257 </conditional> <!-- sParams -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
258 </when> <!-- single -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
259 <when value="paired">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
260 <param format="fastqsanger" name="input2" type="data" label="RNA-Seq FASTQ file" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
261 <param name="mate_inner_distance" type="integer" value="20" label="Mean Inner Distance between Mate Pairs" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
262 <conditional name="pParams">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
263 <param name="pSettingsType" type="select" label="TopHat settings to use" help="Use the Full parameter list to change default settings.">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
264 <option value="preSet">Default settings</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
265 <option value="full">Full parameter list</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
266 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
267 <when value="preSet" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
268 <!-- Full/advanced params. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
269 <when value="full">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
270 <param name="library_type" type="select" label="Library Type" help="TopHat will treat the reads as strand specific. Every read alignment will have an XS attribute tag. Consider supplying library type options below to select the correct RNA-seq protocol.">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
271 <option value="fr-unstranded">FR Unstranded</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
272 <option value="fr-firststrand">FR First Strand</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
273 <option value="fr-secondstrand">FR Second Strand</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
274 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
275 <param name="mate_std_dev" type="integer" value="20" label="Std. Dev for Distance between Mate Pairs" help="The standard deviation for the distribution on inner distances between mate pairs."/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
276 <param name="anchor_length" type="integer" value="8" label="Anchor length (at least 3)" help="Report junctions spanned by reads with at least this many bases on each side of the junction." />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
277 <param name="splice_mismatches" type="integer" value="0" label="Maximum number of mismatches that can appear in the anchor region of spliced alignment" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
278 <param name="min_intron_length" type="integer" value="70" label="The minimum intron length" help="TopHat will ignore donor/acceptor pairs closer than this many bases apart." />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
279 <param name="max_intron_length" type="integer" value="500000" label="The maximum intron length" help="When searching for junctions ab initio, TopHat will ignore donor/acceptor pairs farther than this many bases apart, except when such a pair is supported by a split segment alignment of a long read." />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
280 <expand macro="indel_searchConditional" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
281 <param name="max_multihits" type="integer" value="20" label="Maximum number of alignments to be allowed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
282 <param name="min_segment_intron" type="integer" value="50" label="Minimum intron length that may be found during split-segment (default) search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
283 <param name="max_segment_intron" type="integer" value="500000" label="Maximum intron length that may be found during split-segment (default) search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
284 <param name="initial_read_mismatches" type="integer" min="0" value="2" label="Number of mismatches allowed in the initial read mapping" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
285 <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" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
286 <param name="seg_length" type="integer" value="25" label="Minimum length of read segments" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
287 <!-- Options for supplying own junctions. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
288 <expand macro="own_junctionsConditional" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
289 <!-- Closure search. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
290 <conditional name="closure_search">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
291 <param name="use_search" type="select" label="Use Closure Search">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
292 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
293 <option value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
294 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
295 <when value="Yes">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
296 <param name="min_closure_exon" type="integer" value="50" label="During closure search for paired end reads, exonic hops in the potential splice graph must be at least this long. The default is 50." />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
297 <param name="min_closure_intron" type="integer" value="50" label="Minimum intron length that may be found during closure search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
298 <param name="max_closure_intron" type="integer" value="5000" label="Maximum intron length that may be found during closure search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
299 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
300 <when value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
301 </conditional>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
302 <!-- Coverage search. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
303 <conditional name="coverage_search">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
304 <param name="use_search" type="select" label="Use Coverage Search">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
305 <option selected="true" value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
306 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
307 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
308 <when value="Yes">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
309 <param name="min_coverage_intron" type="integer" value="50" label="Minimum intron length that may be found during coverage search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
310 <param name="max_coverage_intron" type="integer" value="20000" label="Maximum intron length that may be found during coverage search" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
311 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
312 <when value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
313 </conditional>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
314 <param name="microexon_search" type="select" label="Use Microexon Search" help="With this option, the pipeline will attempt to find alignments incident to microexons. Works only for reads 50bp or longer.">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
315 <option value="No">No</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
316 <option value="Yes">Yes</option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
317 </param>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
318 </when> <!-- full -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
319 </conditional> <!-- pParams -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
320 </when> <!-- paired -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
321 </conditional>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
322 </inputs>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
323
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
324 <outputs>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
325 <data format="bed" name="insertions" label="${tool.name} on ${on_string}: insertions" from_work_dir="tophat_out/insertions.bed">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
326 <expand macro="dbKeyActions" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
327 </data>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
328 <data format="bed" name="deletions" label="${tool.name} on ${on_string}: deletions" from_work_dir="tophat_out/deletions.bed">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
329 <expand macro="dbKeyActions" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
330 </data>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
331 <data format="bed" name="junctions" label="${tool.name} on ${on_string}: splice junctions" from_work_dir="tophat_out/junctions.bed">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
332 <expand macro="dbKeyActions" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
333 </data>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
334 <data format="bam" name="accepted_hits" label="${tool.name} on ${on_string}: accepted_hits" from_work_dir="tophat_out/accepted_hits.bam">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
335 <expand macro="dbKeyActions" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
336 </data>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
337 </outputs>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
338 <macros>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
339 <import>tophat_macros.xml</import>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
340 <macro name="dbKeyActions">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
341 <actions>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
342 <conditional name="refGenomeSource.genomeSource">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
343 <when value="indexed">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
344 <action type="metadata" name="dbkey">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
345 <option type="from_data_table" name="tophat_indexes" column="1" offset="0">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
346 <filter type="param_value" column="0" value="#" compare="startswith" keep="False"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
347 <filter type="param_value" ref="refGenomeSource.index" column="0"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
348 </option>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
349 </action>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
350 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
351 <when value="history">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
352 <action type="metadata" name="dbkey">
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
353 <option type="from_param" name="refGenomeSource.ownFile" param_attribute="dbkey" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
354 </action>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
355 </when>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
356 </conditional>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
357 </actions>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
358 </macro>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
359 </macros>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
360 <tests>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
361 <!-- Test base-space single-end reads with pre-built index and preset parameters -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
362 <test>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
363 <!-- TopHat commands:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
364 tophat -o tmp_dir -p 1 tophat_in1 test-data/tophat_in2.fastqsanger
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
365 Rename the files in tmp_dir appropriately
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
366 -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
367 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
368 <param name="genomeSource" value="indexed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
369 <param name="index" value="tophat_test" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
370 <param name="sPaired" value="single" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
371 <param name="sSettingsType" value="preSet" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
372 <output name="junctions" file="tophat_out1j.bed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
373 <output name="accepted_hits" file="tophat_out1h.bam" compare="sim_size" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
374 </test>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
375 <!-- Test using base-space test data: paired-end reads, index from history. -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
376 <test>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
377 <!-- TopHat commands:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
378 bowtie-build -f test-data/tophat_in1.fasta tophat_in1
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
379 tophat -o tmp_dir -p 1 -r 20 tophat_in1 test-data/tophat_in2.fastqsanger test-data/tophat_in3.fastqsanger
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
380 Rename the files in tmp_dir appropriately
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
381 -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
382 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
383 <param name="genomeSource" value="history" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
384 <param name="ownFile" ftype="fasta" value="tophat_in1.fasta" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
385 <param name="sPaired" value="paired" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
386 <param name="input2" ftype="fastqsanger" value="tophat_in3.fastqsanger" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
387 <param name="mate_inner_distance" value="20" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
388 <param name="pSettingsType" value="preSet" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
389 <output name="junctions" file="tophat_out2j.bed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
390 <output name="accepted_hits" file="tophat_out2h.bam" compare="sim_size" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
391 </test>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
392 <!-- Test base-space single-end reads with user-supplied reference fasta and full parameters -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
393 <test>
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
394 <!-- TopHat commands:
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
395 bowtie-build -f test-data/tophat_in1.fasta tophat_in1
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
396 tophat -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 +closure-search +min-closure-exon 50 +min-closure-intron 50 +max-closure-intro 5000 +microexon-search tophat_in1 test-data/tophat_in2.fastqsanger
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
397 Replace the + with double-dash
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
398 Rename the files in tmp_dir appropriately
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
399 -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
400 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
401 <param name="genomeSource" value="history"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
402 <param name="ownFile" value="tophat_in1.fasta"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
403 <param name="sPaired" value="single"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
404 <param name="sSettingsType" value="full"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
405 <param name="library_type" value="FR Unstranded"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
406 <param name="anchor_length" value="8"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
407 <param name="splice_mismatches" value="0"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
408 <param name="min_intron_length" value="70"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
409 <param name="max_intron_length" value="500000"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
410 <param name="max_multihits" value="40"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
411 <param name="min_segment_intron" value="50" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
412 <param name="max_segment_intron" value="500000" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
413 <param name="seg_mismatches" value="2"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
414 <param name="seg_length" value="25"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
415 <param name="allow_indel_search" value="Yes"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
416 <param name="max_insertion_length" value="3"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
417 <param name="max_deletion_length" value="3"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
418 <param name="use_junctions" value="Yes" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
419 <param name="use_annotations" value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
420 <param name="use_juncs" value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
421 <param name="no_novel_juncs" value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
422 <param name="use_search" value="Yes" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
423 <param name="min_closure_exon" value="50" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
424 <param name="min_closure_intron" value="50" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
425 <param name="max_closure_intron" value="5000" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
426 <param name="use_search" value="Yes" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
427 <param name="min_coverage_intron" value="50" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
428 <param name="max_coverage_intron" value="20000" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
429 <param name="microexon_search" value="Yes" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
430 <output name="insertions" file="tophat_out3i.bed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
431 <output name="deletions" file="tophat_out3d.bed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
432 <output name="junctions" file="tophat_out3j.bed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
433 <output name="accepted_hits" file="tophat_out3h.bam" compare="sim_size" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
434 </test>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
435 <!-- Test base-space paired-end reads with user-supplied reference fasta and full parameters -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
436 <test>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
437 <!-- TopHat commands:
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
438 tophat -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 tophat_in1 test-data/tophat_in2.fastqsanger test-data/tophat_in3.fastqsanger
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
439 Replace the + with double-dash
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
440 Rename the files in tmp_dir appropriately
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
441 -->
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
442 <param name="input1" ftype="fastqsanger" value="tophat_in2.fastqsanger"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
443 <param name="genomeSource" value="indexed"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
444 <param name="index" value="tophat_test"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
445 <param name="sPaired" value="paired"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
446 <param name="input2" ftype="fastqsanger" value="tophat_in3.fastqsanger"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
447 <param name="mate_inner_distance" value="20"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
448 <param name="pSettingsType" value="full"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
449 <param name="library_type" value="FR Unstranded"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
450 <param name="mate_std_dev" value="20"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
451 <param name="anchor_length" value="8"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
452 <param name="splice_mismatches" value="0"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
453 <param name="min_intron_length" value="70"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
454 <param name="max_intron_length" value="500000"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
455 <param name="max_multihits" value="40"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
456 <param name="min_segment_intron" value="50" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
457 <param name="max_segment_intron" value="500000" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
458 <param name="seg_mismatches" value="2"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
459 <param name="seg_length" value="25"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
460 <param name="allow_indel_search" value="No"/>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
461 <param name="use_junctions" value="Yes" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
462 <param name="use_annotations" value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
463 <param name="use_juncs" value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
464 <param name="no_novel_juncs" value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
465 <param name="use_search" value="No" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
466 <param name="microexon_search" value="Yes" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
467 <output name="junctions" file="tophat_out4j.bed" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
468 <output name="accepted_hits" file="tophat_out4h.bam" compare="sim_size" />
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
469 </test>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
470 </tests>
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
471
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
472 <help>
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
473 **TopHat Overview**
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
474
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
475 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, and then analyzes the mapping results to identify splice junctions between exons.
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
476
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
477 .. _TopHat: http://ccb.jhu.edu/software/tophat/
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
478
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
479 ------
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
480
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
481 **Know what you are doing**
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
482
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
483 .. class:: warningmark
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
484
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
485 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.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
486
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
487 .. __: http://ccb.jhu.edu/software/tophat/manual.shtml
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
488
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
489 ------
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
490
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
491 **Input formats**
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
492
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
493 TopHat accepts files in Sanger FASTQ format. Use the FASTQ Groomer to prepare your files.
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
494
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
495 ------
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
496
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
497 **Outputs**
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
498
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
499 TopHat produces two output files:
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
500
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
501 - 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.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
502 - accepted_hits -- A list of read alignments in BAM_ format.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
503
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
504 .. _BED: http://genome.ucsc.edu/FAQ/FAQformat.html#format1
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
505 .. _BAM: http://samtools.sourceforge.net/
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
506
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
507 Two other possible outputs, depending on the options you choose, are insertions and deletions, both of which are in BED format.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
508
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
509 -------
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
510
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
511 **TopHat settings**
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
512
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
513 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
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
514
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
515 ------
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
516
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
517 **TopHat parameter list**
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
518
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
519 This is a list of implemented TopHat options::
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
520
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
521 -r This is the expected (mean) inner distance between mate pairs. For, example, for paired end runs with fragments
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
522 selected at 300bp, where each end is 50bp, you should set -r to be 200. There is no default, and this parameter
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
523 is required for paired end runs.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
524 --mate-std-dev INT The standard deviation for the distribution on inner distances between mate pairs. The default is 20bp.
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
525 -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
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
526 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
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
527 read with this many bases on each side. This must be at least 3 and the default is 8.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
528 -m/--splice-mismatches INT The maximum number of mismatches that may appear in the "anchor" region of a spliced alignment. The default is 0.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
529 -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.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
530 -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.
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
531 -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
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
532 alignments. The default is 40.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
533 -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.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
534 -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.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
535 -no-novel-juncs Only look for junctions indicated in the supplied GFF file. (ignored without -G)
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
536 --no-closure-search Disables the mate pair closure-based search for junctions. Currently, has no effect - closure search is off by default.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
537 --closure-search Enables the mate pair closure-based search for junctions. Closure-based search should only be used when the expected inner distance between mates is small (about or less than 50bp)
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
538 --no-coverage-search Disables the coverage based search for junctions.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
539 --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.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
540 --microexon-search With this option, the pipeline will attempt to find alignments incident to microexons. Works only for reads 50bp or longer.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
541 --butterfly-search TopHat will use a slower but potentially more sensitive algorithm to find junctions in addition to its standard search. Consider using this if you expect that your experiment produced a lot of reads from pre-mRNA, that fall within the introns of your transcripts.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
542 --segment-mismatches Read segments are mapped independently, allowing up to this many mismatches in each segment alignment. The default is 2.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
543 --segment-length Each read is cut up into segments, each at least this long. These segments are mapped independently. The default is 25.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
544 --min-closure-exon During closure search for paired end reads, exonic hops in the potential splice graph must be at least this long. The default is 50.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
545 --min-closure-intron The minimum intron length that may be found during closure search. The default is 50.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
546 --max-closure-intron The maximum intron length that may be found during closure search. The default is 5000.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
547 --min-coverage-intron The minimum intron length that may be found during coverage search. The default is 50.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
548 --max-coverage-intron The maximum intron length that may be found during coverage search. The default is 20000.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
549 --min-segment-intron The minimum intron length that may be found during split-segment search. The default is 50.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
550 --max-segment-intron The maximum intron length that may be found during split-segment search. The default is 500000.
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
551 </help>
1
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
552 <citations>
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
553 <citation type="doi">10.1093/bioinformatics/btp120</citation>
af089ca8b4ee planemo upload commit 33927a87ba2eee9bf0ecdd376a66241b17b3d734
devteam
parents: 0
diff changeset
554 </citations>
0
51c6602b46b9 Imported from capsule None
devteam
parents:
diff changeset
555 </tool>