Mercurial > repos > iuc > dexseq

comparison dexseq_count.xml @ 0:276833129f97 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/dexseq commit eecb633ff51d61e8f94f580bb96053434029ab78-dirty
author iuc
date Thu, 24 Sep 2015 17:16:56 -0400
parents
children a26a82bed63f
comparison
equal deleted inserted replaced
-1:000000000000 0:276833129f97
1 <tool id="dexseq_count" name="DEXSeq-Count" version="1.0">
2 <description>Prepare and count exon abundancies from RNA-seq data</description>
3 <requirements>
4 <requirement type="package" version="3.1.0">R</requirement>
5 <requirement type="package" version="1.14.2">dexseq</requirement>
6 <requirement type="package" version="0.6.1">htseq</requirement>
7 <requirement type="package" version="0.7.6">pysam</requirement>
8 </requirements>
9 <command>
10 <![CDATA[
11 #if $mode.mode_select == "prepare":
12 python \$DEXSEQ_ROOT/python_scripts/dexseq_prepare_annotation.py
13 -r $mode.aggregate
14 $mode.gtffile
15 $flattened_gtf_out
16 #elif $mode.mode_select == "count":
17 python \$DEXSEQ_ROOT/python_scripts/dexseq_count.py
18 -f bam
19 -p $mode.paired
20 -s $mode.stranded
21 -a $mode.qual
22 -r $mode.order
23 $mode.flattened_gtf_in
24 $mode.bamfile
25 $counts_file
26 #end if
27 ]]>
28 </command>
29 <stdio>
30 <!-- Anything other than zero is an error -->
31 <exit_code range="1:" />
32 <exit_code range=":-1" />
33 <!-- In case the return code has not been set propery check stderr too -->
34 <regex match="Error:" />
35 <regex match="Exception:" />
36 </stdio>
37 <inputs>
38 <conditional name="mode">
39 <param name="mode_select" type="select" label="Mode of operation">
40 <option value="prepare">Prepare annotation</option>
41 <option value="count">Count reads</option>
42 </param>
43 <when value="prepare">
44 <param name="gtffile" type="data" format="gff" label="GTF file"/>
45 <param name="aggregate" type="boolean" checked="True" truevalue="yes" falsevalue="no"
46 label="Aggretare genes with exons?"
47 help="Indicates whether two or more genes sharing an exon should be merged into an 'aggregate gene'. If 'no', the exons that can not be assiged to a single gene are ignored."/>
48 </when>
49 <when value="count">
50 <param name="bamfile" type="data" format="bam" label="Input bam file"/>
51 <param name="flattened_gtf_in" type="data" format="gff" label="DEXSeq compatible GTF file" help="Created by prepare mode"/>
52 <param name="paired" type="boolean" checked="True" truevalue="yes" falsevalue="no" label="Is libray paired end?"/>
53 <param name="stranded" type="select" label="Is library strand specific?">
54 <option value="no">No</option>
55 <option value="yes">Yes</option>
56 <option value="reverse">Yes, but reverse</option>
57 </param>
58 <param name="qual" type="integer" value="10" label="Skip all reads with alignment quality lower than the given minimum value"/>
59 <param name="order" type="select" label="Sorting order of alignments" help="If you generated your alignments using tophat, they are by default position sorted. Ignored for single-end data">
60 <option value="pos">By position</option>
61 <option value="name">By name</option>
62 </param>
63 </when>
64 </conditional>
65 </inputs>
66
67 <outputs>
68 <data format="tabular" name="counts_file" label="DEXSeq count reads on ${on_string}">
69 <filter>(mode['mode_select'] == 'count')</filter>
70 </data>
71 <data format="gff" name="flattened_gtf_out" label="DEXSeq prepare annotation ${on_string}">
72 <filter>(mode['mode_select'] == 'prepare')</filter>
73 </data>
74 </outputs>
75
76 <tests>
77 <test>
78 <param name="mode_select" value="prepare" />
79 <param name="gtffile" ftype="gff" value="original.gtf"/>
80 <param name="aggregate" value="True"/>
81 <output name="flattened_gtf_out" file="flattened.gtf" ftype="gff"/>
82 </test>
83 </tests>
84
85 <help>
86 <![CDATA[
87 .. class:: infomark
88
89 **What it does**
90
91 The main goal of this tol is to count the number of reads/fragments per exon of each gene in RNA-seq sample. In addition it also prepares your annotation gtf file compatible for counting.
92
93
94 **Inputs**
95
96 Mode-preprare: Takes a normal gtf file as input. For example from ensembl database.
97 Mode-count: Inputs are flattened gtf file and BAM file. The flattened gtf file can be generated from 'prepare' mode of this tool.
98
99 **Output**
100
101 Mode-preprare: Flattened gtf file that contains only exons with corresponding gene ids from given gtf file. Sometimes two or more genes sharing an exon will be merged into an 'aggregate gene' if the aggregate option was used.
102 Mode-count: Two column tab-delimeted file with exon ids and their read counts.
103
104 .. _DEXSeq: http://master.bioconductor.org/packages/release/bioc/html/DEXSeq.html
105
106 ]]>
107 </help>
108 <citations>
109 <citation type="doi">10.1101/gr.133744.111</citation>
110 </citations>
111 </tool>