Mercurial > repos > rnateam > segemehl

comparison segemehl.xml @ 11:d82af79d5f27 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/segemehl commit 3accce77966950e310dbfa7576d5b3c828b9c60c
author bgruening
date Thu, 23 Aug 2018 13:27:21 -0400
parents 90b8941e5115
children
comparison
equal deleted inserted replaced
10:ee6a5ec6963c 11:d82af79d5f27
1 <tool id="segemehl" name="segemehl" version="0.2.0.3"> 1 <tool id="segemehl" name="segemehl" version="0.2.0.4">
2 <description>short read mapping with gaps</description> 2 <description>short read mapping with gaps</description>
3 <requirements> 3 <requirements>
4 <requirement type="package" version="0.2.0">segemehl</requirement> 4 <requirement type="package" version="0.2.0">segemehl</requirement>
5 </requirements> 5 </requirements>
6 <stdio> 6 <stdio>
8 source="both" 8 source="both"
9 level="fatal" 9 level="fatal"
10 description="Execution halted." /> 10 description="Execution halted." />
11 </stdio> 11 </stdio>
12 <command> 12 <command>
13 <!-- 13 <!--
14 ## check for single/pair-end 14 ## check for single/pair-end
15 #if str( $library.type ) == "single": 15 #if str( $library.type ) == "single":
16 #set $query_list = list() 16 #set $query_list = list()
17 ## prepare inputs 17 ## prepare inputs
18 #for $fastq in $library.input_query: 18 #for $fastq in $library.input_query:
31 -q #echo ','.join($mate1) 31 -q #echo ','.join($mate1)
32 -p #echo ','.join($mate2) 32 -p #echo ','.join($mate2)
33 33
34 -I $library.maxinsertsize 34 -I $library.maxinsertsize
35 #end if 35 #end if
36 --> 36 -->
37 <![CDATA[ 37 <![CDATA[
38 ## prepare segemehl index if no reference genome is supplied 38 ## UNIMPLEMENTED
39 ## [SEEDEXTENSIONPARAMS]
40 ## -e, --extensionscore <n> score of a match during extension (default:2)
41 ## -n, --extensionpenalty <n> penalty for a mismatch during extension (default:4)
42 ## -X, --dropoff <n> dropoff parameter for extension (default:8)
43 ## --showalign
44
45 ## prepare segemehl index if no reference genome is supplied
39 #if $refGenomeSource.genomeSource == "history": 46 #if $refGenomeSource.genomeSource == "history":
40 mkdir ./temp_index/ && 47 mkdir ./temp_index/ &&
41 #set $temp_index = './temp_index/temp.idx' 48 #set $temp_index = './temp_index/temp.idx'
42 segemehl.x -x $temp_index -d $refGenomeSource.own_reference_genome && 49 segemehl.x -x $temp_index -d $refGenomeSource.own_reference_genome &&
43 #else: 50 #else:
58 65
59 -i $temp_index 66 -i $temp_index
60 67
61 ## check for single/pair-end 68 ## check for single/pair-end
62 #if str( $library.type ) == "single": 69 #if str( $library.type ) == "single":
63 ## prepare inputs 70 ## prepare inputs
64 -q ${library.input_query} 71 -q ${library.input_query}
65 #else 72 #else
66 -q ${mate_pair.first_strand_query} 73 -q ${mate_pair.first_strand_query}
67 -p ${mate_pair.second_strand_query} 74 -p ${mate_pair.second_strand_query}
68
69 -I ${library.maxinsertsize} 75 -I ${library.maxinsertsize}
76 #end if
77 #if str( $bisulfite ) != "0":
78 -F $bisulfite
70 #end if 79 #end if
71 -m $minsize 80 -m $minsize
72 -A $accuracy 81 -A $accuracy
73 -H $hitstrategy 82 -H $hitstrategy
74 #if str( $prime5 ).strip(): 83 #if str( $prime5 ).strip():
75 -P "$prime5" 84 -P "$prime5"
76 #end if 85 #end if
77 #if str( $prime3 ).strip(): 86 #if str( $prime3 ).strip():
78 -Q "$prime3" 87 -Q "$prime3"
79 #end if 88 #end if
80 $polyA 89 -R $clipacc
81 $autoclip 90 $polyA
82 $hardclip 91 $autoclip
83 $order 92 $hardclip
93 $order
84 #if $maxout: 94 #if $maxout:
85 --maxout $maxout 95 --maxout $maxout
86 #end if 96 #end if
87 #if str( $splitreads.splits ) == "splits": 97 #if str( $splitreads.splits ) == "splits":
88 --splits 98 --splits
89 --minsplicecover $splitreads.minsplicecover 99 --minsplicecover $splitreads.minsplicecover
90 --minfragscore $splitreads.minfragscore 100 --minfragscore $splitreads.minfragscore
91 --minfraglen $splitreads.minfraglen 101 --minfraglen $splitreads.minfraglen
92 --splicescorescale $splitreads.splicescorescale 102 --splicescorescale $splitreads.splicescorescale
103 --maxsplitevalue $splitreads.maxsplitevalue
93 #end if 104 #end if
94 -M $maxinterval 105 -M $maxinterval
95 -E $evalue 106 -E $evalue
96 -D $differences 107 -D $differences
108 -J $jump
97 -s 109 -s
98 -o '$segemehl_out' 110 -o '$segemehl_out'
111 #if str( $nomatchfilename ) == 'yes':
112 -u '$segemehl_outunmatched'
113 #end if
99 ]]> 114 ]]>
100 </command> 115 </command>
101 <inputs> 116 <inputs>
102 <conditional name="refGenomeSource"> 117 <conditional name="refGenomeSource">
103 <param name="genomeSource" type="select" label="Will you select a reference genome from your history or use a built-in index?" help="Built-ins were indexed using default options"> 118 <param name="genomeSource" type="select" label="Will you select a reference genome from your history or use a built-in index?" help="Built-ins were indexed using default options">
135 <!-- ToDo paired coolections --> 150 <!-- ToDo paired coolections -->
136 <repeat name="mate_list" title="Paired End Pairs" min="1"> 151 <repeat name="mate_list" title="Paired End Pairs" min="1">
137 <param name="first_strand_query" type="data" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads from first strand" /> 152 <param name="first_strand_query" type="data" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads from first strand" />
138 <param name="second_strand_query" type="data" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads from second strand" /> 153 <param name="second_strand_query" type="data" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads from second strand" />
139 </repeat> 154 </repeat>
140 <param name="maxinsertsize" type="integer" value="5000" label="Maximum size of the inserts (paired end)" help="default: 5000 (-I)" /> 155 <param argument="--maxinsertsize" type="integer" value="5000" label="Maximum size of the inserts (paired end)" help="default: 5000" />
141 </when> 156 </when>
142 </conditional> 157 </conditional>
143 158
144 <conditional name="splitreads"> 159 <conditional name="splitreads">
145 <param name="splits" type="select" label="Detect split/spliced reads" help="(--splits)"> 160 <param argument="splits" type="select" label="Detect split/spliced reads">
146 <option value="nosplit">No splits</option> 161 <option value="nosplit">No splits</option>
147 <option value="splits">Split reads</option> 162 <option value="splits">Split reads</option>
148 </param> 163 </param>
149 <when value="splits"> 164 <when value="splits">
150 <param name="minsplicecover" type="integer" value="80" label="Min coverage for spliced transcripts" help="(--minsplicecover)" /> 165 <param argument="--minsplicecover" type="integer" value="80" label="Min coverage for spliced transcripts" />
151 <param name="minfragscore" type="integer" value="18" label="Min coverage for spliced transcripts" help="(--minfragscore)" /> 166 <param argument="--minfragscore" type="integer" value="18" label="Min coverage for spliced transcripts" />
152 <param name="minfraglen" type="integer" value="20" label="Min length of a spliced fragment" help="(--minfraglen)" /> 167 <param argument="--minfraglen" type="integer" value="20" label="Min length of a spliced fragment" />
153 <param name="splicescorescale" type="float" value="1.0" label="Report spliced alignment with score greater than this scale times the score" 168 <param argument="--splicescorescale" type="float" value="1.0" label="Report spliced alignment with score greater than this scale times the score"
154 help="Report only if this value x score is larger than next best spliced alignment (--splicescorescale)" /> 169 help="Report only if this value x score is larger than next best spliced alignment" />
155 <param name="sevalue" type="float" min="0" value="50.000000" label="max split evalue" help="(--maxsplitevalue)"/> 170 <param argument="--maxsplitevalue" type="float" min="0" value="50.000000" label="max evalue for splits"/>
156 </when> 171 </when>
157 <when value="nosplit"> 172 <when value="nosplit">
158 </when> 173 </when>
159 </conditional> 174 </conditional>
160 175 <param argument="--bisulfite" type="select" label="Bisulfite mapping">
161 <param name="minsize" type="integer" value="12" min="1" label="Minimum size of queries" help="(-m)" /> 176 <option value="0">No bisulfite mapping</option>
162 <param name="maxout" type="integer" min="0" value="0" optional="True" 177 <option value="1">bisulfite mapping with methylC-seq/Lister et al.</option>
163 label="Maximum number of alignments that will be reported" help="(--maxout)" /> 178 <option value="2">bs-seq/Cokus et al. protocol</option>
164 <param name="accuracy" type="integer" value="85" min="1" max="100" label="Min percentage of matches per read in semi-global alignment" help="(-A)" /> 179 </param>
165 <param name="hitstrategy" type="select" label="Hits to report?" help="(-H)"> 180 <param argument="--minsize" type="integer" value="12" min="1" label="Minimum size of queries" />
181 <param argument="--maxout" type="integer" min="0" value="0" optional="True"
182 label="Maximum number of alignments that will be reported"/>
183 <param argument="--accuracy" type="integer" value="85" min="1" max="100" label="Min percentage of matches per read in semi-global alignment" />
184 <param argument="--hitstrategy" type="select" label="Hits to report?">
166 <option value="1">report only best scoring hits</option> 185 <option value="1">report only best scoring hits</option>
167 <option value="0">report all scoring hits</option> 186 <option value="0">report all scoring hits</option>
168 </param> 187 </param>
169 <param name="prime5" type="text" label="add 5' adapter" help="default: none (-Q)" /> 188 <param argument="--prime5" type="text" label="add 5' adapter" help="default: none" />
170 <param name="prime3" type="text" label="add 3' adapter" help="default: none (-P)"/> 189 <param argument="--prime3" type="text" label="add 3' adapter" help="default: none"/>
171 <param name="polyA" type="boolean" truevalue="--polyA" falsevalue="" checked="false" label="Clip polyA tail" help="(-T)"/> 190 <param argument="--clipacc" value="70" type="integer" label="clipping accuracy" />
172 <param name="autoclip" type="boolean" truevalue="--autoclip" falsevalue="" checked="false" label="Autoclip unknown 3prime adapter" help="(-Y)"/> 191 <param argument="--polyA" type="boolean" truevalue="--polyA" falsevalue="" checked="false" label="Clip polyA tail" />
173 <param name="hardclip" type="boolean" truevalue="--hardclip" falsevalue="" checked="false" label="Enable hard clipping" help="(-C)"/> 192 <param argument="--autoclip" type="boolean" truevalue="--autoclip" falsevalue="" checked="false" label="Autoclip unknown 3prime adapter"/>
174 <param name="order" type="boolean" truevalue="--order" falsevalue="" checked="false" label="Sorts the output by chromsome and position" help="(-O)"/> 193 <param argument="--hardclip" type="boolean" truevalue="--hardclip" falsevalue="" checked="false" label="Enable hard clipping"/>
175 <param name="differences" type="integer" min="0" value="1" label="search seeds initially with n differences" help="(--differences)"/> 194 <param argument="--order" type="boolean" truevalue="--order" falsevalue="" checked="false" label="Sorts the output by chromsome and position" />
176 <param name="evalue" type="float" min="0" value="5.000000" label="max evalue" help="(--evalue)"/> 195 <param argument="--differences" type="integer" min="0" value="1" label="search seeds initially with n differences"/>
177 <param name="maxinterval" type="integer" min="1" value="100" label="maximum width of a suffix array interval, i.e. a query seed will be omitted if it matches more than n times" help="(--maxinterval)"/> 196 <param argument="--jump" type="integer" value="0" min="0" label="search seeds with jump size" help="(0=automatic) (default:0)?"/>
197 <param argument="--evalue" type="float" min="0" value="5.000000" label="max evalue"/>
198 <param argument="--maxinterval" type="integer" min="1" value="100" label="maximum width of a suffix array interval, i.e. a query seed will be omitted if it matches more than n times"/>
199 <param argument="--nomatchfilename" type="boolean" truevalue="yes" falsevalue="no" checked="false" label="Output unmatched reads"/>
178 </inputs> 200 </inputs>
179 <outputs> 201 <outputs>
180 <data format="sam" name="segemehl_out" label="Read alignments on ${on_string}"/> 202 <data format="sam" name="segemehl_out" label="${tool.name} on ${on_string}"/>
203 <data format="fastq" name="segemehl_outunmatched" label="${tool.name} unaligned reads ${on_string}">
204 <filter>output_unmatched</filter>
205 </data>
181 </outputs> 206 </outputs>
182 <tests> 207 <tests>
183 <test> 208 <test>
184 <param name="genomeSource" value="history" /> 209 <param name="genomeSource" value="history" />
185 <param name="own_reference_genome" value="chr1.fa" /> 210 <param name="own_reference_genome" value="chr1.fa" />
192 <param name="genomeSource" value="history" /> 217 <param name="genomeSource" value="history" />
193 <param name="own_reference_genome" value="chr1.fa" /> 218 <param name="own_reference_genome" value="chr1.fa" />
194 <param name="library" value="single" /> 219 <param name="library" value="single" />
195 <param name="input_query" value="test.fastq" /> 220 <param name="input_query" value="test.fastq" />
196 <param name="splits" value="splits" /> 221 <param name="splits" value="splits" />
197 <param name="minsplicecover" value="40" /> 222 <param name="minsplicecover" value="40" />
223 <param name="nomatchfilename" value="yes" />
198 <output name="segemehl_out" file="testmap2.sam" lines_diff="2" /> 224 <output name="segemehl_out" file="testmap2.sam" lines_diff="2" />
225 <output name="segemehl_outunmatched" file="testmap2.fastq" />
199 </test> 226 </test>
200 </tests> 227 </tests>
201 <help> 228 <help>
202 <![CDATA[ 229 <![CDATA[
203 230