diff dunovo.xml @ 19:675a8370675b draft

planemo upload for repository https://github.com/galaxyproject/dunovo commit b'4303231da9e48b2719b4429a29b72421d24310f4\n'-dirty
author nick
date Thu, 02 Feb 2017 19:14:13 -0500
parents
children 53cf867b305b
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/dunovo.xml	Thu Feb 02 19:14:13 2017 -0500
@@ -0,0 +1,67 @@
+<?xml version="1.0"?>
+<tool id="duplex" name="Du Novo: Make consensus reads" version="0.5">
+  <description>from duplex sequencing alignments</description>
+  <requirements>
+    <requirement type="package" version="0.5">duplex</requirement>
+    <requirement type="set_environment">DUPLEX_DIR</requirement>
+    <!-- TODO: require Python 2.7 -->
+  </requirements>
+  <command detect_errors="exit_code"><![CDATA[
+    python "\$DUPLEX_DIR/dunovo.py" -r $min_reads -q $qual_thres -F $qual_format '$input'
+    #if $keep_sscs:
+      --sscs-file sscs.fa
+    #end if
+    > duplex.fa
+    && python "\$DUPLEX_DIR/utils/outconv.py" duplex.fa -1 '$dcs1' -2 '$dcs2'
+    #if $keep_sscs:
+      && python "\$DUPLEX_DIR/utils/outconv.py" sscs.fa -1 '$sscs1' -2 '$sscs2'
+    #end if
+  ]]>
+  </command>
+  <inputs>
+    <param name="input" type="data" format="tabular" label="Aligned input reads" />
+    <param name="min_reads" type="integer" value="3" min="1" label="Minimum reads per family" help="Single-strand families with fewer than this many reads will be skipped."/>
+    <param name="qual_thres" type="integer" value="25" min="1" label="Minimum base quality" help="Bases with a PHRED score less than this will not be counted in the consensus making."/>
+    <param name="qual_format" type="select" label="FASTQ format" help="Solexa should also work for Illumina 1.3+ and 1.5+, and Sanger should work for Illumina 1.8+">
+      <option value="sanger" selected="true">Sanger (PHRED 0 = &quot;!&quot;)</option>
+      <option value="solexa">Solexa (PHRED 0 = &quot;@&quot;)</option>
+    </param>
+    <param name="keep_sscs" type="boolean" truevalue="true" falsevalue="" label="Output single-strand consensus sequences as well" />
+  </inputs>
+  <outputs>
+    <data name="dcs1" format="fasta" label="$tool.name on $on_string (mate 1)"/>
+    <data name="dcs2" format="fasta" label="$tool.name on $on_string (mate 2)"/>
+    <data name="sscs1" format="fasta" label="$tool.name on $on_string (SSCS mate 1)">
+      <filter>keep_sscs</filter>
+    </data>
+    <data name="sscs2" format="fasta" label="$tool.name on $on_string (SSCS mate 2)">
+      <filter>keep_sscs</filter>
+    </data>
+  </outputs>
+  <tests>
+    <test>
+      <param name="input" value="families.msa.tsv"/>
+      <output name="dcs1" file="families.cons_1.fa"/>
+      <output name="dcs2" file="families.cons_2.fa"/>
+    </test>
+  </tests>
+  <help>
+
+**What it does**
+
+This is for processing duplex sequencing data. It creates single-strand and duplex consensus reads from aligned read families.
+
+-----
+
+**Input**
+
+This expects the output format of the "Align families" tool.
+
+-----
+
+**Output**
+
+This will output final, duplex consensus reads in two FASTA files (first and second reads in the pairs). Optionally, you can save the single-strand reads too, in a separate FASTA file.
+
+    </help>
+</tool>