Skip to content

Commit

Permalink
GTDB ar122_taxonomy does not exist anymore, replace with different file
Browse files Browse the repository at this point in the history
  • Loading branch information
milot-mirdita committed Jul 4, 2022
1 parent 7be78c8 commit 678c82a
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions data/workflow/databases.sh
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,7 @@ case "${SELECTION}" in
downloadFile "https://data.ace.uq.edu.au/public/gtdb/data/releases/latest/VERSION" "${TMP_PATH}/version"
downloadFile "https://data.ace.uq.edu.au/public/gtdb/data/releases/latest/genomic_files_reps/gtdb_proteins_aa_reps.tar.gz" "${TMP_PATH}/gtdb.tar.gz"
downloadFile "https://data.ace.uq.edu.au/public/gtdb/data/releases/latest/bac120_taxonomy.tsv" "${TMP_PATH}/bac120_taxonomy.tsv"
downloadFile "https://data.ace.uq.edu.au/public/gtdb/data/releases/latest/ar122_taxonomy.tsv" "${TMP_PATH}/ar122_taxonomy.tsv"
downloadFile "https://data.ace.uq.edu.au/public/gtdb/data/releases/latest/ar53_taxonomy.tsv" "${TMP_PATH}/ar53_taxonomy.tsv"
touch "${TMP_PATH}/download.done"
fi
INPUT_TYPE="GTDB"
Expand Down Expand Up @@ -460,7 +460,7 @@ if [ -n "${TAXONOMY}" ] && notExists "${OUTDB}_mapping"; then
printf("%s\t%s\n", $1, ids[$NF]) > taxdir"/mapping_genomes";
}'
mkdir -p "${TMP_PATH}/taxonomy"
awk -v taxdir="${TMP_PATH}/taxonomy" "$CMD" "${TMP_PATH}/bac120_taxonomy.tsv" "${TMP_PATH}/ar122_taxonomy.tsv"
awk -v taxdir="${TMP_PATH}/taxonomy" "$CMD" "${TMP_PATH}/bac120_taxonomy.tsv" "${TMP_PATH}/ar53_taxonomy.tsv"
touch "${TMP_PATH}/taxonomy/merged.dmp"
touch "${TMP_PATH}/taxonomy/delnodes.dmp"
# shellcheck disable=SC2086
Expand Down

0 comments on commit 678c82a

Please sign in to comment.