18
|
1 <tool id="RSEM_abundance_estimation" name="RSEM_abundance_estimation" version="0.0.1">
|
|
2
|
|
3 <description>run RSEM to estimate transcript abundances</description>
|
|
4 <requirements>
|
|
5 <requirement type="package">trinity</requirement>
|
|
6 </requirements>
|
|
7 <command interpreter="python">
|
|
8
|
|
9 trinityToolWrapper.py util/align_and_estimate_abundance.pl --transcripts $transcripts --est_method RSEM --aln_method bowtie --trinity_mode --prep_reference
|
|
10
|
|
11 ## Inputs.
|
|
12 #if str($inputs.paired_or_single) == "paired":
|
|
13 --left $inputs.left_input --right $inputs.right_input
|
|
14 #if $inputs.left_input.ext == 'fa':
|
|
15 --seqType fa
|
|
16 #else:
|
|
17 --seqType fq
|
|
18 #end if
|
|
19 #if str($inputs.library_type) != "None":
|
|
20 --SS_lib_type $inputs.library_type
|
|
21 #end if
|
|
22
|
|
23 #else:
|
|
24 --single $inputs.input
|
|
25 #if str($inputs.input.ext) == 'fa':
|
|
26 --seqType fa
|
|
27 #else:
|
|
28 --seqType fq
|
|
29 #end if
|
|
30 #if str($inputs.library_type) != "None":
|
|
31 --SS_lib_type $inputs.library_type
|
|
32 #end if
|
|
33 #end if
|
|
34
|
|
35
|
|
36 </command>
|
|
37 <inputs>
|
|
38 <param format="fasta" name="transcripts" type="data" label="transcripts_fasta" help="Fasta sequences for which reads are aligned." />
|
|
39
|
|
40 <conditional name="inputs">
|
|
41 <param name="paired_or_single" type="select" label="Paired or Single-end data?">
|
|
42 <option value="paired">Paired</option>
|
|
43 <option value="single">Single</option>
|
|
44 </param>
|
|
45 <when value="paired">
|
|
46 <param format="fasta,fastq" name="left_input" type="data" label="Left/Forward strand reads" help=""/>
|
|
47 <param format="fasta,fastq" name="right_input" type="data" label="Right/Reverse strand reads" help=""/>
|
|
48 <param name="library_type" type="select" label="Strand-specific Library Type">
|
|
49 <option value="None">None</option>
|
|
50 <option value="FR">FR</option>
|
|
51 <option value="RF">RF</option>
|
|
52 </param>
|
|
53
|
|
54 </when>
|
|
55 <when value="single">
|
|
56 <param format="fasta,fastq" name="input" type="data" label="Single-end reads" help=""/>
|
|
57 <param name="library_type" type="select" label="Strand-specific Library Type">
|
|
58 <option value="None">None</option>
|
|
59 <option value="F">F</option>
|
|
60 <option value="R">R</option>
|
|
61 </param>
|
|
62
|
|
63 </when>
|
|
64 </conditional>
|
|
65
|
|
66
|
|
67 </inputs>
|
|
68 <outputs>
|
|
69 <data format="text" name="transcript_counts" label="${tool.name} on ${on_string}: Isoform Counts" from_work_dir="RSEM.isoforms.results"/>
|
|
70 <data format="text" name="gene_counts" label="${tool.name} on ${on_string}: Gene counts" from_work_dir="RSEM.genes.results"/>
|
|
71
|
|
72
|
|
73 </outputs>
|
|
74 <tests>
|
|
75
|
|
76
|
|
77 <test>
|
|
78 <param name="target" value="trinity/Trinity.fasta" />
|
|
79 <param name="aligner" value="bowtie" />
|
|
80 <param name="paired_or_single" value="single" />
|
|
81 <param name="library_type" value="None" />
|
|
82 <param name="input" value="trinity/reads.left.fq" />
|
|
83 </test>
|
|
84
|
|
85
|
|
86 </tests>
|
|
87 <help>
|
|
88 .. _Trinity: http://trinityrnaseq.sourceforge.net
|
|
89 </help>
|
|
90 </tool>
|