diff --git a/SimiPlot.sh b/SimiPlot.sh index 2a4de599a67e81ef32c1cccc69022db6ee01a865..d3a36c080948e6651b066f079afe6e4cfffa66c3 100755 --- a/SimiPlot.sh +++ b/SimiPlot.sh @@ -32,7 +32,11 @@ # = VERSIONS = # # ============ # # # - VERSION=1.0.220321ac # + VERSION=1.1.230220ac # +# + updating finalizers for BLAST+ version >= 2.13.0 # +# + fixed bug when the reference contains multiple sequences # +# # +# VERSION=1.0.220321ac # # # ############################################################################################################## @@ -264,7 +268,7 @@ TMP=$(randfile $OUTFILE.tmp); ## defining traps ############################################################################################ finalizer() { - rm -f $DB $DB.ndb $DB.nhr $DB.nin $DB.not $DB.nsq $DB.ntf $DB.nto ; + rm -f $DB $DB.ndb $DB.nhr $DB.nin $DB.njs $DB.not $DB.nsq $DB.ntf $DB.nto ; rm -f $FRAG $TMP ; } echoxit() { @@ -311,7 +315,9 @@ fi echo " window size (-w) $WSIZE" ; ## formatting reference sequences ############################################################################ -$MKBNDB -in $REFFILE -out $DB &>/dev/null ; +echo ">ref" > $FRAG ; +grep -v "^>" $REFFILE >> $FRAG ; +$MKBNDB -in $FRAG -out $DB &>/dev/null ; ## checking no. threads ###################################################################################### [ $NTHREADS -lt 1 ] && NTHREADS=1; @@ -341,7 +347,7 @@ do INFO=$(randfile $FASTA.info); INFOFILES="$INFOFILES $INFO"; finalizer() { - rm -f $DB $DB.ndb $DB.nhr $DB.nin $DB.not $DB.nsq $DB.ntf $DB.nto ; + rm -f $DB $DB.ndb $DB.nhr $DB.nin $DB.njs $DB.not $DB.nsq $DB.ntf $DB.nto ; rm -f $FRAG $TMP ; rm -f $INFOFILES ; }