|
1
|
1 <tool id="monorail" name="Run the Monorail RNA-seq analysis pipeline" version="0.1.0">
|
|
15
|
2 <!-- much of this was based on https://github.com/galaxyproject/tools-iuc/blob/master/tools/rgrnastar/rg_rnaStar.xml -->
|
|
1
|
3 <requirements>
|
|
|
4 <requirement type="package" version="1.9">samtools</requirement>
|
|
17
|
5 <requirement type="package" version="2.7.3a">star</requirement>
|
|
|
6 <requirement type="package" version="0.4.0">bamcount</requirement>
|
|
16
|
7 <requirement type="package" version="5.4.0">snakemake-minimal</requirement>
|
|
|
8 <requirement type="package" version="1.3.3">zstd</requirement>
|
|
|
9 <requirement type="package" version="1.3">seqtk</requirement>
|
|
|
10 <!--
|
|
|
11 <requirement type="package" version="2.1.0">hisat2</requirement>
|
|
|
12 <requirement type="package" version="0.8">fastq-tools</requirement>
|
|
|
13 <requirement type="package" version="0.12.0">salmon</requirement>
|
|
|
14 <requirement type="package" version="0.5.0">regtools</requirement>
|
|
|
15 <requirement type="package" version="2.9.1">sra-tools</requirement>
|
|
|
16 -->
|
|
1
|
17 </requirements>
|
|
|
18 <!-- /bin/bash -x monorail.slim.sh ../ath10 4 10 ../ath10/gtf/exons.bed ./tmp2 ../fastqs/SRR8505407_1_100.fastq.gz ../fastqs/SRR8505407_2_100.fastq.gz -->
|
|
8
|
19 <command detect_errors="aggressive"><![CDATA[
|
|
9
|
20 snakemake --snakefile "$__tool_directory__/Snakefile"
|
|
8
|
21 -j \${GALAXY_SLOTS:-4}
|
|
|
22 --config
|
|
1
|
23 #if str($singlePaired.sPaired) == "paired_collection"
|
|
8
|
24 inputs="$singlePaired.input.forward,$singlePaired.input.reverse"
|
|
15
|
25 #if $singlePaired.input.forward.is_of_type("fastq.gz"):
|
|
|
26 compressed=1
|
|
|
27 #end if
|
|
1
|
28 #else
|
|
8
|
29 inputs="$singlePaired.input1"
|
|
1
|
30 #if str($singlePaired.sPaired) == "paired"
|
|
8
|
31 inputs="$singlePaired.input1,$singlePaired.input2"
|
|
1
|
32 #end if
|
|
15
|
33 #if $singlePaired.input1.is_of_type("fastq.gz"):
|
|
|
34 compressed=1
|
|
|
35 #end if
|
|
1
|
36 #end if
|
|
17
|
37 ref="$refGenomeSource.genomeDir.fields.path"
|
|
|
38 exon_bed="$refGenomeSource.genomeDir.fields.exons_path"
|
|
12
|
39 output="."
|
|
8
|
40 temp="./tmp"
|
|
1
|
41 ]]></command>
|
|
|
42 <inputs>
|
|
|
43 <!-- FASTQ input(s) and options specifically for paired-end data. -->
|
|
|
44 <conditional name="singlePaired">
|
|
|
45 <param name="sPaired" type="select" label="Single-end or paired-end reads">
|
|
|
46 <option value="single" selected="true">Single-end</option>
|
|
|
47 <option value="paired">Paired-end (as individual datasets)</option>
|
|
|
48 <option value="paired_collection">Paired-end (as collection)</option>
|
|
|
49 </param>
|
|
|
50 <when value="single">
|
|
|
51 <param format="fastq,fastq.gz" name="input1" type="data" label="RNA-Seq FASTQ file"/>
|
|
|
52 </when>
|
|
|
53 <when value="paired">
|
|
|
54 <param format="fastq,fastq.gz" name="input1" type="data" label="RNA-Seq FASTQ file, forward reads"/>
|
|
|
55 <param format="fastq,fastq.gz" name="input2" type="data" label="RNA-Seq FASTQ file, reverse reads"/>
|
|
|
56 </when>
|
|
|
57 <when value="paired_collection">
|
|
|
58 <param format="fastq,fastq.gz" name="input" type="data_collection" collection_type="paired" label="RNA-Seq FASTQ paired reads"/>
|
|
|
59 </when>
|
|
|
60 </conditional>
|
|
|
61
|
|
|
62 <!--
|
|
|
63 <param name="refGenomeSource" type="select" label="Select reference genome" help="If your genome of interest is not listed, contact the Galaxy team">
|
|
|
64 <options from_data_table="rnastar_index2">
|
|
|
65 <filter type="static_value" column="4" value="0"/>
|
|
|
66 <filter type="sort_by" column="2" />
|
|
|
67 <validator type="no_options" message="No indexes are available for the selected input dataset"/>
|
|
|
68 </options>
|
|
|
69 </param>
|
|
|
70 -->
|
|
|
71
|
|
|
72 <!-- Genome source. -->
|
|
|
73 <conditional name="refGenomeSource">
|
|
17
|
74 <param name="genomeDir" argument="--genomeDir" type="select" label="Select reference genome index set" help="If your genome of interest is not listed, contact the Galaxy team">
|
|
|
75 <options from_data_table="monorail_index">
|
|
1
|
76 <validator type="no_options" message="No indexes are available for the selected input dataset"/>
|
|
|
77 </options>
|
|
|
78 </param>
|
|
|
79 </conditional>
|
|
|
80 </inputs>
|
|
8
|
81 <outputs>
|
|
|
82 <data format="txt" name="auc" label="${tool.name} on ${on_string}: AUC" from_work_dir="bamcount_auc.tsv"/>
|
|
|
83 <!--
|
|
6
|
84 <data format="txt" name="output_log" label="${tool.name} on ${on_string}: log" from_work_dir="Log.final.out"/>
|
|
7
|
85 <data format="txt" name="auc" label="${tool.name} on ${on_string}: AUC" from_work_dir="bc.auc.tsv"/>
|
|
|
86 <data format="txt" name="bc_log" label="${tool.name} on ${on_string}: bamcount log" from_work_dir="bc.log"/>
|
|
6
|
87 <data format="interval" name="chimeric_junctions" label="${tool.name} on ${on_string}: chimeric junctions" from_work_dir="Chimeric.out.junction"/>
|
|
|
88 <data format="interval" name="splice_junctions" label="${tool.name} on ${on_string}: splice junctions.bed" from_work_dir="SJ.out.tab"/>
|
|
|
89 <data name="mapped_reads" format="bam" label="${tool.name} on ${on_string}: mapped.bam" from_work_dir="sorted.bam"/>
|
|
|
90 <data name="mapped_reads_index" format="bai" label="${tool.name} on ${on_string}: mapped.bam.bai" from_work_dir="sorted.bam.bai"/>
|
|
|
91 -->
|
|
1
|
92 </outputs>
|
|
|
93 <help>
|
|
|
94 Run the Monorail RNA-seq analysis pipeline
|
|
|
95 </help>
|
|
|
96 </tool>
|
