diff --git a/detectors/cds/bin/do_exonerate.csh b/detectors/cds/bin/do_exonerate.csh index 316299e..9962ebb 100755 --- a/detectors/cds/bin/do_exonerate.csh +++ b/detectors/cds/bin/do_exonerate.csh @@ -110,8 +110,21 @@ set base = $OutDir/$GenoName.$ProtName # # Alias the genome filename to a short name # to circumvent a bug in exonerate -set ShortLink = $OutDir/genome.fasta -ln -s $GenoFile $ShortLink + +echo $GenoFile | grep '^/' > /dev/null +if ( $status == 1 ) then + set AbsGenoFile = `pwd`/$GenoFile + set DirGenoFile = `dirname $AbsGenoFile` + set DirGenoFile = `(cd $DirGenoFile;pwd)` + set AbsGenoFile = $DirGenoFile/`basename $AbsGenoFile` +else + set AbsGenoFile = $GenoFile +endif + +set ShortDir = `mktemp -d` +set ShortLink = $ShortDir/genome.fasta +ln -s $AbsGenoFile $ShortLink +Notify " Building input shortcut $AbsGenoFile --> $ShortLink" # # skip exonerate calculations if already done @@ -232,7 +245,8 @@ $AwkCmd -f $LIB_DIR/toEmbl.awk $base.iff |\ # end # -rm -f $ShortLink +rm -rf $ShortDir + Notify " output file: $base.res" (\rm -f ?_$$) >> /dev/null diff --git a/detectors/cds/bin/go_cds.csh b/detectors/cds/bin/go_cds.csh index 78f42dc..14b3591 100755 --- a/detectors/cds/bin/go_cds.csh +++ b/detectors/cds/bin/go_cds.csh @@ -63,6 +63,8 @@ foreach dir ("core" "shell" "dust") endif end +cp $temp/ $Genome.cds.fasta $Genome.cds.fasta + # # pass2: transsplicing # @@ -71,6 +73,8 @@ end # pass3: prokov # +$PROG_DIR/do_prokov.sh $Fasta $Genome.cds.fasta $temp + # # end : output on stdout #