Workflow: retrieve sequence and perform pairwise alignment (sub-workflow process)

Fetched 2026-10-11 08:18:46 GMT

\"Perform pairwise alignment of protein sequences for pairs identified by structural similarity search. Step 1: blastdbcmd: ../Tools/14_blastdbcmd.cwl Step 2: seqretsplit: ../Tools/15_seqretsplit.cwl Step 3: needle (Global alignment): ../Tools/16_needle.cwl Step 4: water (Local alignment): ../Tools/16_water.cwl \"

children parents
Workflow as SVG
  • Selected
  • Default Values
  • Nested Workflows
  • Tools
  • Inputs/Outputs

Inputs

ID Type Title Doc
BLAST_INDEX_FILES File blast index files (blastdbcmd)

blast index files for blastdbcmd

FOLDSEEK_EXTRACT_TSV File [TSV] foldseek extract tsv (foldseek easy-search)

foldseek extract tsv

WATER_RESULT_DIR_NAME String [Directory name] water result directory name (water)

water result directory name

NEEDLE_RESULT_DIR_NAME String [Directory name] needle result directory name (needle)

needle result directory name

ENTRY_BATCH_HIT_SPECIES File entry batch file (blastdbcmd)

entry batch file for blastdbcmd

ENTRY_BATCH_QUERY_SPECIES File entry batch file (blastdbcmd)

entry batch file for blastdbcmd

ALIGNMENT_QUERY_COLUMN_NUMBER Integer alignment query column number (needle and water)

alignment column number (query species) for needle and water. Extract columns describing UniProt ID pairs (query IDs) from the TSV file read with the FOLDSEEK_EXTRACT_TSV parameter

ALIGNMENT_TARGET_COLUMN_NUMBER Integer alignment target column number (needle and water)

alignment column number (target species) for needle and water. Extract columns describing UniProt ID pairs (hit IDs) from the TSV file read with the FOLDSEEK_EXTRACT_TSV parameter

BLASTDBCMD_LOGFILE_NAME_HIT_SPECIES String [File name] logfile name (blastdbcmd)

logfile name.

BLASTDBCMD_LOGFILE_NAME_QUERY_SPECIES String [File name] logfile name (blastdbcmd)

logfile name.

BLASTDBCMD_RESULT_FILE_NAME_HIT_SPECIES String [File name] blastdbcmd result file name (blastdbcmd)

blastdbcmd result file name.

SEQRETSPLIT_OUTPUT_DIR_NAME_HIT_SPECIES String [Directory name] output directory name (seqretsplit)

output directory name for seqretsplit

BLASTDBCMD_RESULT_FILE_NAME_QUERY_SPECIES String [File name] blastdbcmd result file name (blastdbcmd)

blastdbcmd result file name.

SEQRETSPLIT_OUTPUT_DIR_NAME_QUERY_SPECIES String [Directory name] output directory name (seqretsplit)

output directory name for seqretsplit

Steps

ID Runs Label Doc
blastdbcmd_hit_species
../Tools/14_blastdbcmd.cwl (CommandLineTool)
seqretsplit_hit_species
../Tools/15_seqretsplit.cwl (CommandLineTool)
seqretsplit command for split fasta file

seqretsplit command for split fasta file which is created by /Tools/14_blastdbcmd.cwl. Before executing, make sure the blastdbcmd result file is already created by /Tools/14_blastdbcmd.cwl

blastdbcmd_query_species
../Tools/14_blastdbcmd.cwl (CommandLineTool)
seqretsplit_query_species
../Tools/15_seqretsplit.cwl (CommandLineTool)
seqretsplit command for split fasta file

seqretsplit command for split fasta file which is created by /Tools/14_blastdbcmd.cwl. Before executing, make sure the blastdbcmd result file is already created by /Tools/14_blastdbcmd.cwl

local_alignment_using_water
../Tools/16_water.cwl (CommandLineTool)
global_alignment_using_needle
../Tools/16_needle.cwl (CommandLineTool)
needle command for global alignment

needle command for global alignment using the foldseek easy-search result TSV file and the split FASTA files.

Outputs

ID Type Label Doc
output_dir_hit_species Directory output directory (hit species)
output_water_result_dir Directory water result directory
output_dir_query_species Directory output directory (query species)
output_needle_result_dir Directory needle result directory
output_water_result_file File[] water result file
output_needle_result_file File[] needle result file
output_logfile_hit_species File logfile (hit species)
output_logfile_query_species File logfile (query species)
output_blastdbcmd_result_hit_species File blastdbcmd result (hit species)
output_split_fasta_files_hit_species File[] split fasta files (hit species)
output_blastdbcmd_result_query_species File blastdbcmd result (query species)
output_split_fasta_files_query_species File[] split fasta files (query species)
Permalink: https://w3id.org/cwl/view/git/62a2b6767507ea65ddd36b853b539060840bde48/Workflow/11_retrieve_sequence_swf.cwl