Repository revision
11:8f9f06995f98

Repository 'unicycler'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/unicycler

Create assemblies with Unicycler tool metadata
Miscellaneous
pipeline for bacterial genomes
unicycler
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.5.0+galaxy1
0.5.0+galaxy1
None
True
Version lineage of this tool (guids ordered most recent to oldest)
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.5.0+galaxy1 (this tool)
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.5.0+galaxy0
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.4.8.0
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.4.7.0
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.4.6.0
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.4.4.0
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.4.1.1
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.4.1
toolshed.g2.bx.psu.edu/repos/iuc/unicycler/unicycler/0.2.0
unicycler
Requirements (dependencies defined in the <requirements> tag set)
name version type
unicycler 0.5.0 package
samtools 1.15.1 package
Additional information about this tool
#for r in $reuse
    ln -s $r.reuse_file ${r.reuse_step}.gfa &&
#end for

## Preparing files
#set $uncompressed = ('fastqsanger','fastq')
#set $compressed = ('fastqsanger.gz','fastq.gz')
#if str( $paired_unpaired.fastq_input_selector ) == "paired"
    #if $paired_unpaired.fastq_input1.file_ext in $uncompressed
        #set fq1 = "fq1.fastq"
    #elif $paired_unpaired.fastq_input1.file_ext in $compressed
        #set fq1 = "fq1.fastq.gz"
    #end if
    #if $paired_unpaired.fastq_input2.file_ext in $uncompressed
        #set fq2 = "fq2.fastq"
    #elif $paired_unpaired.fastq_input2.file_ext in $compressed
        #set fq2 = "fq2.fastq.gz"
    #end if
    ln -s '${paired_unpaired.fastq_input1}' $fq1 &&
    ln -s '${paired_unpaired.fastq_input2}' $fq2 &&
#elif str( $paired_unpaired.fastq_input_selector ) == "paired_collection"
    #if $paired_unpaired.fastq_input1.forward.file_ext in $uncompressed
        #set fq1 = "fq1.fastq"
    #elif $paired_unpaired.fastq_input1.forward.file_ext in $compressed
        #set fq1 = "fq1.fastq.gz"
    #end if
    #if $paired_unpaired.fastq_input1.reverse.file_ext in $uncompressed
        #set fq2 = "fq2.fastq"
    #elif $paired_unpaired.fastq_input1.reverse.file_ext in $compressed
        #set fq2 = "fq2.fastq.gz"
    #end if
    ln -s '${paired_unpaired.fastq_input1.forward}' $fq1 &&
    ln -s '${paired_unpaired.fastq_input1.reverse}' $fq2 &&
#elif str( $paired_unpaired.fastq_input_selector ) == "single"
    #if $paired_unpaired.fastq_input1.file_ext in $uncompressed
        #set fq = "fq.fastq"
    #elif $paired_unpaired.fastq_input1.file_ext in $compressed
        #set fq = "fq.fastq.gz"
    #end if
    ln -s '${paired_unpaired.fastq_input1}' '$fq' &&
#end if
#if $long
    #if $long.file_ext in $uncompressed
        #set lr = "lr.fastq"
    #elif $long.file_ext in $compressed
        #set lr = "lr.fastq.gz"
    #elif $long.is_of_type('fasta')
        #set lr = "lr.fasta"
    #end if
    ln -s '${long}' '$lr' &&
#end if
## Running Unicycler
unicycler -t "\${GALAXY_SLOTS:-4}"
-o ./
--verbosity 3
#if str( $paired_unpaired.fastq_input_selector ) == "paired"
    -1 '$fq1'
    -2 '$fq2'
#elif str( $paired_unpaired.fastq_input_selector ) == "paired_collection"
    -1 '$fq1'
    -2 '$fq2'
#elif str( $paired_unpaired.fastq_input_selector ) == "single"
    -s '$fq'
#end if
#if $long
    -l $lr
#end if
## General Unicycler Options section
## ----------------------------------------------------------
--mode '$mode'
--min_fasta_length '$min_fasta_length'
--linear_seqs '$linear_seqs'
#if str($min_anchor_seg_len) != ''
--min_anchor_seg_len '$min_anchor_seg_len'
#end if
## Spades Options section
## ----------------------------------------------------------
--min_kmer_frac '$spades.min_kmer_frac'
--max_kmer_frac '$spades.max_kmer_frac'
#if str($spades.kmers) != ''
--kmers '$spades.kmers'
#end if
--kmer_count '$spades.kmer_count'
--depth_filter '$spades.depth_filter'
#if $spades.largest_component
    --largest_component
#end if
## Rotation Options section
## ----------------------------------------------------------
$rotation.no_rotate
#if $rotation.start_genes
    --start_genes '$rotation.start_genes'
#end if
--start_gene_id '$rotation.start_gene_id'
--start_gene_cov '$rotation.start_gene_cov'
## Graph cleaning Options sdection
## ----------------------------------------------------------
--min_component_size '$graph_clean.min_component_size'
--min_dead_end_size '$graph_clean.min_dead_end_size'
## Long Read Alignment Options
## ----------------------------------------------------------
#if $lr_align.contamination
    --contamination '$lr_align.contamination'
#end if
--scores '${lr_align.scores}'
#if str($lr_align.low_score) != ''
    --low_score '$lr_align.low_score'
#end if
$lr_align.no_simple_bridges
--keep $keep
#if $keep != '0'
    && mkdir 'spades_graphs'
    && mv 00*gfa './spades_graphs/'
#end if
#if $keep == '2' and $long
        && samtools view -@ "\${GALAXY_SLOTS:-4}" -u 'read_alignment/long_read_alignments.sam' | samtools sort -@ "\${GALAXY_SLOTS:-4}" -o 'read_alignment/long_read_alignments.bam'
#end if
    
None
False
Functional tests
name inputs outputs required files
Test-1 paired_unpaired|fastq_input1: phix_f.fq.gz
paired_unpaired|fastq_input2: phix_r.fq.gz
paired_unpaired|fastq_input_selector: paired
mode: normal
min_fasta_length: 100
linear_seqs: 0
spades|min_kmer_frac: 0.2
spades|max_kmer_frac: 0.95
spades|kmer_count: 10
spades|depth_filter: 0.25
rotation|no_rotate: False
rotation|start_gene_id: 90
rotation|start_gene_cov: 95
graph_clean|min_component_size: 1000
graph_clean|min_dead_end_size: 1000
lr_align|scores: 3,-6,-5,-2
keep: 0
name: value
name: value
phix_f.fq.gz
phix_r.fq.gz
value
Test-2 paired_unpaired|fastq_input1: phix_f.fq.gz
paired_unpaired|fastq_input2: phix_r.fq.gz
paired_unpaired|fastq_input_selector: paired
long: onp.fa
mode: normal
min_fasta_length: 100
linear_seqs: 0
spades|min_kmer_frac: 0.2
spades|max_kmer_frac: 0.95
spades|kmer_count: 10
spades|depth_filter: 0.25
rotation|no_rotate: False
rotation|start_gene_id: 90
rotation|start_gene_cov: 95
graph_clean|min_component_size: 1000
graph_clean|min_dead_end_size: 1000
lr_align|scores: 3,-6,-5,-2
keep: 0
name: value
name: value
phix_f.fq.gz
phix_r.fq.gz
onp.fa
value
Test-3 paired_unpaired|fastq_input1: paired collection
paired_unpaired|fastq_input_selector: paired_collection
mode: normal
min_fasta_length: 100
linear_seqs: 0
spades|min_kmer_frac: 0.2
spades|max_kmer_frac: 0.95
spades|kmer_count: 10
spades|depth_filter: 0.25
rotation|no_rotate: False
rotation|start_gene_id: 90
rotation|start_gene_cov: 95
graph_clean|min_component_size: 1000
graph_clean|min_dead_end_size: 1000
lr_align|scores: 3,-6,-5,-2
keep: 0
name: value
name: value
phix_f.fq.gz
phix_r.fq.gz
value
Test-4 paired_unpaired|fastq_input_selector: none
long: only_long.fasta
min_anchor_seg_len: 10
spades|kmers: 21,23
keep: 0
name: value
name: value
only_long.fasta
value
Test-5 paired_unpaired|fastq_input1: paired collection
paired_unpaired|fastq_input_selector: paired_collection
long: only_long.fasta
keep: 0
reuse_0|reuse_file: phix__spades_graph.gfa1
reuse_0|reuse_step: 002_depth_filter
name: value
name: value
phix_f.fq.gz
phix_r.fq.gz
only_long.fasta
phix__spades_graph.gfa1
value
Test-6 paired_unpaired|fastq_input1: phix_f.fq.gz
paired_unpaired|fastq_input2: phix_r.fq.gz
paired_unpaired|fastq_input_selector: paired
mode: normal
keep: 1
name: value
name: value
phix_f.fq.gz
phix_r.fq.gz
value
Test-7 paired_unpaired|fastq_input1: phix_f.fq.gz
paired_unpaired|fastq_input2: phix_r.fq.gz
paired_unpaired|fastq_input_selector: paired
long: onp.fa
mode: normal
keep: 2
name: value
name: value
name: value
phix_f.fq.gz
phix_r.fq.gz
onp.fa
value
Test-8 paired_unpaired|fastq_input1: phix_f.fq.gz
paired_unpaired|fastq_input2: phix_r.fq.gz
paired_unpaired|fastq_input_selector: paired
long: onp.fa
mode: normal
lr_align|no_simple_bridges: True
keep: 0
name: value
name: value
phix_f.fq.gz
phix_r.fq.gz
onp.fa
value