Repository revision
36:fc35ffc8c548

Repository 'ncbi_blast_plus'
hg clone https://toolshed.g2.bx.psu.edu/repos/devteam/ncbi_blast_plus

NCBI BLAST+ makeblastdb tool metadata
Miscellaneous
Make BLAST database
ncbi_makeblastdb
toolshed.g2.bx.psu.edu/repos/devteam/ncbi_blast_plus/ncbi_makeblastdb/2.16.0+galaxy0
2.16.0+galaxy0
makeblastdb -version
True
Version lineage of this tool (guids ordered most recent to oldest)
toolshed.g2.bx.psu.edu/repos/devteam/ncbi_blast_plus/ncbi_makeblastdb/2.16.0+galaxy0 (this tool)
toolshed.g2.bx.psu.edu/repos/devteam/ncbi_blast_plus/ncbi_makeblastdb/0.0.1
ncbi_makeblastdb
Requirements (dependencies defined in the <requirements> tag set)
name version type
blast 2.16.0 package
python 3.9 package
Additional information about this tool
#set $inputs = []
#set $input_compression = []
#for r in $input.selection:
    #if $input.type == "protein":
        #silent $inputs.append($r.input_file)
        #silent $input_compression.append($r.input_file.is_of_type('fasta.gz'))
    #elif $r.nuc_choice.source == "history":
        #silent $inputs.append($r.nuc_choice.input_file)
        #silent $input_compression.append($r.nuc_choice.input_file.is_of_type('fasta.gz'))
    #else:
        #silent $inputs.append($r.nuc_choice.input_file.fields.path)
        #silent $input_compression.append(False)
    #end if
#end for

python $__tool_directory__/check_no_duplicates.py
##First check for duplicates (since BLAST+ 2.2.28 fails to do so)
##and abort (via the ampersand ampersand trick) if any are found.
#for i in $inputs#'$i' #end for#
&&
##makeblastdb does not like input redirects of the sort
##makeblastdb -in <(gunzip -c gzipped_fasta_file)
##therefore we're cramming everything
##into a single cat command below
cat
#for i, is_gzipped in zip($inputs, $input_compression):
    #if $is_gzipped:
        <(gunzip -c '$i')
    #else:
        '$i'
    #end if
#end for
| makeblastdb -out '${os.path.join($outfile.files_path, "blastdb")}'
-blastdb_version 4
$parse_seqids
$hash_index
-in -
#if $title:
-title '${title}'
#else:
##Would default to being based on the cryptic Galaxy filenames, which is unhelpful
-title 'BLAST Database'
#end if
-dbtype
#if $input.type == "protein":
    prot
#else:
    nucl
#end if
## --------------------------------------------------------------------
## Masking
## --------------------------------------------------------------------
## HACK: If no mask files, evaluates as a list with just None in it:
## See Trello issue https://trello.com/c/lp5YmA1O
#if ' '.join( map(str, $mask_data_file) ) != 'None':
#for i in $mask_data_file:
-mask_data '${i}'
#end for
#end if
## --------------------------------------------------------------------
## Taxonomy
## --------------------------------------------------------------------
#if $tax.taxselect == 'id':
-taxid $tax.taxid
## #else if $tax.taxselect == 'map':
## -taxid_map $tax.taxmap
#end if
## --------------------------------------------------------------------
## Capture the stdout log information to the primary file (plain text):
> '$outfile'
    
None
False
Functional tests
name inputs outputs required files
Test-1 input|selection_0|input_file: four_human_proteins.fasta
input|type: protein
title: Just 4 human proteins
parse_seqids: False
hash_index: True
name: value
four_human_proteins.fasta
value
Test-2 input|selection_0|input_file: four_human_proteins.fasta
input|type: protein
title: Just 4 human proteins
parse_seqids: False
hash_index: True
tax|taxid: 9606
tax|taxselect: id
name: value
four_human_proteins.fasta
value
Test-3 input|selection_0|input_file: four_human_proteins.fasta
input|type: protein
title: Just 4 human proteins
parse_seqids: False
hash_index: True
mask_data_file: segmasker_four_human.maskinfo-asn1
name: value
four_human_proteins.fasta
segmasker_four_human.maskinfo-asn1
value
Test-4 input|selection_0|nuc_choice|input_file: three_human_mRNA.fasta.gz
input|selection_0|nuc_choice|source: history
input|type: nucleotide
title: Just 3 human mRNA sequences
parse_seqids: False
hash_index: True
tax|taxid: 9606
tax|taxselect: id
name: value
three_human_mRNA.fasta.gz
value
Test-5 input|selection_0|nuc_choice|input_file: three_human_mRNA
input|selection_0|nuc_choice|source: cached
input|type: nucleotide
title: Just 3 human mRNA sequences
parse_seqids: False
hash_index: True
tax|taxid: 9606
tax|taxselect: id
name: value
value