|
5
|
1 <tool id="aurora_htseq" name="aurora_htseq" version="2.0.0">
|
|
4
|
2 <description>Counting reads in features.</description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package" version="1.15.0.6-0">pandoc</requirement>
|
|
|
5 <requirement type="package" version="1.20.0">r-getopt</requirement>
|
|
|
6 <requirement type="package" version="1.6">r-rmarkdown</requirement>
|
|
|
7 <requirement type="package" version="0.9.1">htseq</requirement>
|
|
|
8 <requirement type="package" version="0.2">r-dt</requirement>
|
|
|
9 </requirements>
|
|
|
10 <stdio>
|
|
|
11 <regex match="XXX" source="stderr" level="warning"
|
|
|
12 description="Check the warnings_and_errors.txt file for more details."/>
|
|
|
13 </stdio>
|
|
5
|
14 <command><![CDATA[
|
|
|
15
|
|
|
16 ######### each aurora tool generates a html file and have an files path directory associated with it.
|
|
|
17 mkdir -p $report.files_path &&
|
|
0
|
18
|
|
5
|
19 ######### three important paths:
|
|
|
20 ######### 1. path to tool installation directory
|
|
|
21 ######### 2. path to report html
|
|
|
22 ######### 3. path to files_path directory associated with the report output.
|
|
|
23 export TOOL_INSTALL_DIR='${__tool_directory__}' &&
|
|
|
24 export REPORT='$report' &&
|
|
|
25 export REPORT_FILES_PATH='$report.files_path' &&
|
|
|
26
|
|
|
27 ############ create a hidden file to store r markdown rendering log
|
|
|
28 echo '------ Rmd rendering log -------\n' > $report.files_path/.r_rendering.log.txt &&
|
|
|
29
|
|
|
30 echo echo 'R markdown rendering started.' &&
|
|
|
31
|
|
|
32 Rscript '${__tool_directory__}/htseq_count_render.R'
|
|
|
33
|
|
0
|
34 -o $report
|
|
|
35 -d $report.files_path
|
|
5
|
36
|
|
2
|
37 -A '$alignment_files'
|
|
3
|
38 -B '$sample_names'
|
|
2
|
39 -G $gff
|
|
1
|
40 -f $format
|
|
|
41 -r $order
|
|
|
42 -S $stranded
|
|
|
43 -a $minaqual
|
|
|
44 -T $feature_type
|
|
|
45 -i $idattr
|
|
2
|
46 -m $mode
|
|
4
|
47 -c $count
|
|
5
|
48 -O $count_rdata &&
|
|
|
49
|
|
|
50 echo 'R markdown rendering finished.'
|
|
|
51
|
|
|
52 ]]></command>
|
|
4
|
53 <inputs>
|
|
|
54 <param type="data" name="alignment_files" label="Alignment file (SAM/BAM)" optional="False" format="sam,bam"
|
|
|
55 multiple="True"/>
|
|
|
56 <param type="text" name="sample_names" label="Sample names"
|
|
|
57 help="sample names for each input SAM/BAM file. Use comma to separate multiple names, for example, sample_1, sample_2, sample_3. The number of sample names should be the same as the number of input SAM/BAM files."
|
|
|
58 optional="False" size="5x25"/>
|
|
|
59 <param type="select" name="format" argument="-f" label="Input data format" optional="False" multiple="False">
|
|
|
60 <option value="sam" selected="false">SAM</option>
|
|
|
61 <option value="bam" selected="true">BAM</option>
|
|
|
62 </param>
|
|
|
63 <param type="data" name="gff" label="GFF file" optional="False" format="gff, gff3"/>
|
|
|
64 <param type="select" name="order" argument="-r" label="Order"
|
|
|
65 help="Use this option to indicate how the input data has been sorted. The default is name. "
|
|
|
66 optional="False" multiple="False">
|
|
|
67 <option value="name" selected="false">Name</option>
|
|
|
68 <option value="pos" selected="false">Position</option>
|
|
|
69 </param>
|
|
|
70 <param type="select" name="stranded" argument="-s" label="Stranded" optional="False" multiple="False">
|
|
|
71 <option value="yes" selected="true">Yes</option>
|
|
|
72 <option value="no" selected="false">No</option>
|
|
|
73 <option value="reverse" selected="false">Reverse</option>
|
|
|
74 </param>
|
|
|
75 <param type="integer" name="minaqual" argument="-a" label="Alignment quality "
|
|
|
76 help="Skip all reads with alignment quality lower than the given minimum value (default: 10)"
|
|
|
77 optional="False" value="10" min="0"/>
|
|
|
78 <param type="text" name="feature_type" argument="-t" label="Feature type"
|
|
|
79 help="Feature type (3rd column in GFF file) to be used, all features of other type are ignored (default, suitable for RNA-Seq analysis using an Ensembl GTF file: exon)"
|
|
|
80 optional="False" value="exon"/>
|
|
|
81 <param type="text" name="idattr" argument="-i" label="ID attribute"
|
|
|
82 help="GFF attribute to be used as feature ID. Several GFF lines with the same feature ID will be considered as parts of the same feature. The feature ID is used to identity the counts in the output table. The default, suitable for RNA-Seq analysis using an Ensembl GTF file, is gene_id."
|
|
|
83 optional="False" value="gene_id"/>
|
|
|
84 <param type="select" name="mode" argument="-m" label="Mode"
|
|
|
85 help="Mode to handle reads overlapping more than one feature. Possible values for <mode> are union, intersection-strict and intersection-nonempty (default: union)"
|
|
|
86 optional="False" multiple="False">
|
|
|
87 <option value="union" selected="true">Union</option>
|
|
|
88 <option value="intersection-strict" selected="false">Intersection (strict)</option>
|
|
|
89 <option value="intersection-nonempty" selected="false">Intersection (nonempty)</option>
|
|
|
90 </param>
|
|
|
91 </inputs>
|
|
|
92 <outputs>
|
|
5
|
93 <data name="report" format="html" label="${tool.name} report on ${on_string}" hidden="false"/>
|
|
4
|
94 <data name="count" format="txt" label="Aurora htseq-count on ${on_string}" hidden="false"/>
|
|
5
|
95 <data name="count_rdata" format="rdata" label="${tool.name} count data on ${on_string}"/>
|
|
4
|
96 </outputs>
|
|
|
97 <citations>
|
|
0
|
98 <citation type="bibtex"><![CDATA[
|
|
|
99 @article{allaire2016rmarkdown,
|
|
|
100 title={rmarkdown: Dynamic Documents for R, 2016},
|
|
|
101 author={Allaire, J and Cheng, Joe and Xie, Yihui and McPherson, Jonathan and Chang, Winston and Allen, Jeff
|
|
|
102 and Wickham, Hadley and Atkins, Aron and Hyndman, Rob},
|
|
|
103 journal={R package version 0.9},
|
|
|
104 volume={6},
|
|
|
105 year={2016}
|
|
|
106 }
|
|
4
|
107 ]]></citation>
|
|
|
108 <citation type="bibtex"><![CDATA[
|
|
0
|
109 @book{xie2015dynamic,
|
|
|
110 title={Dynamic Documents with R and knitr},
|
|
|
111 author={Xie, Yihui},
|
|
|
112 volume={29},
|
|
|
113 year={2015},
|
|
|
114 publisher={CRC Press}
|
|
|
115 }
|
|
4
|
116 ]]></citation>
|
|
|
117 <citation type="bibtex"><![CDATA[@article{anders2015htseq,
|
|
|
118 title={HTSeq—a Python framework to work with high-throughput sequencing data},
|
|
|
119 author={Anders, Simon and Pyl, Paul Theodor and Huber, Wolfgang},
|
|
|
120 journal={Bioinformatics},
|
|
|
121 volume={31},
|
|
|
122 number={2},
|
|
|
123 pages={166--169},
|
|
|
124 year={2015},
|
|
|
125 publisher={Oxford University Press}
|
|
|
126 }]]></citation>
|
|
|
127 </citations>
|
|
0
|
128 </tool>
|
