Mercurial > repos > jjohnson > gmap
annotate gsnap.xml @ 6:4f358603ee12 draft default tip
Uploaded v3.0.1c with table fix
| author | peterjc |
|---|---|
| date | Fri, 21 Oct 2016 11:15:09 -0400 |
| parents | 14561eb803a5 |
| children |
| rev | line source |
|---|---|
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
1 <tool id="gsnap" name="GSNAP" version="3.0.1"> |
| 0 | 2 <description>Genomic Short-read Nucleotide Alignment Program</description> |
| 3 <requirements> | |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
4 <requirement type="package" version="2013-05-09">gmap</requirement> |
| 0 | 5 </requirements> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
6 <version_command>gsnap --version</version_command> |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
7 <command detect_errors="exit_code"><![CDATA[ |
| 0 | 8 #import os.path, re |
| 9 gsnap | |
| 10 --nthreads="4" --ordered | |
| 11 #if $refGenomeSource.genomeSource == "gmapdb": | |
| 12 --dir=$refGenomeSource.gmapdb.extra_files_path --db=$refGenomeSource.gmapdb.metadata.db_name | |
| 13 #else: | |
| 14 --dir=$os.path.dirname($refGenomeSource.gmapindex.value) --db=$os.path.basename($refGenomeSource.gmapindex.value) | |
| 15 #end if | |
| 16 #if $refGenomeSource.kmer != None and len($refGenomeSource.kmer.__str__) == 2: | |
| 17 --kmer=$refGenomeSource.kmer | |
| 18 #end if | |
| 19 #if $refGenomeSource.use_splicing.src == 'gmapdb': | |
| 20 #if $refGenomeSource.use_splicing.splicemap != None and len($refGenomeSource.use_splicing.splicemap.__str__) > 0: | |
| 21 -s $refGenomeSource.use_splicing.splicemap.value | |
| 22 #if $computation.trim_mismatch_score.__str__ == '0': | |
| 23 $ambig_splice_noclip | |
| 24 #end if | |
| 25 #end if | |
| 26 #elif $refGenomeSource.use_splicing.src == 'history': | |
| 27 #if $refGenomeSource.use_splicing.splicemap != None and len($refGenomeSource.use_splicing.splicemap.__str__) > 0: | |
| 28 -S $os.path.dirname($refGenomeSource.use_splicing.splicemap) -s $os.path.basename($refGenomeSource.use_splicing.splicemap) | |
| 29 #if $computation.trim_mismatch_score.__str__ == '0': | |
| 30 $ambig_splice_noclip | |
| 31 #end if | |
| 32 #end if | |
| 33 #end if | |
| 34 #if $refGenomeSource.use_snps.src == 'gmapdb': | |
| 35 #if $refGenomeSource.use_snps.snpindex != None and len($refGenomeSource.use_snps.snpindex.__str__) > 0: | |
| 36 -v $refGenomeSource.use_snps.snpindex.value | |
| 37 #end if | |
| 38 #elif $refGenomeSource.use_snps.src == 'history': | |
| 39 #if $refGenomeSource.use_snps.snpindex != None and len($refGenomeSource.use_snps.snpindex.__str__) > 0: | |
| 40 -V $refGenomeSource.use_snps.snpindex.extra_files_path -v $refGenomeSource.use_snps.snpindex.metadata.snps_name | |
| 41 #end if | |
| 42 #end if | |
| 43 #if $refGenomeSource.mode.__str__ != '': | |
| 44 --mode=$refGenomeSource.mode | |
| 45 #end if | |
| 46 #* ## No longer in options as of version 2011-11-30 | |
| 47 #if $mapq_unique_score.__str__ != '': | |
| 48 --mapq-unique-score=$mapq_unique_score | |
| 49 #end if | |
| 50 *# | |
| 51 #if $computation.options == "advanced": | |
| 52 #if $computation.max_mismatches.__str__ != '': | |
| 53 --max-mismatches=$computation.max_mismatches | |
| 54 #end if | |
| 55 $computation.query_unk_mismatch | |
| 56 $computation.genome_unk_mismatch | |
| 57 #if $computation.terminal_threshold.__str__ != '': | |
| 58 --terminal-threshold=$computation.terminal_threshold | |
| 59 #end if | |
| 60 #if $computation.indel_penalty.__str__ != '': | |
| 61 --indel-penalty=$computation.indel_penalty | |
| 62 #end if | |
| 63 #if $computation.indel_endlength.__str__ != '': | |
| 64 --indel-endlength=$computation.indel_endlength | |
| 65 #end if | |
| 66 #if $computation.max_middle_insertions.__str__ != '': | |
| 67 --max-middle-insertions=$computation.max_middle_insertions | |
| 68 #end if | |
| 69 #if $computation.max_middle_deletions.__str__ != '': | |
| 70 --max-middle-deletions=$computation.max_middle_deletions | |
| 71 #end if | |
| 72 #if $computation.max_end_insertions.__str__ != '': | |
| 73 --max-end-insertions=$computation.max_end_insertions | |
| 74 #end if | |
| 75 #if $computation.max_end_deletions.__str__ != '': | |
| 76 --max-end-deletions=$computation.max_end_deletions | |
| 77 #end if | |
| 78 #if $computation.suboptimal_levels.__str__ != '': | |
| 79 --suboptimal-levels=$computation.suboptimal_levels | |
| 80 #end if | |
| 81 #if $computation.adapter_strip.__str__ != '': | |
| 82 --adapter-strip=$computation.adapter_strip | |
| 83 #end if | |
| 84 #if $computation.trim_mismatch_score.__str__ != '': | |
| 85 --trim-mismatch-score=$computation.trim_mismatch_score | |
| 86 #end if | |
| 87 #if $computation.trim_indel_score.__str__ != '': | |
| 88 --trim-indel-score=$computation.trim_indel_score | |
| 89 #end if | |
| 90 ## TODO - do we need these options (Is it tally XOR runlength?): | |
| 91 ## --tallydir= --use-tally=tally | |
| 92 ## --runlengthdir --use-runlength=runlength | |
| 93 #if $computation.use_tally != None and len($computation.use_tally.__str__) > 0: | |
| 94 ##--tallydir $os.path.dirname($computation.use_tally) --use-tally $os.path.basename($computation.use_tally) | |
| 95 --use-tally=$computation.use_tally | |
| 96 #end if | |
| 97 ## gmap options | |
| 98 #if $computation.gmap_mode.__str__ != '' and $computation.gmap_mode.__str__ != 'None': | |
| 99 --gmap-mode='$computation.gmap_mode' | |
| 100 #end if | |
| 101 #if $computation.trigger_score_for_gmap.__str__ != '': | |
| 102 --trigger-score-for-gmap=$computation.trigger_score_for_gmap | |
| 103 #end if | |
| 104 #if $computation.max_gmap_pairsearch.__str__ != '' and $re.search("pairsearch",$computation.gmap_mode): | |
| 105 --max-gmap-pairsearch=$computation.max_gmap_pairsearch | |
| 106 #end if | |
| 107 #if $computation.max_gmap_terminal.__str__ != '' and $re.search("terminal",$computation.gmap_mode): | |
| 108 --max-gmap-terminal=$computation.max_gmap_terminal | |
| 109 #end if | |
| 110 #if $computation.max_gmap_improvement.__str__ != '' and $re.search("improv",$computation.gmap_mode): | |
| 111 --max-gmap-improvement=$computation.max_gmap_improvement | |
| 112 #end if | |
| 113 #if $computation.microexon_spliceprob.__str__ != '': | |
| 114 --microexon-spliceprob=$computation.microexon_spliceprob | |
| 115 #end if | |
| 116 #end if | |
| 117 #if $splicing.options == "advanced": | |
| 118 $splicing.novelsplicing | |
| 119 #if $splicing.localsplicedist.__str__ != '': | |
| 120 --localsplicedist=$splicing.localsplicedist | |
| 121 #end if | |
| 122 #if $splicing.local_splice_penalty.__str__ != '': | |
| 123 --local-splice-penalty=$splicing.local_splice_penalty | |
| 124 #end if | |
| 125 #if $splicing.distant_splice_penalty.__str__ != '': | |
| 126 --distant-splice-penalty=$splicing.distant_splice_penalty | |
| 127 #end if | |
| 128 #if $splicing.local_splice_endlength.__str__ != '': | |
| 129 --local-splice-endlength=$splicing.local_splice_endlength | |
| 130 #end if | |
| 131 #if $splicing.distant_splice_endlength.__str__ != '': | |
| 132 --distant-splice-endlength=$splicing.distant_splice_endlength | |
| 133 #end if | |
| 134 #if $splicing.distant_splice_identity.__str__ != '': | |
| 135 --distant-splice-identity=$splicing.distant_splice_identity | |
| 136 #end if | |
| 137 #end if | |
| 138 #if $output.options == "advanced": | |
| 139 #if $output.npath.__str__ != '': | |
| 140 --npath=$output.npath | |
| 141 #end if | |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
142 #if $output.maxsearch.__str__ != '': |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
143 --maxsearch=$output.maxsearch |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
144 #end if |
| 0 | 145 $output.quiet_if_excessive |
| 146 $output.show_refdiff | |
| 147 $output.clip_overlap | |
| 148 #end if | |
| 149 #if $result.format == "sam": | |
| 150 --format=sam | |
| 151 $result.no_sam_headers | |
| 152 #if $result.read_group_id.__str__.strip != '': | |
| 153 --read-group-id='$result.read_group_id' | |
| 154 #end if | |
| 155 #if $result.read_group_name.__str__ != '': | |
| 156 --read-group-name='$result.read_group_name' | |
| 157 #end if | |
| 158 #if $result.read_group_library.__str__ != '': | |
| 159 --read-group-library='$result.read_group_library' | |
| 160 #end if | |
| 161 #if $result.read_group_platform.__str__ != '': | |
| 162 --read-group-platform='$result.read_group_platform' | |
| 163 #end if | |
| 164 #if $result.quality_shift.__str__ != '': | |
| 165 --quality-shift=$result.quality_shift | |
| 166 #end if | |
| 167 #elif $result.format == "goby": | |
| 168 #if $result.goby_output.__str__ != '': | |
| 169 --goby-output='$result.goby_output' | |
| 170 #end if | |
| 171 #if $result.creads_window_start.__str__ != '': | |
| 172 --creads-window-start=$result.creads_window_start | |
| 173 #end if | |
| 174 #if $result.creads_window_end.__str__ != '': | |
| 175 --creads-window-end=$result.creads_window_end | |
| 176 #end if | |
| 177 $result.creads_complement | |
| 178 #end if | |
| 179 #if $results.split_output == 'yes': | |
| 180 --split-output=gsnap_out | |
| 181 #if $results.fails.choice == 'nofails': | |
| 182 --nofails | |
| 183 #elif $results.fails.choice == 'failsonly': | |
| 184 --failsonly | |
| 185 #end if | |
| 186 $results.fails_as_input | |
| 187 #else | |
| 188 #if $results.fails.choice == 'nofails': | |
| 189 --nofails | |
| 190 #elif $results.fails.choice == 'failsonly': | |
| 191 --failsonly | |
| 192 $results.fails.fails_as_input | |
| 193 #end if | |
| 194 #end if | |
| 195 #if $seq.format == "gsnap_fasta": | |
| 196 $seq.circularinput $seq.gsnap | |
| 197 #else if $seq.format == "fastq": | |
| 198 #if $seq.barcode_length.__str__ != '': | |
| 199 --barcode-length=$seq.barcode_length | |
| 200 #end if | |
| 201 #if $seq.fastq_id_start.__str__ != '': | |
| 202 --fastq-id-start=$seq.fastq_id_start | |
| 203 #end if | |
| 204 #if $seq.fastq_id_end.__str__ != '': | |
| 205 --fastq-id-end=$seq.fastq_id_end | |
| 206 #end if | |
| 207 #if $seq.filter_chastity.__str__ != 'off': | |
| 208 --filter-chastity=$seq.filter_chastity | |
| 209 #end if | |
| 210 #if $seq.paired.ispaired.__str__ == 'yes': | |
| 211 #if $seq.paired.pairmax_dna.__str__ != '': | |
| 212 --pairmax-dna=$seq.paired.pairmax_dna | |
| 213 #end if | |
| 214 #if $seq.paired.pairmax_rna.__str__ != '': | |
| 215 --pairmax-rna=$seq.paired.pairmax_rna | |
| 216 #end if | |
| 217 #if $seq.paired.pairexpect.__str__ != '': | |
| 218 --pairexpect=$seq.paired.pairexpect | |
| 219 #end if | |
| 220 #if $seq.paired.pairdev.__str__ != '': | |
| 221 --pairdev=$seq.paired.pairdev | |
| 222 #end if | |
| 223 $seq.fastq $seq.paired.fastq | |
| 224 #else | |
| 225 $seq.fastq | |
| 226 #end if | |
| 227 #end if | |
| 228 #if $results.split_output == 'yes': | |
| 229 2> $gsnap_stderr | |
| 230 #else: | |
| 231 #if $results.fails.choice.__str__ == 'failsonly' and $results.fails.fails_as_input.__str__ != '': | |
| 232 2> $gsnap_stderr > $gsnap_fq | |
| 233 #else | |
| 234 2> $gsnap_stderr > $gsnap_out | |
| 235 #end if | |
| 236 #end if | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
237 ]]></command> |
| 0 | 238 <inputs> |
| 239 <!-- Input data --> | |
| 240 <conditional name="seq"> | |
| 241 <param name="format" type="select" label="<H2>Input Sequences</H2>Select the input format" help=""> | |
| 242 <option value="fastq">Fastq</option> | |
| 243 <!-- | |
| 244 <option value="goby">Goby compact-reads</option> | |
| 245 --> | |
| 246 <option value="gsnap_fasta">GNSAP fasta</option> | |
| 247 </param> | |
| 248 <when value="fastq"> | |
| 249 <param name="fastq" type="data" format="fastq" label="Select a fastq dataset" /> | |
| 250 <conditional name="paired"> | |
| 251 <param name="ispaired" type="boolean" truevalue="yes" falsevalue="no" checked="false" label="Use Paired Reads?"/> | |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
252 <when value="no"> |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
253 <!-- |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
254 Could allow multiple input fastq and set the force-single-end flag |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
255 force-single-end When multiple FASTQ files are provided on the command line, GSNAP assumes |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
256 they are match paired-end files. This flag treats each file as single-end. |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
257 --> |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
258 </when> |
| 0 | 259 <when value="yes"> |
| 260 <param name="fastq" type="data" format="fastq" label="Select the paired reads reverse dataset" /> | |
| 261 <param name="orientation" type="select" label="Orientation of paired-end reads" help=""> | |
| 262 <option value="FR">fwd-rev, typical Illumina default</option> | |
| 263 <option value="RF">rev-fwd, for circularized inserts</option> | |
| 264 <option value="FF">fwd-fwd, same strand</option> | |
| 265 </param> | |
| 266 <param name="pairmax_dna" type="integer" value="" optional="true" label="Max total genomic length for DNA-Seq paired reads, or other reads without splicing (default 1000)." help="Used if no splice file is provided and novelsplicing is off."/> | |
| 267 <param name="pairmax_rna" type="integer" value="" optional="true" label="Max total genomic length for RNA-Seq paired reads, or other reads that could have a splice (default 200000)." help="Used when novel splicing is specified or a splice file is provided. Should probably match the value for localsplicedist."/> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
268 <param name="pairexpect" type="integer" value="" optional="true" label="Expected paired-end length" |
| 0 | 269 help="Used for calling splices in medial part of paired-end reads (default 200)"/> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
270 <param name="pairdev" type="integer" value="" optional="true" label="Allowable deviation from expected paired-end length" |
| 0 | 271 help="Used for calling splices in medial part of paired-end reads (default 25)"/> |
| 272 </when> | |
| 273 </conditional> | |
| 274 <param name="barcode_length" type="integer" value="" optional="true" label="Amount of barcode to remove from start of read (default 0)" /> | |
| 275 <param name="fastq_id_start" type="integer" value="" optional="true" label="Starting field of identifier in FASTQ header, whitespace-delimited, starting from 1" /> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
276 <param name="fastq_id_end" type="integer" value="" optional="true" label="Ending field of identifier in FASTQ header, whitespace-delimited, starting from 1" |
| 0 | 277 help="Examples: |
| 278 <br>@HWUSI-EAS100R:6:73:941:1973#0/1 | |
| 279 <br> . start=1, end=1 (default) => identifier is HWUSI-EAS100R:6:73:941:1973#0/1 | |
| 280 <br>@SRR001666.1 071112_SLXA-EAS1_s_7:5:1:817:345 length=36 | |
| 281 <br> . start=1, end=1 => identifier is SRR001666.1 | |
| 282 <br> . start=2, end=2 => identifier is 071112_SLXA-EAS1_s_7:5:1:817:345 | |
| 283 <br> . start=1, end=2 => identifier is SRR001666.1 071112_SLXA-EAS1_s_7:5:1:817:345" | |
| 284 /> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
285 <param name="filter_chastity" type="select" label="Skip reads marked by the Illumina chastity program" |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
286 help="String after the accession having a 'Y' after the first colon, like this: |
| 0 | 287 <br>@accession 1:Y:0:CTTGTA |
| 288 <br>where the 'Y' signifies filtering by chastity. | |
| 289 <br> For 'either', a 'Y' on either end of a paired-end read will be filtered. | |
| 290 <br> For 'both', a 'Y' is required on both ends of a paired-end read (or on the only end of a single-end read)" | |
| 291 > | |
| 292 <option value="off">off - no filtering</option> | |
| 293 <option value="either">either - a 'Y' on either end of a paired-end read</option> | |
| 294 <option value="both">both - a 'Y' is required on both ends of a paired-end read or the only end of a single-end read</option> | |
| 295 </param> | |
| 296 </when> | |
| 297 <!-- | |
| 298 <when value="goby"> | |
| 299 </when> | |
| 300 --> | |
| 301 <when value="gsnap_fasta"> | |
| 302 <param name="gsnap" type="data" format="fasta" label="Select a single-end dataset" help="GSNAP fasta must have the sequence entirely on one line, a second line is interpreted as the paired-end sequence"/> | |
| 303 <param name="circularinput" type="boolean" checked="false" truevalue="--circular-input=true" falsevalue="" label="Circular-end data (paired reads are on same strand)"/> | |
| 304 </when> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
305 |
| 0 | 306 </conditional> |
| 307 <!-- No longer in options as of version 2011-11-30 | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
308 <param name="mapq_unique_score" type="integer" value="" optional="true" label="MAPQ score threshold" |
| 0 | 309 help="For multiple results, consider as a unique result if only one of the results has a MAPQ score equal or greater than this |
| 310 (if not selected, then reports all multiple results, up to npaths)" /> | |
| 311 --> | |
| 312 | |
| 313 <!-- GMAPDB for alignment --> | |
| 314 <conditional name="refGenomeSource"> | |
| 315 <param name="genomeSource" type="select" label="<HR><H2>Align To</H2>Will you select a reference genome from your history or use a built-in index?" help="Built-ins were indexed using default options"> | |
| 316 <option value="indexed">Use a built-in index</option> | |
| 317 <option value="gmapdb">Use a gmapdb from your history</option> | |
| 318 </param> | |
| 319 <when value="indexed"> | |
| 320 <param name="gmapindex" type="select" label="Select a reference genome" help="if your genome of interest is not listed - contact Galaxy team"> | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
321 <options from_data_table="gmap_indices"> |
| 0 | 322 <column name="uid" index="0" /> |
| 323 <column name="dbkey" index="1" /> | |
| 324 <column name="name" index="2" /> | |
| 325 <column name="kmers" index="3" /> | |
| 326 <column name="maps" index="4" /> | |
| 327 <column name="snps" index="5" /> | |
| 328 <column name="value" index="6" /> | |
| 329 </options> | |
| 330 </param> | |
| 331 | |
| 332 <param name="kmer" type="select" data_ref="gmapindex" label="kmer size" help="Defaults to highest available kmer size"> | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
333 <options from_data_table="gmap_indices"> |
| 0 | 334 <column name="name" index="3"/> |
| 335 <column name="value" index="3"/> | |
| 336 <filter type="param_value" ref="gmapindex" column="6"/> | |
| 337 <filter type="multiple_splitter" column="3" separator=","/> | |
| 338 <filter type="add_value" name="" value=""/> | |
| 339 <filter type="sort_by" column="3"/> | |
| 340 </options> | |
| 341 </param> | |
| 342 | |
| 343 <param name="mode" type="select" label="Alignment mode" help="Assumes cmetindex and atoiindex were run on the gmap datatbase."> | |
| 344 <option value="">standard</option> | |
| 345 <option value="cmet-stranded">cmet-stranded for bisulfite-treated DNA reads (tolerance to C-to-T changes)</option> | |
| 346 <option value="cmet-nonstranded">cmet-nonstranded for bisulfite-treated DNA reads (tolerance to C-to-T changes)</option> | |
| 347 <option value="atoi-stranded">atoi-stranded for RNA-editing tolerance (A-to-G changes)</option> | |
| 348 <option value="atoi-nonstranded">atoi-nonstranded for RNA-editing tolerance (A-to-G changes)</option> | |
| 349 </param> | |
| 350 | |
| 351 <conditional name="use_splicing"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
352 <param name="src" type="select" label="<HR>Known Splicesite and Introns" |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
353 help="Look for splicing involving known sites or known introns at short or long distances |
| 0 | 354 See README instructions for the distinction between known sites and known introns"> |
| 355 <option value="none" selected="true">None</option> | |
| 356 <option value="gmapdb">From the GMAP Database</option> | |
| 357 <option value="history">A Map in your history</option> | |
| 358 </param> | |
| 359 <when value="none"/> | |
| 360 <when value="history"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
361 <param name="splicemap" type="data" format="splicesites.iit,introns.iit" label="Select a splicesite map" |
| 0 | 362 help="built with GMAP IIT"/> |
| 363 </when> | |
| 364 <when value="gmapdb"> | |
| 365 <param name="splicemap" type="select" data_ref="gmapindex" label="Use map for splicing involving known sites or known introns" help=""> | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
366 <options from_data_table="gmap_indices"> |
| 0 | 367 <column name="name" index="4"/> |
| 368 <column name="value" index="4"/> | |
| 369 <filter type="param_value" ref="gmapindex" column="6"/> | |
| 370 <filter type="multiple_splitter" column="4" separator=","/> | |
| 371 <filter type="add_value" name="" value=""/> | |
| 372 <filter type="sort_by" column="4"/> | |
| 373 </options> | |
| 374 </param> | |
| 375 </when> | |
| 376 </conditional> | |
| 377 | |
| 378 <conditional name="use_snps"> | |
| 379 <param name="src" type="select" label="<HR>Known SNPs" help="for SNP tolerant alignments"> | |
| 380 <option value="none" selected="true">None</option> | |
| 381 <option value="gmapdb">From the GMAP Database</option> | |
| 382 <option value="history">A SNP Index in your history</option> | |
| 383 </param> | |
| 384 <when value="none"/> | |
| 385 <when value="history"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
386 <param name="snpindex" type="data" format="gmapsnpindex" label="Select a snpindex" |
| 0 | 387 help="built with GMAP SNP Index"/> |
| 388 </when> | |
| 389 <when value="gmapdb"> | |
| 390 <param name="snpindex" type="select" data_ref="gmapindex" label="Use database containing known SNPs" help=""> | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
391 <options from_data_table="gmap_indices"> |
| 0 | 392 <column name="name" index="5"/> |
| 393 <column name="value" index="5"/> | |
| 394 <filter type="param_value" ref="gmapindex" column="6"/> | |
| 395 <filter type="multiple_splitter" column="5" separator=","/> | |
| 396 <filter type="add_value" name="" value=""/> | |
| 397 <filter type="sort_by" column="5"/> | |
| 398 </options> | |
| 399 </param> | |
| 400 </when> | |
| 401 </conditional> | |
| 402 | |
| 403 </when> | |
| 404 <when value="gmapdb"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
405 <param name="gmapdb" type="data" format="gmapdb" label="Select a gmapdb" |
| 0 | 406 help="A GMAP database built with GMAP Build"/> |
| 407 <param name="kmer" type="select" data_ref="gmapdb" label="kmer size" help="Defaults to highest available kmer size"> | |
| 408 <options> | |
| 409 <filter type="data_meta" ref="gmapdb" key="kmers" multiple="True" separator=","/> | |
| 410 </options> | |
| 411 </param> | |
| 412 | |
| 413 <param name="mode" type="select" label="Alignment mode" help="Assumes cmetindex and atoiindex were run on the gmap datatbase."> | |
| 414 <option value="">standard</option> | |
| 415 <option value="cmet-stranded">cmet-stranded for bisulfite-treated DNA reads (tolerance to C-to-T changes)</option> | |
| 416 <option value="cmet-nonstranded">cmet-nonstranded for bisulfite-treated DNA reads (tolerance to C-to-T changes)</option> | |
| 417 <option value="atoi-stranded">atoi-stranded for RNA-editing tolerance (A-to-G changes)</option> | |
| 418 <option value="atoi-nonstranded">atoi-nonstranded for RNA-editing tolerance (A-to-G changes)</option> | |
| 419 </param> | |
| 420 | |
| 421 <conditional name="use_splicing"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
422 <param name="src" type="select" label="<HR>Known Splicesite and Introns" |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
423 help="Look for splicing involving known sites or known introns at short or long distances |
| 0 | 424 See README instructions for the distinction between known sites and known introns"> |
| 425 <option value="none" selected="true">None</option> | |
| 426 <option value="gmapdb">From the GMAP Database</option> | |
| 427 <option value="history">A Map in your history</option> | |
| 428 </param> | |
| 429 <when value="none"/> | |
| 430 <when value="history"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
431 <param name="splicemap" type="data" format="splicesites.iit,introns.iit" label="Select a splicesite map" |
| 0 | 432 help="built with GMAP IIT"/> |
| 433 <param name="ambig_splice_noclip" type="boolean" checked="false" truevalue="--ambig-splice-noclip" falsevalue="" label="Do not clip at ambiguous splice sites" | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
434 help="For ambiguous known splicing at ends of the read, do not clip at the splice site, but extend instead into the intron. |
| 0 | 435 This flag makes sense only if you are trying to eliminate all soft clipping with --trim-mismatch-score=0"/> |
| 436 </when> | |
| 437 <when value="gmapdb"> | |
| 438 <param name="splicemap" type="select" data_ref="gmapdb" label="Use map for splicing involving known sites or known introns" help=""> | |
| 439 <options> | |
| 440 <filter type="data_meta" ref="gmapdb" key="maps" multiple="True"/> | |
| 441 </options> | |
| 442 </param> | |
| 443 <param name="ambig_splice_noclip" type="boolean" checked="false" truevalue="--ambig-splice-noclip" falsevalue="" label="Do not clip at ambiguous splice sites" | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
444 help="For ambiguous known splicing at ends of the read, do not clip at the splice site, but extend instead into the intron. |
| 0 | 445 This flag makes sense only if you are trying to eliminate all soft clipping with --trim-mismatch-score=0"/> |
| 446 </when> | |
| 447 </conditional> | |
| 448 | |
| 449 <conditional name="use_snps"> | |
| 450 <param name="src" type="select" label="<HR>Known SNPs" help="for SNP tolerant alignments"> | |
| 451 <option value="none" selected="true">None</option> | |
| 452 <option value="gmapdb">From the GMAP Database</option> | |
| 453 <option value="history">A SNP Index in your history</option> | |
| 454 </param> | |
| 455 <when value="none"/> | |
| 456 <when value="history"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
457 <param name="snpindex" type="data" format="gmapsnpindex" label="Select a snpindex" |
| 0 | 458 help="built with GMAP SNP Index"/> |
| 459 </when> | |
| 460 <when value="gmapdb"> | |
| 461 <param name="snpindex" type="select" data_ref="gmapdb" label="Use database containing known SNPs" help=""> | |
| 462 <options> | |
| 463 <filter type="data_meta" ref="gmapdb" key="snps" multiple="True" separator=","/> | |
| 464 </options> | |
| 465 </param> | |
| 466 </when> | |
| 467 </conditional> | |
| 468 | |
| 469 </when> | |
| 470 </conditional> | |
| 471 | |
| 472 <!-- Computation options --> | |
| 473 <conditional name="computation"> | |
| 474 <param name="options" type="select" label="<HR>Computational Settings" help=""> | |
| 475 <option value="default">Use default settings</option> | |
| 476 <option value="advanced">Set Computation Options</option> | |
| 477 </param> | |
| 478 <when value="default"/> | |
| 479 <when value="advanced"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
480 <param name="max_mismatches" type="float" value="" optional="true" label="Maximum number of mismatches allowed (uses default when negative)" |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
481 help="Maximum number of mismatches allowed (if not specified, then |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
482 defaults to the ultrafast level of ((readlength+index_interval-1)/kmer - 2)) |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
483 (By default, the genome index interval is 3, but this can be changed |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
484 by providing a different value for -q to gmap_build when processing the genome.) |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
485 If specified between 0.0 and 1.0, then treated as a fraction |
| 0 | 486 of each read length. Otherwise, treated as an integral number |
| 487 of mismatches (including indel and splicing penalties) | |
| 488 For RNA-Seq, you may need to increase this value slightly | |
| 489 to align reads extending past the ends of an exon."> | |
| 490 <validator type="in_range" message="The mismatches must >= 0." min="0."/> | |
| 491 </param> | |
| 492 <param name="query_unk_mismatch" type="boolean" checked="false" truevalue="--query-unk-mismatch=1" falsevalue="" label="Count unknown (N) characters in the query as a mismatch"/> | |
| 493 <param name="genome_unk_mismatch" type="boolean" checked="true" truevalue="" falsevalue="--genome-unk-mismatch=0" label="Count unknown (N) characters in the genome as a mismatch"/> | |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
494 <param name="maxsearch" type="integer" value="" optional="true" label="Maximum number of alignments to find (default 1000)" |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
495 help="Must be larger than paths, which is the number to report. |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
496 Keeping this number large will allow for random selection among multiple alignments. |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
497 Reducing this number can speed up the program. "/> |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
498 |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
499 <param name="terminal_threshold" type="integer" value="" optional="true" label="Threshold for searching for a terminal alignment" |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
500 help="Threshold for searching for a terminal alignment (from one end of the |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
501 read to the best possible position at the other end) (default 2 |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
502 for standard, atoi-stranded, and atoi-nonstranded mode; default 100 |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
503 for cmet-stranded and cmet-nonstranded mode). |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
504 For example, if this value is 2, then if GSNAP finds an exact or |
| 0 | 505 1-mismatch alignment, it will not try to find a terminal alignment. |
| 506 Note that this default value may not be low enough if you want to | |
| 507 obtain terminal alignments for very short reads, although such reads | |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
508 probably don't have enough specificity for terminal alignments anyway. |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
509 To turn off terminal alignments, set this to a high value, greater |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
510 than the value for max-mismatches. |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
511 "/> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
512 <param name="indel_penalty" type="integer" value="" optional="true" label="Penalty for an indel (default 2)" |
| 0 | 513 help="Counts against mismatches allowed. To find indels, make indel-penalty less than or equal to max-mismatches. A value < 2 can lead to false positives at read ends" /> |
| 514 <param name="indel_endlength" type="integer" value="" optional="true" label="Minimum length at end required for indel alignments (default 4)" /> | |
| 515 <param name="max_middle_insertions" type="integer" value="" optional="true" label="Maximum number of middle insertions allowed (default 9)" /> | |
| 516 <param name="max_middle_deletions" type="integer" value="" optional="true" label="Maximum number of middle deletions allowed (default 30)" /> | |
| 517 <param name="max_end_insertions" type="integer" value="" optional="true" label="Maximum number of end insertions allowed (default 3)" /> | |
| 518 <param name="max_end_deletions" type="integer" value="" optional="true" label="Maximum number of end deletions allowed (default 6)" /> | |
| 519 <param name="suboptimal_levels" type="integer" value="" optional="true" label="Report suboptimal hits beyond best hit (default 0)" | |
| 520 help="All hits with best score plus suboptimal-levels are reported" /> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
521 <param name="adapter_strip" type="select" label="Method for removing adapters from reads" |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
522 help="Default is 'off'. To turn on, specify 'paired', which removes adapters |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
523 from paired-end reads if they appear to be present."> |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
524 <option value="paired">paired</option> |
|
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
525 <option value="off" selected="true">off</option> |
| 0 | 526 </param> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
527 <param name="trim_mismatch_score" type="integer" value="" optional="true" label="Score to use for mismatches when trimming at ends (default is -3)" |
| 0 | 528 help="to turn off trimming, specify 0 (Warning: turning trimming off will give false positive mismatches at the ends of reads)"/> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
529 <param name="trim_indel_score" type="integer" value="" optional="true" label="Score to use for indels when trimming at ends (default is -4)" |
| 0 | 530 help="to turn off trimming, specify 0 (Warning: turning trimming off will give false positive indels at the ends of reads)"/> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
531 <param name="use_tally" type="data" format="tally.iit" optional="true" label="Select a tally IIT file to resolve concordant multiple results" |
| 0 | 532 help="generated by gsnap_tally and iit_store"/> |
| 533 | |
| 534 <!-- | |
| 535 tallydir=STRING Directory for tally IIT file to resolve concordant multiple results (default is | |
| 536 location of genome index files specified using -D and -d). Note: can | |
| 537 just give full path name to use-tally instead. | |
| 538 use-tally=STRING Use this tally IIT file to resolve concordant multiple results | |
| 539 runlengthdir=STRING Directory for runlength IIT file to resolve concordant multiple results (default is | |
| 540 location of genome index files specified using -D and -d). Note: can | |
| 541 just give full path name to use-runlength instead. | |
| 542 use-runlength=STRING Use this runlength IIT file to resolve concordant multiple results | |
| 543 --> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
544 |
| 0 | 545 <!-- Options for GMAP alignment within GSNAP --> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
546 <param name="gmap_mode" type="select" multiple="true" optional="true" display="checkboxes" label="Cases to use GMAP for complex alignments containing multiple splices or indels" |
| 0 | 547 help="Default: pairsearch,terminal,improve"> |
| 548 <option value="pairsearch" selected="true">pairsearch</option> | |
|
2
f6ba0f12cca2
Untested work-in-progress GMAP wrappers v3.0.0, from JJ back in June 2013
peterjc
parents:
1
diff
changeset
|
549 <option value="indel_knownsplice" selected="true">indel_knownsplice</option> |
| 0 | 550 <option value="terminal" selected="true">terminal</option> |
| 551 <option value="improve" selected="true">improve</option> | |
| 552 </param> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
553 <param name="trigger_score_for_gmap" type="integer" value="" optional="true" label="GMAP pairsearch threshold (default 5)" |
| 0 | 554 help="Try GMAP pairsearch on nearby genomic regions if best score (the total of both ends if paired-end) exceeds this value (default 5)" /> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
555 <param name="max_gmap_pairsearch" type="integer" value="" optional="true" label="GMAP pairsearch threshold (default 3)" |
| 0 | 556 help="Perform GMAP pairsearch on nearby genomic regions up to this many candidate ends (default 3)." /> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
557 <param name="max_gmap_terminal" type="integer" value="" optional="true" label="GMAP terminal threshold (default 3)" |
| 0 | 558 help="Perform GMAP terminal on nearby genomic regions up to this many candidate ends (default 3)." /> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
559 <param name="max_gmap_improvement" type="integer" value="" optional="true" label="GMAP improvement threshold (default 3)" |
| 0 | 560 help="Perform GMAP improvement on nearby genomic regions up to this many candidate ends (default 3)." /> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
561 <param name="microexon_spliceprob" type="float" value="" optional="true" label="GMAP microexons threshold (default .90)" |
| 0 | 562 help="Allow microexons only if one of the splice site probabilities is greater than this value." > |
| 563 <validator type="in_range" message="The microexons probability must be between 0. and 1." min="0." max="1."/> | |
| 564 </param> | |
| 565 </when> | |
| 566 </conditional> | |
| 567 | |
| 568 <conditional name="splicing"> | |
| 569 <param name="options" type="select" label="<HR>Splicing options for RNA-Seq" help=""> | |
| 570 <option value="default">Use default settings</option> | |
| 571 <option value="advanced">Set Splicing Options</option> | |
| 572 </param> | |
| 573 <when value="default"/> | |
| 574 <when value="advanced"> | |
| 575 <!-- Splicing options for RNA-Seq --> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
576 <!-- use-splicing This should be either a select list from the gmapdb maps or a data type using splicesdir and use-splicing --> |
| 0 | 577 <!-- Neither novel splicing (-N) nor known splicing (-s) turned on => assume reads are DNA-Seq (genomic) --> |
| 578 <param name="novelsplicing" type="boolean" checked="false" truevalue="--novelsplicing=1" falsevalue="" label="Look for novel splicing "/> | |
| 579 <param name="localsplicedist" type="integer" value="" optional="true" label="Definition of local novel splicing event (default 200000)"/> | |
| 580 <param name="local_splice_penalty" type="integer" value="" optional="true" label="Penalty for a local splice (default 0). Counts against mismatches allowed"/> | |
| 581 <param name="distant_splice_penalty" type="integer" value="" optional="true" label="Penalty for a distant splice (default 3). Counts against mismatches allowed" | |
| 582 help="A distant splice is one where the intron length exceeds the value of localsplicedist or is an | |
| 583 inversion, scramble, or translocation between two different chromosomes. Counts against mismatches allowed"/> | |
| 584 <param name="distant_splice_endlength" type="integer" value="" optional="true" label="Minimum length at end required for distant spliced alignments" | |
| 585 help="(default 16, min is the kmer length)"/> | |
| 586 <param name="shortend_splice_endlength" type="integer" value="" optional="true" label="Minimum length at end required for short-end spliced alignments" | |
| 587 help="(default 2, but unless known splice sites are provided, GSNAP may still need the end length to be the value of kmer size to find a given splice"/> | |
| 588 <param name="distant_splice_identity" type="float" value="" optional="true" label="Minimum identity at end required for distant spliced alignments (default 0.95)"/> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
589 <param name="antistranded_penalty" type="integer" value="" optional="true" label="Penalty for antistranded splicing when using stranded RNA-Seq protocols" |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
590 help="A positive value, such as 1, expects antisense on the first read and sense on the second read. |
| 0 | 591 Default is 0, which treats sense and antisense equally well"/> |
| 592 </when> | |
| 593 </conditional> | |
| 594 | |
| 595 <!-- Output data --> | |
| 596 <conditional name="output"> | |
| 597 <param name="options" type="select" label="<HR><H2>Output</H2>Output options for RNA-Seq" help=""> | |
| 598 <option value="default">Use default settings</option> | |
| 599 <option value="advanced">Set Output Options</option> | |
| 600 </param> | |
| 601 <when value="default"/> | |
| 602 <when value="advanced"> | |
| 603 <param name="npath" type="integer" value="" optional="true" label="Maximum number of paths to print (default 100)"/> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
604 <param name="quiet_if_excessive" type="boolean" checked="false" truevalue="--quiet-if-excessive" falsevalue="" label="Quiet if Excessive" |
| 0 | 605 help="If more than maximum number of paths are found, then nothing is printed."/> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
606 <param name="show_refdiff" type="boolean" checked="false" truevalue="--show-refdiff" falsevalue="" label="Show SNP-tolerant alignment" |
| 0 | 607 help="For GSNAP output in SNP-tolerant alignment, shows all differences relative to the reference genome as lower case (otherwise, it shows all differences relative to both the reference and alternate genome)"/> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
608 <param name="clip_overlap" type="boolean" checked="false" truevalue="--clip-overlap" falsevalue="" label="Clip Overlap" |
| 0 | 609 help="For paired-end reads whose alignments overlap, clip the overlapping region."/> |
| 610 </when> | |
| 611 </conditional> | |
| 612 <conditional name="result"> | |
| 613 <param name="format" type="select" label="Select the output format" help=""> | |
| 614 <option value="sam">SAM</option> | |
| 615 <!-- goby should only be an option if the input is in goby format | |
| 616 <option value="goby">Goby</option> | |
| 617 --> | |
| 618 <option value="gsnap">GSNAP default output</option> | |
| 619 </param> | |
| 620 <when value="gsnap"> | |
| 621 </when> | |
| 622 <when value="sam"> | |
| 623 <param name="no_sam_headers" type="boolean" truevalue="--no-sam-headers" falsevalue="" checked="false" label="Do not print headers beginning with '@'"/> | |
| 624 <param name="read_group_id" type="text" value="" optional="true" label="Value to put into read-group id (RG-ID) field"/> | |
| 625 <param name="read_group_name" type="text" value="" optional="true" label="Value to put into read-group name (RG-SM) field"/> | |
| 626 <param name="read_group_library" type="text" value="" optional="true" label="Value to put into read-group library (RG-LB) field"/> | |
| 627 <param name="read_group_platform" type="text" value="" optional="true" label="Value to put into read-group library platform (RG-PL) field"/> | |
| 628 <param name="quality_shift" type="integer" value="" optional="true" label="Shift FASTQ quality scores by this amount in SAM output (default -31)"/> | |
| 629 </when> | |
| 630 <!-- | |
| 631 <when value="goby"> | |
| 632 <param name="goby_output" type="text" value="" label="Basename for Goby output files"/> | |
| 633 <param name="creads_window_start" type="integer" value="" optional="true" label="Compact reads window start (default: 0=start of file)"/> | |
| 634 <param name="creads_window_end" type="integer" value="" optional="true" label="Compact reads window end (default: 0=end of file)"/> | |
| 635 <param name="creads_complement" type="boolean" truevalue="-\-creads-complement" falsevalue="" checked="false" label="Complement read sequences (without reversing)"/> | |
| 636 </when> | |
| 637 --> | |
| 638 </conditional> | |
| 639 <!-- TODO combine fails and split_output --> | |
| 640 | |
| 641 <conditional name="results"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
642 <param name="split_output" type="select" label="<HR>Split outputs" |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
643 help="Separate outputs for: nomapping, halfmapping_uniq, halfmapping_mult, unpaired_uniq, unpaired_mult, paired_uniq, paired_mult, concordant_uniq, and concordant_mult results"> |
| 0 | 644 <option value="no">no</option> |
| 645 <option value="yes">yes</option> | |
| 646 </param> | |
| 647 <when value="no"> | |
| 648 <conditional name="fails"> | |
| 649 <param name="choice" type="select" label="How to deal with fails" help=""> | |
| 650 <option value="default">default - include them in results</option> | |
| 651 <option value="nofails">nofails - exclude fails from results</option> | |
| 652 <option value="failsonly">failsonly - only output failing results</option> | |
| 653 </param> | |
| 654 <when value="default"/> | |
| 655 <when value="nofails"/> | |
| 656 <when value="failsonly"> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
657 <param name="fails_as_input" type="boolean" truevalue="--fails-as-input" falsevalue="" checked="false" label="Print completely failed alignments as input FASTA or FASTQ format" |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
658 help=""/> |
| 0 | 659 </when> |
| 660 </conditional> | |
| 661 </when> | |
| 662 <when value="yes"> | |
| 663 <conditional name="fails"> | |
| 664 <param name="choice" type="select" label="How to deal with fails" help=""> | |
| 665 <option value="default">default - include them in results</option> | |
| 666 <option value="nofails">nofails - exclude fails from results</option> | |
| 667 <option value="failsonly">failsonly - only output failing results</option> | |
| 668 </param> | |
| 669 <when value="default"/> | |
| 670 <when value="nofails"/> | |
| 671 <when value="failsonly"/> | |
| 672 </conditional> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
673 <param name="fails_as_input" type="boolean" truevalue="--fails-as-input" falsevalue="" checked="false" label="Print completely failed alignments as input FASTA or FASTQ format" |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
674 help=""/> |
| 0 | 675 </when> |
| 676 </conditional> | |
| 677 | |
| 678 </inputs> | |
| 679 <outputs> | |
| 680 <data format="txt" name="gsnap_stderr" label="${tool.name} on ${on_string}: gsnap.log"/> | |
| 681 | |
| 682 <data format="txt" name="gsnap_out" label="${tool.name} on ${on_string} ${result.format}" > | |
| 683 <filter>(results['split_output'] == 'no' and (results['fails']['choice'] != 'failsonly' or results['fails']['fails_as_input'] == False))</filter> | |
| 684 <change_format> | |
| 685 <when input="result['format']" value="sam" format="sam"/> | |
| 686 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 687 </change_format> | |
| 688 </data> | |
| 689 | |
| 690 <data format="fastq" name="gsnap_fq" label="${tool.name} on ${on_string} fails.fq" > | |
| 691 <filter>(results['split_output'] == 'no' and results['fails']['choice'] == 'failsonly' and results['fails']['fails_as_input'] == True)</filter> | |
| 692 </data> | |
| 693 | |
| 694 <!-- nomapping, halfmapping_uniq, halfmapping_mult, unpaired_uniq, unpaired_mult, paired_uniq, paired_mult, concordant_uniq, concordant_mult --> | |
| 695 | |
| 696 <data format="txt" name="unpaired_mult" label="${tool.name} on ${on_string} unpaired_mult.${result.format}" from_work_dir="gsnap_out.unpaired_mult"> | |
| 697 <filter>(results['split_output'] == 'yes')</filter> | |
| 698 <change_format> | |
| 699 <when input="result['format']" value="sam" format="sam"/> | |
| 700 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 701 </change_format> | |
| 702 </data> | |
| 703 <data format="txt" name="unpaired_uniq" label="${tool.name} on ${on_string} unpaired_uniq.${result.format}" from_work_dir="gsnap_out.unpaired_uniq"> | |
| 704 <filter>(results['split_output'] == 'yes')</filter> | |
| 705 <change_format> | |
| 706 <when input="result['format']" value="sam" format="sam"/> | |
| 707 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 708 </change_format> | |
| 709 </data> | |
| 710 <data format="txt" name="unpaired_transloc" label="${tool.name} on ${on_string} unpaired_transloc.${result.format}" from_work_dir="gsnap_out.unpaired_transloc"> | |
| 711 <filter>(results['split_output'] == 'yes')</filter> | |
| 712 <change_format> | |
| 713 <when input="result['format']" value="sam" format="sam"/> | |
| 714 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 715 </change_format> | |
| 716 </data> | |
| 717 <data format="txt" name="halfmapping_mult" label="${tool.name} on ${on_string} halfmapping_mult.${result.format}" from_work_dir="gsnap_out.halfmapping_mult"> | |
| 718 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 719 <change_format> | |
| 720 <when input="result['format']" value="sam" format="sam"/> | |
| 721 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 722 </change_format> | |
| 723 </data> | |
| 724 <data format="txt" name="halfmapping_uniq" label="${tool.name} on ${on_string} halfmapping_uniq.${result.format}" from_work_dir="gsnap_out.halfmapping_uniq"> | |
| 725 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 726 <change_format> | |
| 727 <when input="result['format']" value="sam" format="sam"/> | |
| 728 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 729 </change_format> | |
| 730 </data> | |
| 731 <data format="txt" name="halfmapping_transloc" label="${tool.name} on ${on_string} halfmapping_transloc.${result.format}" from_work_dir="gsnap_out.halfmapping_transloc"> | |
| 732 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 733 <change_format> | |
| 734 <when input="result['format']" value="sam" format="sam"/> | |
| 735 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 736 </change_format> | |
| 737 </data> | |
| 738 <data format="txt" name="paired_mult" label="${tool.name} on ${on_string} paired_mult.${result.format}" from_work_dir="gsnap_out.paired_mult"> | |
| 739 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 740 <change_format> | |
| 741 <when input="result['format']" value="sam" format="sam"/> | |
| 742 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 743 </change_format> | |
| 744 </data> | |
| 745 <data format="txt" name="paired_uniq" label="${tool.name} on ${on_string} paired_uniq.${result.format}" from_work_dir="gsnap_out.paired_uniq"> | |
| 746 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 747 <change_format> | |
| 748 <when input="result['format']" value="sam" format="sam"/> | |
| 749 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 750 </change_format> | |
| 751 </data> | |
| 752 <data format="txt" name="paired_transloc" label="${tool.name} on ${on_string} paired_transloc.${result.format}" from_work_dir="gsnap_out.paired_transloc"> | |
| 753 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 754 <change_format> | |
| 755 <when input="result['format']" value="sam" format="sam"/> | |
| 756 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 757 </change_format> | |
| 758 </data> | |
| 759 | |
| 760 <data format="txt" name="concordant_mult" label="${tool.name} on ${on_string} concordant_mult.${result.format}" from_work_dir="gsnap_out.concordant_mult"> | |
| 761 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 762 <change_format> | |
| 763 <when input="result['format']" value="sam" format="sam"/> | |
| 764 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 765 </change_format> | |
| 766 </data> | |
| 767 <data format="txt" name="concordant_uniq" label="${tool.name} on ${on_string} concordant_uniq.${result.format}" from_work_dir="gsnap_out.concordant_uniq"> | |
| 768 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 769 <change_format> | |
| 770 <when input="result['format']" value="sam" format="sam"/> | |
| 771 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 772 </change_format> | |
| 773 </data> | |
| 774 <data format="txt" name="concordant_transloc" label="${tool.name} on ${on_string} concordant_transloc.${result.format}" from_work_dir="gsnap_out.concordant_transloc"> | |
| 775 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 776 <change_format> | |
| 777 <when input="result['format']" value="sam" format="sam"/> | |
| 778 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 779 </change_format> | |
| 780 </data> | |
| 781 | |
| 782 <data format="txt" name="nomapping" label="${tool.name} on ${on_string} nomapping.${result.format}" from_work_dir="gsnap_out.nomapping"> | |
| 783 <filter>(results['split_output'] == 'yes' and results['fails_as_input'] == False)</filter> | |
| 784 <change_format> | |
| 785 <when input="result['format']" value="sam" format="sam"/> | |
| 786 <when input="result['format']" value="gsnap" format="gsnap"/> | |
| 787 </change_format> | |
| 788 </data> | |
| 789 | |
| 790 <data format="fastq" name="nomapping_fq" label="${tool.name} on ${on_string} nomapping.fq" from_work_dir="gsnap_out.nomapping.fq"> | |
| 791 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == False)</filter> | |
| 792 </data> | |
| 793 | |
| 794 <data format="fastq" name="nomapping_1_fq" label="${tool.name} on ${on_string} nomapping.1.fq" from_work_dir="gsnap_out.nomapping.1.fq"> | |
| 795 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 796 </data> | |
| 797 | |
| 798 <data format="fastq" name="nomapping_2_fq" label="${tool.name} on ${on_string} nomapping.2.fq" from_work_dir="gsnap_out.nomapping.2.fq"> | |
| 799 <filter>(results['split_output'] == 'yes' and seq['format'] == 'fastq' and seq['paired']['ispaired'] == True)</filter> | |
| 800 </data> | |
| 801 | |
| 802 <!-- Will problay need wrapper code to generate composite datatype for goby alignment | |
| 803 <data format="gobyalignment" name="goby_alignment" label="${tool.name} on ${on_string} uniq.${result.format}" from_work_dir="gsnap_out.nomapping"> | |
| 804 <filter>result['format'] == 'goby'</filter> | |
| 805 </data> | |
| 806 --> | |
| 807 | |
| 808 </outputs> | |
| 809 <tests> | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
810 </tests> |
| 0 | 811 |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
812 <help><![CDATA[ |
| 0 | 813 |
| 814 **What it does** | |
| 815 | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
816 GSNAP_ (Genomic Short-read Nucleotide Alignment Program) is a short read aligner which can align both single- and paired-end reads as short as 14nt and of arbitrarily long length. It can detect short- and long-distance splicing, including interchromosomal splicing, in individual reads, using probabilistic models or a database of known splice sites. Our program also permits SNP-tolerant alignment to a reference space of all possible combinations of major and minor alleles, and can align reads from bisulfite-treated DNA for the study of methylation state. It is developed by Thomas D. Wu of Genentech, Inc. |
| 0 | 817 Publication_ citation: Thomas D. Wu, Serban Nacu "Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010 Apr 1;26(7):873-81. Epub 2010 Feb 10. |
| 818 | |
| 819 .. _GSNAP: http://research-pub.gene.com/gmap/ | |
| 820 .. _Publication: http://bioinformatics.oupjournals.org/cgi/content/full/26/7/873 | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
821 https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2844994/?tool=pubmed |
| 0 | 822 |
| 823 ------ | |
| 824 | |
| 825 **Know what you are doing** | |
| 826 | |
| 827 .. class:: warningmark | |
| 828 | |
| 829 You will want to read the README_ | |
| 830 | |
| 831 .. _README: http://research-pub.gene.com/gmap/src/README | |
| 832 | |
| 833 ------ | |
| 834 | |
| 835 **Input formats** | |
| 836 | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
837 Input to GSNAP should be either in FASTQ or FASTA format. |
| 0 | 838 |
| 839 The FASTQ input may include quality scores, which will then be included in SAM | |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
840 output, if that output format is selected. |
| 0 | 841 |
| 842 For FASTA format, you should include one line per read (or end of a | |
| 843 paired-end read). The same FASTA file can have a mixture of | |
| 844 single-end and paired-end reads of varying lengths, if desired. | |
| 845 | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
846 *Single-end reads*: |
| 0 | 847 |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
848 Each FASTA entry should contain one short read per line, like this:: |
| 0 | 849 |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
850 >Header information |
|
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
851 AAAACATTCTCCTCCGCATAAGCCTGCGTCAGATTA |
| 0 | 852 |
| 853 Each short read can have a different length. However, the entire read | |
| 854 needs to be on a single line, and may not wrap around multiple lines. | |
| 855 If it extends to a second line, GSNAP will think that the read is | |
| 856 paired-end. | |
| 857 | |
| 858 | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
859 *Paired-end reads*: |
| 0 | 860 |
| 861 Each FASTA entry should contain two short reads, one per line, like | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
862 this:: |
| 0 | 863 |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
864 >Header information |
|
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
865 AAAACATTCTCCTCCGCATAAGCCTAGTAGATTA |
|
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
866 GGCGTAGGTAGAAGTAGAGGTTAAGGCGCGTCAG |
| 0 | 867 |
| 868 By default, the program assumes that the second end is in the reverse | |
| 869 complement direction compared with the first end. If they are in the | |
| 870 same direction, you may need to use the --circular-input (or -c) flag. | |
| 871 | |
| 872 ( The Galaxy tool: "FASTA Width formatter" can be used to reformat fasta files to have single line sequences. ) | |
| 873 | |
| 874 ------ | |
| 875 | |
| 876 **Output formats in GSNAP** | |
| 877 | |
| 878 SAM output format | |
| 879 | |
| 880 Default GSNAP format | |
|
5
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
881 |
|
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
882 See the README_ |
|
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
883 |
|
14561eb803a5
Uploaded v3.0.1b (still working on this prior to main Tool Shed release)
peterjc
parents:
3
diff
changeset
|
884 ]]></help> |
|
3
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
885 <citations> |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
886 <citation type="doi">10.1093/bioinformatics/btq057</citation> |
|
488e9d642566
GMAP wrappers v3.0.1 after linting and cleanup, still untested work-in-progress
peterjc
parents:
2
diff
changeset
|
887 </citations> |
| 0 | 888 </tool> |
| 889 |