Mercurial > repos > richard-burhans > segalign
changeset 17:52a13f003b2d draft
planemo upload for repository https://github.com/richard-burhans/galaxytools/tree/main/tools/segalign commit b5d08d89e75e2dbfa03f74e45aa6baa465582eee
author | richard-burhans |
---|---|
date | Tue, 30 Jul 2024 17:25:08 +0000 (4 months ago) |
parents | 4966e095e3c3 |
children | a5769923297c |
files | run_segalign_diagonal_partition |
diffstat | 1 files changed, 0 insertions(+), 107 deletions(-) [+] |
line wrap: on
line diff
--- a/run_segalign_diagonal_partition Tue Jul 30 17:19:56 2024 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,107 +0,0 @@ -#!/usr/bin/env bash - -set -o errexit -set -o nounset -set -o pipefail - -## -## parse arguments -## - -SEGALIGN_ARGS=() - -MAX_SEGMENT_SIZE="" -OUTPUT_FILENAME="" -LASTZ_COMMAND_FILE="lastz_commands.txt" -HELP=0 - -while [[ $# -gt 0 ]]; do - case $1 in - --help) - HELP=1 - shift - ;; - --max_segments) - MAX_SEGMENT_SIZE="$2" - shift 2 - ;; - --output) - OUTPUT_FILENAME=$(readlink -f "$2") - shift 2 - ;; - --tool_directory) - TOOL_DIRECTORY="$2" - shift 2 - ;; - *) - SEGALIGN_ARGS+=("$1") - shift - esac -done - -set -- "${SEGALIGN_ARGS[@]}" - -## -## check arguments -## - -if [[ $# == 0 || "$HELP" == "1" ]]; then - segalign --help - exit 0 -fi - -if [[ $# -lt 2 ]]; then - echo "run_segalign_diagonal_partition: missing target and query sequences" 1>&2 - exit 1 -fi - -ref_path=$(readlink -f "$1") -test -e "$ref_path" || { - echo "run_segalign_diagonal_partition: target file \"$ref_path\" does not exist" 1>&2 - exit 1 -} -query_path=$(readlink -f "$2") -test -e "$query_path" || { - echo "run_segalign_diagonal_partition: query file \"$query_path\" does not exist" 1>&2 - exit 1 -} -shift 2 - -DATA_FOLDER="data" -mkdir -p "$DATA_FOLDER" || { - echo "run_segalign_diagonal_partition: cannont create data directory \"$DATA_FOLDER\"" 1>&2 - exit 1 -} - -cd $DATA_FOLDER/.. -echo "" -echo "Converting fasta files to 2bit format" 1>&2 - -## -## convert target and query to 2bit -## -faToTwoBit <(gzip -cdfq "$ref_path") "$DATA_FOLDER/ref.2bit" || { - echo "run_segalign_diagonal_partition: cannot convert \"$ref_path\" to 2bit" 1>&2 - exit 1 -} -faToTwoBit <(gzip -cdfq "$query_path") "$DATA_FOLDER/query.2bit" || { - echo "run_segalign_diagonal_partition: cannont convert \"$ref_path\" to 2bit" 1>&2 - exit 1 -} - - - -time { - while IFS= read -r line; do - python "$TOOL_DIRECTORY/diagonal_partition.py" $MAX_SEGMENT_SIZE $line >> $LASTZ_COMMAND_FILE || { - echo "run_segalign_diagonal_partition: failed: diagonal_partition.py $MAX_SEGMENT_SIZE $line" 1>&2 - exit 1 - } - # segalign sort writes out the partitioned segment files to the working - # directory and prints the modified lastz commands. - done < <(stdbuf -oL segalign $ref_path $query_path "${DATA_FOLDER}/" "$@" ) # segalign begins running in this line, -} 1>&2 # and every newline written to stdout, get assigned to $line which - # gets sent to diagonal_partition for diagonal partitioning - -exit 0 -