Commit ba55f354 by Eric Coissac

Add logs to print the sequence length and if the sequence is reverse

complemented
parent 213735f5
......@@ -23,7 +23,6 @@ source "${THIS_DIR}/../../../scripts/bash_init.sh"
pushTmpDir ORG.its
loginfo "Normalizing nuclear rDNA cistron..."
RRNADB="${NUCRRNA_DATA_DIR}/plants/nuc_RRNA.hmm"
if [[ ! "$1" =~ ^/ ]]; then
......@@ -32,6 +31,7 @@ pushTmpDir ORG.its
QUERY="$1"
fi
loginfo "Sequence length $(seqlength ${QUERY})"
strand=( $(hmmsearch --max ${RRNADB} ${QUERY} | \
$AwkCmd '/Query: / { \
......@@ -55,7 +55,8 @@ pushTmpDir ORG.its
if [[ "${strand[0]}" == "Forward" ]] ; then
cat ${QUERY}
else
fastarevcomp -f ${QUERY}
loginfo "Revert complement rDNA cluster cistron"
fastarevcomp ${QUERY}
fi
else
logerror "Cannot determine the Cistron orientation"
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment