diff vmap_rnamap.xml @ 0:3ba5983012cf draft

Uploaded
author mzeidler
date Tue, 24 Sep 2013 10:19:40 -0400
parents
children a1d290099305
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/vmap_rnamap.xml	Tue Sep 24 10:19:40 2013 -0400
@@ -0,0 +1,186 @@
+<tool id="vmap_rnamap" name="Virana RNAMap">
+
+<description>Map input reads against a STAR index. </description>
+
+	<requirements>
+        <requirement type="package" version="2.3.0">STAR</requirement>
+        <requirement type="package" version="0.1.19">samtools</requirement>
+        <requirement type="package" version="1.0">virana_python</requirement>
+        <requirement type="package" version="1.7.1">numpy</requirement>
+        <requirement type="package" version="1.61">biopython</requirement>
+        <requirement type="package" version="0.5.4">htseq</requirement>
+    </requirements>
+
+    <command interpreter="python">vmap.py rnamap
+    
+    
+    
+    --reads $reads
+    
+    #if $paired.select_paired == "paired"
+        --reads $paired.paired_read
+    #end if
+    
+    #if $refseq
+        #set $groups = str($refseq).split(",")
+        #for $ref in $groups:
+            --virana_hit_filter $ref
+        #end for
+    #end if
+        
+
+        #if $bam
+            --bam $bam_output
+        #end if
+        
+        #if $sam
+            --sam $sam_output
+        #end if
+        
+        #if $output_type.hit
+            --virana_hits $hit_output
+            
+            $output_type.filter_complexity
+            
+            #if $output_type.max_mismatches
+                --max_mismatches $output_type.max_mismatches
+            #end if
+            
+            #if $output_type.max_relative_mismatches
+                --max_relative_mismatches $output_type.max_relative_mismatches
+            #end if
+            
+            #if $output_type.min_alignment_score
+                --min_alignment_score $output_type.min_alignment_score
+            #end if
+            
+            #if $output_type.min_continiously_matching
+                --min_continiously_matching $output_type.min_continiously_matching
+            #end if
+            
+            #if $output_type.min_mapping_score
+                --min_mapping_score $output_type.min_mapping_score
+            #end if
+       #end if
+       #if $qual
+            --qual $qual_output
+       #end if
+        
+       #if $tax
+            --taxonomy $tax_output
+       #end if
+        #if $chimeric
+            --chimeric_mappings $chimeric_output
+        #end if
+        #if $unmapped
+            --unmapped_end_1 $unmapped1_output
+            #if $paired.select_paired == "paired"
+                --unmapped_end_2 $unmapped2_output
+            #end if
+        #end if
+    
+    #if $sample_id:
+        --sample_id $sample_id
+    #end if 
+        
+    #if $splice_junctions:
+        --splice_junctions $splice_junctions
+    #end if
+    
+    --index_dir $index.extra_files_path
+    
+    $sensitive
+    
+    2&gt;&amp;1
+        
+    </command>
+    
+    <inputs>
+        <param name="reads" type="data" format="fastq" label="Read file"/>
+        
+        <conditional name="paired">
+            <param name="select_paired" type="select" label="Paired end reads?" >
+                <option value="single">single end</option>
+                <option value="paired">paired end</option>
+            </param>
+            
+           <when value="paired">
+                <param name="paired_read" type="data" format="fastq" label="Paired end read file" />
+           </when>
+        
+        </conditional>
+        <conditional name="output_type">
+            <param name="hit" type="boolean" label="Hit Output"/>
+                           
+            <when value="hit">
+                <param name="max_mismatches" type="integer" min="0" max="10000000"  optional="true" label="Maximum number of mismatches for saved hits"/>
+                <param name="max_relative_mismatches" type="float" min="0" max="1" optional="true" label="Maximum number of mismatches relative to read length for saved hits " />
+                <param name="min_continiously_matching" type="integer" min="0" max="10000000" optional="true" label="Minimum number of continious matches for saved hits"  />
+                <param name="min_mapping_score" type="integer" min="1" max="255"  optional="true" label="Mimimum mapping score for saved hit"  />        
+                <param name="filter_complexity" type="boolean" truevalue="--filter_complexity" falsevalue="" label="Discard low-complexity reads"/>
+                <param name="min_alignment_score" type="integer" min="0" max="255" optional="true" label="Mimimum alignment score for reported hits"/>
+
+            </when>
+        
+        </conditional>
+        
+        <param name="bam" type="boolean" label="BAM Output"/>
+        <param name="sam" type="boolean" label="SAM Output"/>
+        <param name="qual" type="boolean"  label="Quality Output"/>
+        <param name="tax" type="boolean"  label="Taxonomy Output"/>
+        <param name="chimeric" type="boolean" label="Chimeric Mapping Ouput" />
+        <param name="unmapped" type="boolean" label="Unmapped Reads Output" />
+        
+        <param name="refseq" type="select" label="Set RefSeq taxonomic groups" multiple="true" display="checkboxes" optional="true">
+            <option value="UniVec">UniVec</option>
+            <option value="rRNA">rRNA</option>
+            <option value="Homo_sapiens">Homo sapiens</option>
+            <option value="Fungi">Fungi</option>
+            <option value="Plasmids">Plasmids</option>
+            <option value="Protozoa">Protozoa</option>
+            <option value="Viruses">Viruses</option>
+            <option value="Homo_sapiens_cDNA">Homo sapiens cDNA</option>
+            <option value="Fungi_DRAFT">Fungi DRAFT</option>
+            <option value="Bacteria">Bacteria</option>
+            <option value="Bacteria_DRAFT">Bacteria DRAFT</option>
+        </param>
+        
+        <param name="splice_junctions" type="data" format="txt" optional="true" label="splice junction file"/>
+        
+        <param name="index" type="data" format="star_index" label="STAR Index files"/>
+        
+        <param name="sample_id" type="text" optional="true" label="String used to designate sample information within the hit file"/>
+        
+        <param name="sensitive" type="boolean" truevalue="--sensitive" falsevalue="" label="If given, mapping will process slower and more sensitive"/>
+
+    </inputs>
+    
+    
+    <outputs>
+        <data format="hit" name="hit_output" label="">
+            <filter>output_type['hit']</filter>
+        </data>
+        <data format="bam" name="bam_output" label="" >
+            <filter>bam</filter>
+        </data>
+        <data format="sam" name="sam_output" label="">
+            <filter>sam</filter>
+        </data>
+        <data format="txt" name="qual_output" label="" >
+            <filter>qual</filter>
+        </data>
+        <data format="txt" name="tax_output" label="" >
+            <filter>tax</filter>
+        </data>
+        <data format="sam" name="chimeric_output" label="" >
+            <filter>chimeric</filter>
+        </data>
+        <data format="fastq" name="unmapped1_output" label="" >
+            <filter>unmapped</filter>
+        </data>
+        <data format="fastq" name="unmapped2_output" label="" >
+            <filter>unmapped</filter>
+            <filter> paired['select_paired'] == 'paired'</filter>
+        </data>
+    </outputs>
+</tool>