0
+ − 1 <tool id='id_align_and_estimate_abundance' name='align_and_estimate_abundance' version='1.0'>
+ − 2
+ − 3 <description>align_and_estimate_abundance</description>
+ − 4
+ − 5 <requirements>
+ − 6 <requirement type="package" version="2013_08_14">trinityrnaseq</requirement>
+ − 7 <requirement type="package" version="1.1.17">rsem</requirement>
+ − 8 <requirement type="package" version="1.5.1">eXpress</requirement>
+ − 9 <requirement type="package" version="2.1.0">bowtie2</requirement>
+ − 10 <requirement type="package" version="0.12.7">bowtie</requirement>
+ − 11 <requirement type="package" version="0.1.18">samtools</requirement>
+ − 12 <requirement type="set_environment">TRINITY_HOME</requirement>
+ − 13
+ − 14 </requirements>
+ − 15 <command >perl \$TRINITY_HOME/util/align_and_estimate_abundance.pl
+ − 16 --transcripts $transcripts
+ − 17 --est_method $est_method
+ − 18 --aln_method $aln_method
+ − 19 --prep_reference
+ − 20
+ − 21 ## Inputs.
+ − 22 #if str($reads.paired_or_single) == "paired":
+ − 23
+ − 24 --left $reads.left_input --right $reads.right_input
+ − 25
+ − 26 #if $reads.left_input.ext == 'fa':
+ − 27 --seqType fa
+ − 28 #else:
+ − 29 --seqType fq
+ − 30 #end if
+ − 31
+ − 32 ## Additional parameters.
+ − 33 #if str($reads.optional.use_options) == "yes":
+ − 34
+ − 35 #if str($reads.optional.library_type) != "None":
+ − 36 --SS_lib_type $reads.optional.library_type
+ − 37 #end if
+ − 38
+ − 39
+ − 40
+ − 41 #end if
+ − 42
+ − 43 #else:
+ − 44 --single $reads.input
+ − 45
+ − 46 #if str($reads.input.ext) == 'fa':
+ − 47 --seqType fa
+ − 48 #else:
+ − 49 --seqType fq
+ − 50 #end if
+ − 51
+ − 52 ## Additional parameters.
+ − 53 #if str($reads.optional.use_option) == "yes":
+ − 54
+ − 55 #if str($reads.additional_params.library_type) != "None":
+ − 56 --SS_lib_type $reads.optional.library_type
+ − 57 #end if
+ − 58
+ − 59 #end if
+ − 60 #end if
+ − 61
+ − 62 ## direct to output
+ − 63 > $align_and_estimate_abundance_log 2>&1
+ − 64
+ − 65 </command>
+ − 66
+ − 67 <inputs>
+ − 68
+ − 69 <!-- Required -->
+ − 70 <param name="transcripts" type="data" format="fasta" label="Transcripts [ in fasta format ]" help="--transcripts < string > ; transcript fasta file"/>
+ − 71
+ − 72 <param name="est_method" type="select" display="radio" label="Choose abundance estimation method" help="">
+ − 73 <option value="RSEM" selected="True">RSEM</option>
+ − 74 <option value="eXpress">eXpress</option>
+ − 75 </param>
+ − 76
+ − 77 <param name="aln_method" type="select" display="radio" label="Choose alignment method" help="">
+ − 78 <option value="bowtie" selected="True">bowtie</option>
+ − 79 <option value="bowtie2">bowtie2</option>
+ − 80 </param>
+ − 81
+ − 82 <conditional name="reads">
+ − 83 <param name="paired_or_single" type="select" label="Paired or Single-end reads?">
+ − 84 <option value="paired">Paired</option>
+ − 85 <option value="single">Single</option>
+ − 86 </param>
+ − 87
+ − 88
+ − 89 <when value="paired">
+ − 90 <param format="fastq" name="left_input" type="data" label="Left/Forward strand reads" help=""/>
+ − 91 <param format="fastq" name="right_input" type="data" label="Right/Reverse strand reads" help=""/>
+ − 92
+ − 93 <param name="seqtype" type="select" label="--seqType" help="">
+ − 94 <option value='fq' selected='true'>fastq</option>
+ − 95 <option value='fa'>fasta</option>
+ − 96 </param>
+ − 97
+ − 98 <conditional name='optional'>
+ − 99 <param name='use_options' type='select' label='Use Optional parameters?'>
+ − 100 <option value='no'>NO</option>
+ − 101 <option value='yes'>YES</option>
+ − 102 </param>
+ − 103 <when value='no'/>
+ − 104 <when value='yes'>
+ − 105 <param name="library_type" type="select" label="Strand-specific Library Type">
+ − 106 <option value="None">None</option>
+ − 107 <option value="FR">FR</option>
+ − 108 <option value="RF">RF</option>
+ − 109 </param>
+ − 110
+ − 111
+ − 112
+ − 113 </when>
+ − 114 </conditional>
+ − 115 </when>
+ − 116
+ − 117 <when value="single">
+ − 118 <param format="fastq" name="input" type="data" label="Single-end reads" help=""/>
+ − 119 <param name="seqtype" type="select" label="--seqType" help="">
+ − 120 <option value='fq' selected='true'>fastq</option>
+ − 121 <option value='fa'>fasta</option>
+ − 122 </param>
+ − 123
+ − 124 <conditional name='optional'>
+ − 125 <param name='use_options' type='select' label='Use Optional parameters?'>
+ − 126 <option value='no'>NO</option>
+ − 127 <option value='yes'>YES</option>
+ − 128 </param>
+ − 129 <when value='no'/>
+ − 130 <when value='yes'>
+ − 131 <param name="library_type" type="select" label="Strand-specific Library Type">
+ − 132 <option value="None">None</option>
+ − 133 <option value="F">F</option>
+ − 134 <option value="R">R</option>
+ − 135 </param>
+ − 136 </when>
+ − 137 </conditional>
+ − 138 </when>
+ − 139 </conditional>
+ − 140
+ − 141 </inputs>
+ − 142
+ − 143 <outputs>
+ − 144 <data format="txt" name="align_and_estimate_abundance_log" label="${tool.name} on ${on_string}: log" />
+ − 145 <data format="bam" name="bowtie_bam" label="${tool.name} on ${on_string}: bam" from_work_dir="bowtie.bam" />
+ − 146 <data format="tabular" name="rsem_isoforms" label="${tool.name} on ${on_string}: isoforms" from_work_dir="RSEM.isoforms.results"/>
+ − 147 <data format="tabular" name="rsem_genes" label="${tool.name} on ${on_string}: genes" from_work_dir="RSEM.genes.results"/>
+ − 148 </outputs>
+ − 149 </tool>
+ − 150