Repository revision
2:1ec42f033bb4

Repository 'refseq_masher'
hg clone https://toolshed.g2.bx.psu.edu/repos/nml/refseq_masher

RefSeq Masher Matches tool metadata
Miscellaneous
Find closest matching NCBI RefSeq Genomes to your sequences
refseq_masher_matches
toolshed.g2.bx.psu.edu/repos/nml/refseq_masher/refseq_masher_matches/0.1.2
0.1.2
None
True
Version lineage of this tool (guids ordered most recent to oldest)
toolshed.g2.bx.psu.edu/repos/nml/refseq_masher/refseq_masher_matches/0.1.2 (this tool)
toolshed.g2.bx.psu.edu/repos/nml/refseq_masher/refseq_masher_matches/0.1.1
refseq_masher_matches
Requirements (dependencies defined in the <requirements> tag set)
name version type
refseq_masher 0.1.2 package
Additional information about this tool
#import re
#import os

#if $input.type == 'fasta'
#set $input_files = '"{}.fasta"'.format(os.path.splitext($input.fasta.name)[0])
  ln -s "$input.fasta" $input_files &&
#elif $input.type == 'paired'
#set $_forward_ext = '.fastq.gz' if $re.match(r'.*\.gz$', $input.forward.name) else '.fastq'
#set $_forward = '"{}_1{}"'.format($re.sub(r'_[12]\..+$', '', $input.forward.name), $_forward_ext)
#set $_reverse_ext = '.fastq.gz' if $re.match(r'.*\.gz$', $input.reverse.name) else '.fastq'
#set $_reverse = '"{}_2{}"'.format($re.sub(r'_[12]\..+$', '', $input.reverse.name), $_reverse_ext)
#set $input_files = '{} {}'.format($_forward, $_reverse)
  ln -s "$input.forward" $_forward &&
  ln -s "$input.reverse" $_reverse &&
#elif $input.type == 'single'
#set $input_files = '"{}"'.format($input.single.name)
  ln -s "$input.single" $input_files &&
#elif $input.type == 'paired_collection'
#set $_forward_ext = '.fastq.gz' if $re.match(r'.*\.gz$', str($input.paired_collection.forward)) else '.fastq'
#set $_forward = '"{}_1{}"'.format($input.paired_collection.name, $_forward_ext)
#set $_reverse_ext = '.fastq.gz' if $re.match(r'.*\.gz$', str($input.paired_collection.reverse)) else '.fastq'
#set $_reverse = '"{}_2{}"'.format($input.paired_collection.name, $_reverse_ext)
#set $input_files = '{} {}'.format($_forward, $_reverse)
  ln -s "$input.paired_collection.forward" $_forward &&
  ln -s "$input.paired_collection.reverse" $_reverse &&
#end if

refseq_masher 
  $adv.verbosity 
  matches 
  --output refseq_masher-matches.${adv.output_type}
  --output-type $adv.output_type
  --top-n-results $top_n_results
#if $adv.min_kmer_threshold
  --min-kmer-threshold $adv.min_kmer_threshold
#end if
  -T "\${TMPDIR:-/tmp}"
  $input_files

  
None
False
Functional tests
name inputs outputs required files
Test-1 input|fasta: Se-Enteritidis.fasta
input|type: fasta
top_n_results: 1
adv|output_type: tab
name: value
Se-Enteritidis.fasta
value
Test-2 input|single: SRR1203042_1-head4000.fastq
input|type: single
top_n_results: 1
adv|min_kmer_threshold: 2
adv|output_type: tab
name: value
SRR1203042_1-head4000.fastq
value