|
0
|
1 <tool id="edge_pro" name="EDGE-pro" version="1.0.0">
|
|
|
2 <description>Gene expression in Prokaryotes</description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package" version="1.3.1">edge-pro</requirement>
|
|
|
5 </requirements>
|
|
|
6 <command interpreter="python">
|
|
|
7 edge_pro.py
|
|
|
8 \${EDGE_PRO_SITE_OPTIONS:---num-threads=4}
|
|
|
9
|
|
|
10 ## Mandatory input parameters
|
|
|
11 --genome $genome
|
|
|
12 --ptt $ptt
|
|
|
13 --rnt $rnt
|
|
|
14 ## First input file always required
|
|
|
15 --input1 $singlePaired.input1
|
|
|
16 ## Second input only if input is paired-end
|
|
|
17 #if $singlePaired.sPaired == "paired"
|
|
|
18 --input2=$singlePaired.input2
|
|
|
19 #end if
|
|
|
20
|
|
|
21 ## Optional input parameters
|
|
|
22 #if $params.settingsType == "full"
|
|
|
23 #if str($params.minInsertSize)
|
|
|
24 --minInsertSize=$params.minInsertSize
|
|
|
25 #end if
|
|
|
26 #if str($params.maxInsertSize)
|
|
|
27 --maxInsertSize=$params.maxInsertSize
|
|
|
28 #end if
|
|
|
29 #if str($params.window)
|
|
|
30 --window=$params.window
|
|
|
31 #end if
|
|
|
32 #if str($params.utrSize)
|
|
|
33 --utrSize=$params.utrSize
|
|
|
34 #end if
|
|
|
35 #if str($params.similarity)
|
|
|
36 --similarity=$params.similarity
|
|
|
37 #end if
|
|
|
38 #if str($params.readLength)
|
|
|
39 --readLength=$params.readLength
|
|
|
40 #end if
|
|
|
41 #if str($params.minCoverage)
|
|
|
42 --minCoverage=$params.minCoverage
|
|
|
43 #end if
|
|
|
44 #end if
|
|
|
45
|
|
|
46 ## Outputs
|
|
|
47 --out-aln $out_aln
|
|
|
48 --out-rpkm $out_rpkm
|
|
|
49 --out-log $out_log
|
|
|
50 </command>
|
|
|
51 <inputs>
|
|
|
52 <conditional name="singlePaired">
|
|
|
53 <param name="sPaired" type="select" label="Is this library mate-paired?">
|
|
|
54 <option value="single">Single-end</option>
|
|
|
55 <option value="paired">Paired-end</option>
|
|
|
56 </param>
|
|
|
57 <when value="single">
|
|
|
58 <param format="fastqsanger" name="input1" type="data" label="FASTQ file" help="FASTQ format with Sanger-scaled quality values (Galaxy fastqsanger datatype)"/>
|
|
|
59 </when>
|
|
|
60 <when value="paired">
|
|
|
61 <param format="fastqsanger" name="input1" type="data" label="Forward FASTQ file" help="FASTQ format with Sanger-scaled quality values (Galaxy fastqsanger datatype)" />
|
|
|
62 <param format="fastqsanger" name="input2" type="data" label="Reverse FASTQ file" help="FASTQ format with Sanger-scaled quality values (Galaxy fastqsanger datatype)" />
|
|
|
63 </when>
|
|
|
64 </conditional>
|
|
|
65
|
|
|
66 <param format="fasta" name="genome" type="data" label="Select the reference genome from your history (-g)" help="FASTA format" />
|
|
|
67 <param format="ptt" name="ptt" type="data" label="Coordinates of coding genes (PTT file)" help="PTT file with coordinates of coding genes (-p)" />
|
|
|
68 <param format="rnt" name="rnt" type="data" label="Coordinates of structural RNAs (RNT file)" help="RNT file with coordinates of structural RNA (-r)" />
|
|
|
69
|
|
|
70 <conditional name="params">
|
|
|
71 <param name="settingsType" type="select" label="Parameter settings" help="For most needs, use default settings. If you want full control use Full Parameter List">
|
|
|
72 <option value="preSet">Use Defaults</option>
|
|
|
73 <option value="full">Full parameter list</option>
|
|
|
74 </param>
|
|
|
75 <when value="preSet" />
|
|
|
76 <!-- Full/advanced params. -->
|
|
|
77 <when value="full">
|
|
|
78 <param name="minInsertSize" type="integer" optional="true" value="0" label="Minimun insert size for a read pair (-m)" help="For paired-end reads only" />
|
|
|
79 <param name="maxInsertSize" type="integer" optional="true" value="500" label="Maximun insert size for a read pair (-M)" help="For paired-end reads only" />
|
|
|
80 <param name="window" type="integer" optional="true" value="100" label="Window length for coverage distribution (-w)" help="Used to distribute the coverage between two overlapping genes. See help below for details" />
|
|
|
81 <param name="utrSize" type="integer" optional="true" value="40" label="Size of the untranslated region (-i)" help="Enter the size of the untranslated region between the initial transcription site and the start codon" />
|
|
|
82 <param name="similarity" type="float" optional="true" value="0.15" label="Percentage for similar coverage (-x)" help="Enter the percentage used to determine when two coverage values are considered similar. See help below for details" />
|
|
|
83 <param name="readLength" type="integer" optional="true" value="" label="Read length (-l)" help="If not specified, the first 1000 reads are used to approximate the read length" />
|
|
|
84 <param name="minCoverage" type="integer" optional="true" value="3" label="Minimum average coverage for expressed genes (-c)" help="Coverage less than specified is assumed to be noise and gene is considered to not be expressed" />
|
|
|
85 </when> <!-- full -->
|
|
|
86 </conditional>
|
|
|
87 </inputs>
|
|
|
88
|
|
|
89 <outputs>
|
|
|
90 <data format="sam" name="out_aln" label="${tool.name} on ${on_string}: alignment"/>
|
|
|
91 <data format="tabular" name="out_rpkm" label="${tool.name} on ${on_string}: rpkm"/>
|
|
|
92 <data format="txt" name="out_log" label="${tool.name} on ${on_string}: log"/>
|
|
|
93 </outputs>
|
|
|
94
|
|
|
95 <help>
|
|
|
96
|
|
|
97 **What it does**
|
|
|
98
|
|
|
99 `EDGE-pro`_, Estimated Degree of Gene Expression in PROkaryots is an efficient software system to estimate gene expression levels in prokaryotic genomes from RNA-seq data. EDGE-pro uses Bowtie2 for alignment and then estimates expression directly from the alignment results.
|
|
|
100 EDGE-pro includes routines to assign reads aligning to overlapping gene regions accurately. 15% or more of bacterial genes overlap other genes, making this a significant problem for bacterial RNA-seq, one that is generally ignored by programs designed for eukaryotic RNA-seq experiments.
|
|
|
101
|
|
|
102 **Input files:**
|
|
|
103
|
|
|
104 .. class:: infomark
|
|
|
105
|
|
|
106 Input files with gene coordinates in PTT and RNT format can be retrieved with the Get EDGE-pro Files tool available in Galaxy, or downloaded from the `NCBI ftp repository`_.
|
|
|
107 This tool accepts files in Sanger FASTQ format (Galaxy *fastqsanger* datatype). Use the FASTQ Groomer tool to prepare your files.
|
|
|
108
|
|
|
109 .. _NCBI ftp repository: ftp://ftp.ncbi.nlm.nih.gov/genomes/Bacteria/
|
|
|
110
|
|
|
111 .. class:: warningmark
|
|
|
112
|
|
|
113 All 3 types of files (FASTA reference genome, PTT and RNT) must have the same order of chromosomes/plasmids (e.g. if chr1 is before chr2 in genome.fasta file, then chr1 must be before chr2 in ptt and rnt files as well). If there is no PTT or RNT file for one of chromosomes/plasmids, place this chromosome/plasmid at the end of the file.
|
|
|
114
|
|
|
115 **License and citation**
|
|
|
116
|
|
|
117 This Galaxy tool is Copyright © 2012-2013 `CRS4 Srl.`_ and is released under the `MIT license`_.
|
|
|
118
|
|
|
119 .. _CRS4 Srl.: http://www.crs4.it/
|
|
|
120 .. _MIT license: http://opensource.org/licenses/MIT
|
|
|
121
|
|
|
122 If you use this tool in Galaxy, please cite |Cuccuru2013|_.
|
|
|
123
|
|
|
124 .. |Cuccuru2013| replace:: Cuccuru, G., Orsini, M., Pinna, A., Sbardellati, A., Soranzo, N., Travaglione, A., Uva, P., Zanetti, G., Fotia, G. (2013) Orione, a web-based framework for NGS analysis in microbiology. *Submitted*
|
|
|
125 .. _Cuccuru2013: http://orione.crs4.it/
|
|
|
126
|
|
|
127 This tool uses `EDGE-pro`_, which is licensed separately. Please cite |Magoc2013|_.
|
|
|
128
|
|
|
129 .. _EDGE-pro: http://ccb.jhu.edu/software/EDGE-pro/
|
|
|
130 .. |Magoc2013| replace:: Magoc, T., Wood, D., Salzberg, S. L. (2013) EDGE-pro: Estimated Degree of Gene Expression in Prokaryotic Genomes. *Evol. Bioinform.* 2013:9, 127-136
|
|
|
131 .. _Magoc2013: http://www.la-press.com/edge-pro-estimated-degree-of-gene-expression-in-prokaryotic-genomes-article-a3586
|
|
|
132 </help>
|
|
|
133 </tool>
|
