Skip to content

Commit

Permalink
Warn if taxonomy is set to "auto" when protein fasta is provided
Browse files Browse the repository at this point in the history
- Currently, Bacteria/Archaea can only be predicted based on the
nucleotide genome sequences, since it is based on 16S rRNA gene
similarity.

- While a predictor from protein fasta will be implemented in the
future, for now it will print a warning and assumed a bacterial genome
input.
  • Loading branch information
Waschina committed Sep 29, 2022
1 parent 780e528 commit 706c49b
Showing 1 changed file with 7 additions and 1 deletion.
8 changes: 7 additions & 1 deletion src/gapseq_find.sh
Original file line number Diff line number Diff line change
Expand Up @@ -227,7 +227,7 @@ if [[ "$user_temp" = true ]]; then
tmpdir=$(mktemp -d $user_temp_folder/"$tmp_fasta"_XXXXXX)
else
tmpdir=$(mktemp -d)
#trap 'rm -rf "$tmpdir"' EXIT
trap 'rm -rf "$tmpdir"' EXIT
fi
echo $tmpdir
cd $tmpdir
Expand Down Expand Up @@ -302,6 +302,12 @@ case $pathways in
;;
esac

# Taxonomy prediction does currently work only with nuceotide fasta input
if [ $input_mode == "prot" ] && [ $taxonomy == "auto" ]; then
echo "WARNING: Automated taxonomy prediction is not yet implemented for protein fasta input. Assuming Bacteria..."
taxonomy="Bacteria"
fi

# determine taxonomy
if [ "$taxonomy" == "auto" ]; then
pred_biom=$($dir/predict_biomass_from16S.sh "$fasta")
Expand Down

0 comments on commit 706c49b

Please sign in to comment.