Mercurial > repos > cstrittmatter > skesa
changeset 20:a368d48a4562 draft
planemo upload commit cd1454c40a43ad9da3d59e6ba8359318fc772c43-dirty
author | cstrittmatter |
---|---|
date | Thu, 13 Sep 2018 09:34:22 -0400 |
parents | abb622c228c4 |
children | 8bafd3d18864 |
files | run_skesa_collection.py skesa.xml |
diffstat | 2 files changed, 11 insertions(+), 83 deletions(-) [+] |
line wrap: on
line diff
--- a/run_skesa_collection.py Mon Aug 27 08:59:02 2018 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,43 +0,0 @@ -import sys - -from pathlib import Path - -path1 = sys.argv[1] - -cores = sys.argv[2] - -memory = sys.argv[2] - -print(path1) - -#print(path2) - -pathlist = Path(path1).glob('*fastq*') - -i = 0 -skcmd = [] - -skcmd.append("skesa ") - -for path in pathlist: - # because path is object not string-wq - i += 1 - path_in_str = str(path) - if i == 2: - skcmd.append("," + path_in_str) - #print("," + path_in_str) - i = 0 - if i == 1: - skcmd.append(" --fastq " + path_in_str) - #print(" --fastq " + path_in_str) - - - -if cores != 0: - skcmd.append(" --cores " + cores) - - -skcmd.append(" --memory " + memory) - -skcmds = ''.join(skcmd) -print(skcmds)
--- a/skesa.xml Mon Aug 27 08:59:02 2018 -0400 +++ b/skesa.xml Thu Sep 13 09:34:22 2018 -0400 @@ -1,40 +1,22 @@ <tool id="skesa" name="skesa" version="0.1"> <requirements> - <requirement type="package" version="2.2">skesa</requirement> - <requirement type="package" >python</requirement> + <requirement type="package" version="2.2">skesa</requirement> </requirements> <command detect_errors="exit_code"><![CDATA[ - - - #if $jobtype.select != "cl" + #if $jobtype.select != "cl" skesa - #if $jobtype.select == "asm" + #if $jobtype.select == "asm" --fasta $draft - #else if $jobtype.select == "se" + #else if $jobtype.select == "se" --fastq $fastq1 - #else if $jobtype.select == "pe" - --fastq $fastq1,$fastq2 --use_paired_ends - #else if $jobtype.select == "rp" - --fastq $jobtype.pairedf.forward,$jobtype.pairedf.reverse --use_paired_ends - #if $cores != 0 - --cores $cores - #end if - --memory $memory > results.skesa.fasta + #else if $jobtype.select == "pe" + --fastq $fastq1,$fastq2 --use_paired_ends #end if - - - #if $jobtype.select == "cl" - - #set $pathOuput = 0 - mkdir temp_data; - #for $collection_file in $jobtype.collection_files - #silent sys.stderr.write($collection_file+"\n") - cp $collection_file temp_data/.; - - #end for - - + #if $cores != 0 + --cores $cores #end if + --memory $memory > results.skesa.fasta + #end if ]]></command> <inputs> @@ -43,12 +25,7 @@ <option value="asm">Genome Assembly</option> <option value="se">Single-End Reads</option> <option value="pe">Paired-End Reads (Separate Files)</option> - <option value="rp">Paired-End Reads (Paired Data Set)</option> - <option value="cl">Collection of Reads</option> </param> - <when value="srr"> - <param name="srrnum" type="text" label="Sra run number"/> - </when> <when value="asm"> <param name="draft" type="data" format="fasta" label="FASTA" /> </when> @@ -59,16 +36,10 @@ <param name="fastq1" type="data" format="fastq" label="FASTQ" /> <param name="fastq2" type="data" format="fastq" label="FASTQ" /> </when> - <when value="cl"> - <param name="collection_files" type="data_collection" format="fastq,fastqsanger,fastq.gz,fastqsanger.gz,fastq.bz2,fastqsanger.bz2" collection_type="paired" label="FASTQS: Must be a Data Set list built from multiple fastq files" /> - </when> - <when value="rp"> - <param name="pairedf" type="data_collection" collection_type="paired" format="fastq" label="FASTQS: Must be a paired set of forward and reverse fastq files" /> - </when> </conditional> <param name="memory" type="integer" label="Memory available (GB) [integer]" value="16" /> <param name="cores" type="integer" label="Number of cores to use (default all) [integer]" value="0" /> - + </inputs> <outputs> <data format="fasta" label="skesa Results" name="${input.name}.skesa.fasta" from_work_dir="*.fasta"/>