diff prepare.xml @ 0:e01de823e919 draft default tip

Uploaded
author biotechcoder
date Fri, 01 May 2015 05:41:51 -0400
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/prepare.xml	Fri May 01 05:41:51 2015 -0400
@@ -0,0 +1,174 @@
+<tool id="riboseqr_prepare_input" name="Prepare riboSeqR input" version="0.4.0">
+    <description>
+        (Step 1) Prepare alignment file (SAM format, Ribo-Seq or RNA-Seq alignments)
+        for riboSeqR analysis.
+    </description>
+    <requirements>
+        <requirement type="R-module">riboSeqR</requirement>
+    </requirements>
+    <command interpreter="python">riboseqr/prepare.py
+        #set $files_ribo = []
+        #set $files_rna = []
+        #set $replicates = []
+
+        #for $i, $s in enumerate($rnaseq.datasets)
+            #silent $replicates.append($s.replicate_name.value)
+            #silent $files_ribo.append($s.ribo_files.file_name)
+            #try:
+                #silent $s.rna_files
+            #except:
+                #pass
+            #else:
+                #silent $files_rna.append($s.rna_files.file_name)
+            #end try
+        #end for
+
+        #set $ribofiles = ", ".join($files_ribo)
+        #set $rnafiles = ", ".join($files_rna)
+        #set $replicate_names = ", ".join($replicates)
+
+        --ribo_files "$ribofiles"
+        --rna_files "$rnafiles"
+        --replicate_names "$replicate_names"
+        --seqnames "$seqnames"
+        --rdata_save "$rdata_save"
+        --sam_format
+        --html_file "$html_file"
+        --output_path "$html_file.files_path"
+    </command>
+    <inputs>
+        <conditional name="rnaseq">
+            <param name="have-rnaseq" type="boolean" truevalue="checked"
+                   falsevalue="unchecked" label="Have RNA-Seq data?"/>
+            <when value="checked">
+                <repeat name="datasets" title="Ribo and RNA-Seq alignment"
+                        default="1">
+                    <param type="text" name="replicate_name" size="20"
+                           label="Provide a common name for each replicate group"
+                           value=""/>
+                    <param format="sam" name="ribo_files" type="data"
+                           label="Ribo-Seq alignment file (SAM)"
+                           multiple="false"/>
+                    <param format="sam" name="rna_files" type="data"
+                           label="RNA-Seq alignment file (SAM)"
+                           multiple="false"/>
+                </repeat>
+            </when>
+            <when value="unchecked">
+                <repeat name="datasets" title="Ribo-Seq alignment" default="1">
+                    <param type="text" name="replicate_name" size="20"
+                           label="Provide a common name for each replicate group"
+                           value=""/>
+                    <param format="sam" name="ribo_files" type="data"
+                           label="Ribo-Seq alignment file (SAM)"
+                           multiple="false">
+                    </param>
+                </repeat>
+            </when>
+        </conditional>
+        <param name="seqnames" type="text" area="True" label="Transcript (seqname)
+             names to be read" size="3x60"
+               help="[Optional] - only the mapped footprints for these transcripts will be
+             used in further riboSeqR analysis. Multiple transcript names should
+             be separated by commas. There should be no quotes around transcript names.">
+            <sanitizer>
+                <valid>
+                    <add value="|"/>
+                </valid>
+            </sanitizer>
+        </param>
+    </inputs>
+    <outputs>
+        <data format="rda" name="rdata_save"
+              label="Prepare riboSeqR input (R data file)"/>
+        <data format="html" name="html_file"
+              label="Prepare riboSeqR input (HTML report)"/>
+    </outputs>
+    <tests>
+        <test>
+            <param name="ribo_files" value="12hRPF.sam" />
+            <output name="html_file" file="Prepare_riboSeqR_input_(HTML_report).html">
+                <extra_files type="file" name="RiboSeq file 1" value="RiboSeq file 1"/>
+            </output>
+        </test>
+    </tests>
+    <help>
+About
+-----
+This tool can be used to prepare input data for riboSeqR from SAM
+format alignments of Ribo or RNA-Seq data to a reference transcriptome. You can
+do this alignment manually using bowtie or the
+"Transcriptome Mapping" -> "Align to transcriptome using Bowtie"
+tool on RiboGalaxy_
+
+The required input format for riboSeqR_ is mentioned in the
+**"Getting Data"** section of the documentation_.
+
+riboSeqR version: **1.0.4**.
+
+How to use?
+-----------
+Inputs
+......
+Select SAM format Ribo-Seq alignment files in the input section.
+
+If you have RNA-Seq data, these can be included if the **"Have RNA-Seq data"**
+option is checked.
+
+Outputs
+.......
+The following files will be generated on completion:
+
+#. Prepare riboSeqR input (HTML report) - A HTML file with links to all other
+   output files.
+
+   * Generated riboSeqR format input files of RiboSeq and RNASeq data(if provided).
+     These files are plain text and lines have the following information -
+     strand, transcript name, alignment position, sequence.
+
+     Please note the alignments are made **0-indexed**.
+
+   * R script used in this session.
+
+#. Prepare riboSeqR input (R data file) - used as input for the next step -
+   *Triplet Periodicity*.
+
+How are the SAM alignments processed?
+.....................................
+#. Lines starting with ``@`` are ignored.
+#. Lines having a ``FLAG=0`` are considered as successful alignments. These are
+   considered for the next step.
+#. Alignment start is located on ``column 4``. These are decremented by
+   1 as SAM alignments are 1-indexed.
+#. riboSeqR input file is written with the strand (``+``), transcript name,
+   alignment start and the aligned sequence.
+
+riboSeqR functions used
+.......................
+readRibodata.
+
+For detailed description of the functions and the options used, please consult
+the riboSeqR documentation.
+
+Links
+.....
+* `Bioconductor package information on riboSeqR`__
+* riboSeqR - `Reference manual`_
+* riboSeqR - `Introduction and workflow example`_
+
+.. _RiboGalaxy: http://ribogalaxy.ucc.ie
+
+.. _riboSeqR: http://bioconductor.org/packages/3.0/bioc/html/riboSeqR.html
+
+__ riboSeqR_
+
+.. _manual: http://bioconductor.org/packages/3.0/bioc/manuals/riboSeqR/man/riboSeqR.pdf
+
+.. _`Reference manual`: manual_
+
+.. _documentation: http://bioconductor.org/packages/3.0/bioc/vignettes/riboSeqR/inst/doc/riboSeqR.pdf
+
+.. _`Introduction and workflow example`: documentation_
+
+    </help>
+</tool>