Mercurial > repos > rnateam > segemehl

changeset 5:9ffdddb42700 draft

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 2e93c960ffa75613f0e8621a7e8e40dfeb81703c
author bgruening
date Tue, 25 Jul 2017 10:07:39 -0400
parents d3600fa8d3fd
children c6cef240817e
files segemehl.tar.gz segemehl.xml tool_dependencies.xml
diffstat 3 files changed, 136 insertions(+), 142 deletions(-) [+]
line wrap: on
line diff
Binary file segemehl.tar.gz has changed
--- a/segemehl.xml	Fri Dec 16 07:29:30 2016 -0500
+++ b/segemehl.xml	Tue Jul 25 10:07:39 2017 -0400
@@ -1,171 +1,171 @@
-<tool id="segemehl" name="segemehl" version="0.2.0">
-    <description>based short read aligner</description>
-    <requirements>
-        <requirement type="package" version="0.2.0">segemehl</requirement>
-    </requirements>
-    <stdio>
-        <regex match="Exit forced"
+<tool id="segemehl" name="segemehl" version="0.2.0.1">
+  <description>short read mapping with gaps</description>
+  <requirements>
+    <requirement type="package" version="0.2.0">segemehl</requirement>
+  </requirements>
+  <stdio>
+    <regex match="Exit forced"
            source="both"
            level="fatal"
            description="Execution halted." />
-    </stdio>
-    <command>
-<![CDATA[
+  </stdio>
+  <command>
+    <![CDATA[
         ## prepare segemehl index if no reference genome is supplied
         #if $refGenomeSource.genomeSource == "history":
             mkdir ./temp_index/ &&
-	    #set $temp_index = './temp_index/temp.idx'
-	    segemehl.x -x $temp_index -d $refGenomeSource.own_reference_genome &&
+        #set $temp_index = './temp_index/temp.idx'
+            segemehl.x -x $temp_index -d $refGenomeSource.own_reference_genome &&
         #else:
-            #set $temp_index = $refGenomeSource.index.fields.index_path
+        #set $temp_index = $refGenomeSource.index.fields.index_path
         #end if
-
-
+        
         ## execute segemehl
-        segemehl.x
-
+            segemehl.x
+        
         ## number of threads
-        -t "\${GALAXY_SLOTS:-12}"
-
+            -t "\${GALAXY_SLOTS:-12}"
+        
         #if $refGenomeSource.genomeSource == "history":
-	    -d $refGenomeSource.own_reference_genome
+            -d $refGenomeSource.own_reference_genome
         #else:
             -d ${refGenomeSource.index.fields.db_path}
         #end if
-
-        -i $temp_index
-
+        
+            -i $temp_index
+        
         ## check for single/pair-end
         #if str( $library.type ) == "single":
-            #set $query_list = list()
-            ## prepare inputs
-            #for $fastq in $library.input_query:
-                $query_list.append('%s' % $fastq )
-            #end for
+        #set $query_list = list()
+        ## prepare inputs
+        #for $fastq in $library.input_query:
+            $query_list.append('%s' % $fastq )
+        #end for
             -q "#echo ' '.join( $query_list )#"
         #else
-            ## prepare inputs
-
-            #set $mate1 = list()
-            #set $mate2 = list()
-            #for $mate_pair in $library.mate_list:
-                $mate1.append( str($mate_pair.first_strand_query) )
-                $mate2.append( str($mate_pair.second_strand_query) )
-            #end for
-
+        ## prepare inputs
+        
+        #set $mate1 = list()
+        #set $mate2 = list()
+        #for $mate_pair in $library.mate_list:
+            $mate1.append( str($mate_pair.first_strand_query) )
+            $mate2.append( str($mate_pair.second_strand_query) )
+        #end for
+        
             -q #echo ','.join($mate1)
             -p #echo ','.join($mate2)
-
+        
             -I $library.maxinsertsize
         #end if
-        -m $minsize
-        -A $accuracy
-        -H $hitstrategy
+            -m $minsize
+            -A $accuracy
+            -H $hitstrategy
         #if str( $prime5 ).strip():
             -P "$prime5"
         #end if
         #if str( $prime3 ).strip():
             -Q "$prime3"
         #end if
-        $polyA
-        $autoclip
-        $hardclip
-        $order
-	$splits
+            $polyA
+            $autoclip
+            $hardclip
+            $order
+            $splits
         #if $maxout:
             --maxout $maxout
         #end if
-        -s
-
-        --minsplicecover $minsplicecover
-        --minfragscore $minfragscore
-        --minfraglen $minfraglen
-        --splicescorescale $splicescorescale
+            -s
+        
+            --minsplicecover $minsplicecover
+            --minfragscore $minfragscore
+            --minfraglen $minfraglen
+            --splicescorescale $splicescorescale
+            
+            -o '$segemehl_out'
 
-        -o '$segemehl_out'
-]]>
-    </command>
-    <inputs>
-        <conditional name="refGenomeSource">
-          <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">
-            <option value="indexed">Use a built-in index</option>
-            <option value="history">Use one from the history</option>
-          </param>
-          <when value="indexed">
-            <param name="index" type="select" label="Select a reference genome" help="If your genome of interest is not listed, contact your Galaxy admin">
-              <options from_data_table="segemehl_indexes">
-                <column name="value" index="0"/>
-                <column name="dbkey" index="1"/>
-                <column name="name" index="2"/>
-                <column name="db_path" index="3"/>
-                <column name="index_path" index="4"/>
-                <filter type="sort_by" column="2"/>
-                <validator type="no_options" message="No indexes are available for the selected input dataset"/>
-              </options>
-            </param>
-          </when>  <!-- build-in -->
-          <when value="history">
-              <param name="own_reference_genome" type="data" format="fasta" label="Select the reference genome" />
-          </when>  <!-- history -->
-        </conditional>  <!-- refGenomeSource -->
+    ]]>
+  </command>
+  <inputs>
+    <conditional name="refGenomeSource">
+      <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">
+        <option value="indexed">Use a built-in index</option>
+        <option value="history">Use one from the history</option>
+      </param>
+      <when value="indexed">
+        <param name="index" type="select" label="Select a reference genome" help="If your genome of interest is not listed, contact your Galaxy admin">
+          <options from_data_table="segemehl_indexes">
+            <column name="value" index="0"/>
+            <column name="dbkey" index="1"/>
+            <column name="name" index="2"/>
+            <column name="db_path" index="3"/>
+            <column name="index_path" index="4"/>
+            <filter type="sort_by" column="2"/>
+            <validator type="no_options" message="No indexes are available for the selected input dataset"/>
+          </options>
+        </param>
+      </when>  <!-- build-in -->
+      <when value="history">
+        <param name="own_reference_genome" type="data" format="fasta" label="Select the reference genome" />
+      </when>  <!-- history -->
+    </conditional>  <!-- refGenomeSource -->
 
-        <conditional name="library">
-            <param name="type" type="select" label="Is this library paired-end?">
-                <option value="single">Single-end</option>
-                <option value="paired">Paired-end</option>
-            </param>
-            <when value="single">
-                <param name="input_query" type="data" multiple="True" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads in FASTQ/FASTA files" />
-            </when>
-            <when value="paired">
-                <!-- ToDo paired coolections -->
-                <repeat name="mate_list" title="Paired End Pairs" min="1">
-                    <param name="first_strand_query" type="data" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads from first strand" />
-                    <param name="second_strand_query" type="data" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads from second strand" />
-                </repeat>
-                <param name="maxinsertsize" type="integer" value="5000" label="Maximum size of the inserts (paired end)" help="default: 5000 (-I)" />
-            </when>
-        </conditional>
+    <conditional name="library">
+      <param name="type" type="select" label="Is this library paired-end?">
+        <option value="single">Single-end</option>
+        <option value="paired">Paired-end</option>
+      </param>
+      <when value="single">
+        <param name="input_query" type="data" multiple="True" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads in FASTQ/FASTA files" />
+      </when>
+      <when value="paired">
+        <!-- ToDo paired coolections -->
+        <repeat name="mate_list" title="Paired End Pairs" min="1">
+          <param name="first_strand_query" type="data" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads from first strand" />
+          <param name="second_strand_query" type="data" format="fastqsanger,fastqillumina,fastq,fasta" label="Reads from second strand" />
+        </repeat>
+        <param name="maxinsertsize" type="integer" value="5000" label="Maximum size of the inserts (paired end)" help="default: 5000 (-I)" />
+      </when>
+    </conditional>
 
-        <param name="minsplicecover" type="integer" value="80" label="Min coverage for spliced transcripts" help="(--minsplicecover)" />
-        <param name="minfragscore" type="integer" value="18" label="Min coverage for spliced transcripts" help="(--minfragscore)" />
-        <param name="minfraglen" type="integer" value="20" label="Min length of a spliced fragment" help="(--minfraglen)" />
-        <param name="splicescorescale" type="float" value="1.0" label="Report spliced alignment with score greater than this scale times the score"
-            help="Report only if this value x score is larger than next best spliced alignment (--splicescorescale)" />
+    <param name="minsplicecover" type="integer" value="80" label="Min coverage for spliced transcripts" help="(--minsplicecover)" />
+    <param name="minfragscore" type="integer" value="18" label="Min coverage for spliced transcripts" help="(--minfragscore)" />
+    <param name="minfraglen" type="integer" value="20" label="Min length of a spliced fragment" help="(--minfraglen)" />
+    <param name="splicescorescale" type="float" value="1.0" label="Report spliced alignment with score greater than this scale times the score"
+           help="Report only if this value x score is larger than next best spliced alignment (--splicescorescale)" />
 
-        <param name="minsize" type="integer" value="12" min="1" label="Minimum size of queries" help="(-m)" />
+    <param name="minsize" type="integer" value="12" min="1" label="Minimum size of queries" help="(-m)" />
 
-        <param name="maxout" type="integer" min="0" value="0" optional="True" 
-            label="Maximum number of alignments that will be reported" help="(--maxout)" />
-        <param name="accuracy" type="integer" value="85" min="1" max="100" label="Min percentage of matches per read in semi-global alignment" help="(-A)" />
+    <param name="maxout" type="integer" min="0" value="0" optional="True" 
+           label="Maximum number of alignments that will be reported" help="(--maxout)" />
+    <param name="accuracy" type="integer" value="85" min="1" max="100" label="Min percentage of matches per read in semi-global alignment" help="(-A)" />
 
-        <param name="hitstrategy" type="select" label="Hits to report?" help="(-H)">
-            <option value="1">report only best scoring hits</option>
-            <option value="0">report all scoring hits</option>
-        </param>
-        <param name="prime5" type="text" label="add 5' adapter" help="default: none (-Q)" />
-        <param name="prime3" type="text" label="add 3' adapter" help="default: none (-P)"/>
-        <param name="polyA" type="boolean" truevalue="--polyA" falsevalue="" checked="false" label="Clip polyA tail" help="(-T)"/>
-        <param name="autoclip" type="boolean" truevalue="--autoclip" falsevalue="" checked="false" label="Autoclip unknown 3prime adapter" help="(-Y)"/>
-        <param name="hardclip" type="boolean" truevalue="--hardclip" falsevalue="" checked="false" label="Enable hard clipping" help="(-C)"/>
-        <param name="order" type="boolean" truevalue="--order" falsevalue="" checked="false" label="Sorts the output by chromsome and position" help="(-O)"/>
-        <param name="splits" type="boolean" truevalue="--splits" falsevalue="" checked="false" label="Detect split/spliced reads" help="(--splits)"/>
-    </inputs>
-    <outputs>
-        <data format="sam" name="segemehl_out" label="Read alignments on ${on_string}"/>
-    </outputs>
-    <tests>
-      <test>
-	<param name="genomeSource" value="history" />
-        <param name="own_reference_genome" value="chr1.fa" />
-	<param name="library" value="single" />
-	<param name="input_query" value="test.fastq" />	
-	<param name="splits" value="true" />
-        <output name="segemehl_out" file="testmap.sam" lines_diff="2" />
-      </test>
-    </tests>
-    <help>
-<![CDATA[
+    <param name="hitstrategy" type="select" label="Hits to report?" help="(-H)">
+      <option value="1">report only best scoring hits</option>
+      <option value="0">report all scoring hits</option>
+    </param>
+    <param name="prime5" type="text" label="add 5' adapter" help="default: none (-Q)" />
+    <param name="prime3" type="text" label="add 3' adapter" help="default: none (-P)"/>
+    <param name="polyA" type="boolean" truevalue="--polyA" falsevalue="" checked="false" label="Clip polyA tail" help="(-T)"/>
+    <param name="autoclip" type="boolean" truevalue="--autoclip" falsevalue="" checked="false" label="Autoclip unknown 3prime adapter" help="(-Y)"/>
+    <param name="hardclip" type="boolean" truevalue="--hardclip" falsevalue="" checked="false" label="Enable hard clipping" help="(-C)"/>
+    <param name="order" type="boolean" truevalue="--order" falsevalue="" checked="false" label="Sorts the output by chromsome and position" help="(-O)"/>
+    <param name="splits" type="boolean" truevalue="--splits" falsevalue="" checked="false" label="Detect split/spliced reads" help="(--splits)"/>
+  </inputs>
+  <outputs>
+    <data format="sam" name="segemehl_out" label="Read alignments on ${on_string}"/>
+  </outputs>
+  <tests>
+    <test>
+      <param name="genomeSource" value="history" />
+      <param name="own_reference_genome" value="chr1.fa" />
+      <param name="library" value="single" />
+      <param name="input_query" value="test.fastq" />
+      <param name="splits" value="true" />
+      <output name="segemehl_out" file="testmap.sam" lines_diff="2" />
+    </test>
+  </tests>
+  <help>
+    <![CDATA[
 
 .. class:: infomark
 
@@ -181,9 +181,9 @@
 .. _Segemehl: http://www.bioinf.uni-leipzig.de/Software/segemehl/
 
 
-]]>
-    </help>
-    <citations>
-        <citation type="doi">10.1371/journal.pcbi.1000502</citation>
-    </citations>
+    ]]>
+  </help>
+  <citations>
+    <citation type="doi">10.1371/journal.pcbi.1000502</citation>
+  </citations>
 </tool>
--- a/tool_dependencies.xml	Fri Dec 16 07:29:30 2016 -0500
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,6 +0,0 @@
-<?xml version="1.0"?>
-<tool_dependency>
-    <package name="segemehl" version="0.2.0">
-        <repository changeset_revision="1cea83b71281" name="package_segemehl_0_2_0" owner="iuc" toolshed="https://testtoolshed.g2.bx.psu.edu" />
-    </package>
-</tool_dependency>