# HG changeset patch # User iuc # Date 1580179849 18000 # Node ID e4d0231d8595bf17d59a6a1d392984f49888b5e2 # Parent 32f2211eeec3e3ec6c5eb8ede5724f257335125b "planemo upload commit 02413e02c9484cdb604bda2c09931599f9259792" diff -r 32f2211eeec3 -r e4d0231d8595 snippy-core.xml --- a/snippy-core.xml Sun Dec 01 13:06:05 2019 -0500 +++ b/snippy-core.xml Mon Jan 27 21:50:49 2020 -0500 @@ -1,5 +1,5 @@ - + Combine multiple Snippy outputs into a core SNP alignment @@ -9,14 +9,14 @@ diff -r 32f2211eeec3 -r e4d0231d8595 snippy.xml --- a/snippy.xml Sun Dec 01 13:06:05 2019 -0500 +++ b/snippy.xml Mon Jan 27 21:50:49 2020 -0500 @@ -1,4 +1,4 @@ - + Snippy finds SNPs between a haploid reference genome and your NGS sequence reads. @@ -12,8 +12,19 @@ @REFERENCE_SOURCE_FILE@ + #import re + #if str( $fastq_input.fastq_input_selector ) == "paired" + #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input1.element_identifier) + #elif str( $fastq_input.fastq_input_selector ) == "paired_collection" + #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input.name) + #elif str( $fastq_input.fastq_input_selector ) == "single" + #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input_single.element_identifier) + #elif str( $fastq_input.fastq_input_selector ) == "paired_iv" + #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input_interleaved.element_identifier) + #end if + snippy - --outdir 'out' + --outdir '${dir_name}' --cpus \${GALAXY_SLOTS:-1} --ram \$((\${GALAXY_MEMORY_MB:-4096}/1024)) @REFERENCE_COMMAND@ @@ -40,26 +51,12 @@ --peil '$fastq_input.fastq_input_interleaved' #end if - && - - #import re - #if str( $fastq_input.fastq_input_selector ) == "paired" - #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input1.element_identifier) - #elif str( $fastq_input.fastq_input_selector ) == "paired_collection" - #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input.name) - #elif str( $fastq_input.fastq_input_selector ) == "single" - #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input_single.element_identifier) - #elif str( $fastq_input.fastq_input_selector ) == "paired_iv" - #set $dir_name = re.sub('[^\w_]', '_', $fastq_input.fastq_input_interleaved.element_identifier) + #if "outcon" in str($outputs) and $adv.rename_cons + && sed -i 's/>.*/>${dir_name}/' '${dir_name}/snps.consensus.fa' #end if - mkdir -p ${dir_name} && cp -r out/reference out/snps.tab out/snps.aligned.fa out/snps.vcf ${dir_name}/ && - tar -czf out.tgz ${dir_name} - #if "outcon" in str($outputs) and $adv.rename_cons - && sed -i 's/>.*/>${dir_name}/' out/snps.consensus.fa - #end if - - + && cp -r '${dir_name}' out + && tar -czf out.tgz '${dir_name}' ]]>