diff seqtk_seq.xml @ 0:878ec9a5f8c0 draft default tip

planemo upload
author matthias
date Tue, 08 Jan 2019 09:43:35 -0500
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/seqtk_seq.xml	Tue Jan 08 09:43:35 2019 -0500
@@ -0,0 +1,95 @@
+<tool id="seqtk_seq" name="Convert to FASTA (seqtk)" version="0.1.0">
+    <requirements>
+        <requirement type="package" version="1.2">seqtk</requirement>
+    </requirements>
+    <command detect_errors="exit_code"><![CDATA[
+seqtk seq 
+-a
+$shift_quality 
+-q $quality_min
+-X $quality_max
+#if $mask_regions
+    -M '$maskregions'
+#end if
+#if $sample.sample_selector
+    -f $sample.fraction
+    #if str($sample.seed)!=''
+        -s $sample.seed
+    #end if
+#end if
+
+'$input1' 
+> 
+'$output1'
+    ]]></command>
+    <inputs>
+        <param type="data" name="input1" format="fastq" />
+        <param name="shift_quality" type="boolean" label="Shift quality"
+               truevalue="-V" falsevalue=""
+               help="shift quality by '(-Q) - 33' (-V)" />
+        <param name="quality_min" type="integer" label="Mask bases with quality lower than"
+               value="0" min="0" max="255" help="(-q)" />
+        <param name="quality_max" type="integer" label="Mask bases with quality higher than"
+               value="255" min="0" max="255" help="(-X)" />
+        <param name="mask_regions" type="data" label="Mask regions in BED"
+		format="bed" help="(-M)" optional="true" />
+        <conditional name="sample">
+            <param name="sample_selector" type="boolean" label="Sample fraction of sequences" />
+            <when value="true">
+                <param name="fraction" label="Fraction" type="float" value="1.0"
+                       help="(-f)" />
+                <param name="seed" label="Random seed" type="integer" value=""
+                       help="(-s)" optional="true" />
+            </when>
+            <when value="false">
+            </when>
+        </conditional>
+    </inputs>
+
+    <outputs>
+        <data name="output1" format="fasta" />
+    </outputs>
+    <tests>
+        <test>
+            <param name="input1" value="2.fastq"/>
+            <output name="output1" file="2.fasta"/>
+        </test>
+        <test>
+            <param name="input1" value="2.fastq"/>
+            <param name="shift_quality" value="-V"/>
+            <param name="quality_min" value="30"/>
+            <param name="quality_max" value="31"/>
+            <output name="output1" file="2.fasta"/>
+        </test>
+    </tests>
+    <help><![CDATA[
+        
+Usage:   seqtk seq [options] <in.fq>|<in.fa>
+
+Options: -q INT    mask bases with quality lower than INT [0]
+         -n CHAR   masked bases converted to CHAR; 0 for lowercase [0]
+         -l INT    number of residues per line; 0 for 2^32-1 [0]
+         -Q INT    quality shift: ASCII-INT gives base quality [33]
+         -s INT    random seed (effective with -f) [11]
+         -f FLOAT  sample FLOAT fraction of sequences [1]
+         -M FILE   mask regions in BED or name list FILE [null]
+         -L INT    drop sequences with length shorter than INT [0]
+         -c        mask complement region (effective with -M)
+         -r        reverse complement
+         -A        force FASTA output (discard quality)
+         -C        drop comments at the header lines
+
+
+    ]]></help>
+    <citations>
+        <citation type="bibtex">
+@misc{githubseqtk,
+  author = {LastTODO, FirstTODO},
+  year = {TODO},
+  title = {seqtk},
+  publisher = {GitHub},
+  journal = {GitHub repository},
+  url = {https://github.com/lh3/seqtk},
+}</citation>
+    </citations>
+</tool>