$wrkdir (reference) in projects: dports

Project(s)

Full Search
Definition
Symbol
File Path
History
Type

Searched refs:wrkdir (Results 1 – 25 of 98) sorted by relevance

12 3 4

/dports/biology/infernal/infernal-1.1.3/hmmer/profmark/
H A D	x-psiblast	19 $wrkdir = shift; 37 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 51 $output = `esl-reformat -o $wrkdir/$msaname.pbl psiblast $wrkdir/$msaname.sto`; 60 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 66 $output = `esl-reformat -o $wrkdir/$msaname.fa fasta $wrkdir/$msaname.sto`; 69 $output = `seg $wrkdir/$msaname.fa -x > $wrkdir/$msaname.x.fa`; 75 …$output = `$blastpgp $blastopts1 -d $wrkdir/$msaname.x.fa -i $wrkdir/$msaname.query.fa -C $wrkdir/… 93 unlink "$wrkdir/$msaname.hmm"; 94 unlink "$wrkdir/$msaname.pbl"; 95 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-iterate-jackhmmer	10 $wrkdir = shift; 27 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 43 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 46 …$cmd = "$hmmbuild --hand $wrkdir/$msaname.hmm $wrkdir/$msaname.jck > /dev/null`"; … 47 …$cmd = "$hmmsearch $searchopts --tblout $wrkdir/$msaname.tmp $wrkdir/$msaname.hmm $fafile > /dev/n… 49 …if (! open(OUTPUT, "$wrkdir/$msaname.tmp")) { print "FAILED: to open $wrkdir/$msaname.tmp"; next M… 60 unlink "$wrkdir/$msaname.tmp"; 61 unlink "$wrkdir/$msaname.hmm"; 62 unlink "$wrkdir/$msaname.jck"; 63 unlink "$wrkdir/$msaname.query"; [all …]
H A D	x-sam	17 $wrkdir = shift; 31 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 40 $output = `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 44 $output = `esl-reformat -o $wrkdir/$msaname.a2m a2m $wrkdir/$msaname.sto`; 48 `(cd $wrkdir; $w05 $msaname.a2m $msaname.mod &> /dev/null)`; 73 unlink "$wrkdir/$msaname.1.$msaname.mod.dist"; 74 unlink "$wrkdir/$msaname.dist"; 75 unlink "$wrkdir/$msaname.mlib"; 76 unlink "$wrkdir/$msaname.mod"; 77 unlink "$wrkdir/$msaname.a2m"; [all …]
H A D	x-fps-phmmer	18 $wrkdir = shift; 31 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 44 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 48 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 57 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 58 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 60 …`$phmmer --cpu 1 --tblout $wrkdir/$msaname.tmp $opts $wrkdir/$msaname.query.fa $fafile > /dev/null… 63 …open(OUTPUT, "$wrkdir/$msaname.tmp") \|\| die "FAILED: to open $wrkdir/$msaname.tmp tabular output f… 87 unlink "$wrkdir/$msaname.tmp"; 88 unlink "$wrkdir/$msaname.query.fa"; [all …]
H A D	x-single-phmmer	5 $wrkdir = shift; 17 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 26 …$cmd = "esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname"; … 27 …$cmd = "esl-seqstat --amino -a $wrkdir/$msaname.sto \| grep '^=' \| awk '{print \$2}'"; … 31 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 33 …$cmd = "$phmmer $searchopts --tblout $wrkdir/$msaname.tmp $wrkdir/$msaname.query $fafile > /dev/nu… 35 …if (! open(OUTPUT, "$wrkdir/$msaname.tmp")) { print "FAILED: to open $wrkdir/$msaname.tmp"; next M… 47 unlink "$wrkdir/$msaname.tmp"; 48 unlink "$wrkdir/$msaname.query"; 49 unlink "$wrkdir/$msaname.sto";
H A D	x-fps-fasta	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 59 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 60 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 62 if (! open(FASTA, "$fasta $opts $wrkdir/$msaname.query.fa $fafile \|")) { 63 print "FAILED: $fasta $opts $wrkdir/$msaname.query.fa $fafile \|\n"; next MSA; 92 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-fps-ncbiblast	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 59 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 60 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 62 if (! open(BLASTP, "$blastp -d $fafile -i $wrkdir/$msaname.query.fa $blastopts \|")) { 63 print "FAILED: $blastp -d $fafile -i $wrkdir/$msaname.query.fa $blastopts \|\n"; next MSA; 92 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-fps-ncbiblast+	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 47 if ($?) { print "FAILED: esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname\n"; next MSA; } 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 59 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 60 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 62 if (! open(BLASTP, "$blastp -db $fafile -query $wrkdir/$msaname.query.fa $blastopts \|")) { 92 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-fps-ssearch	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 59 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 60 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 62 if (! open(FASTA, "$fasta $opts $wrkdir/$msaname.query.fa $fafile \|")) { 63 print "FAILED: $fasta $opts $wrkdir/$msaname.query.fa $fafile \|\n"; next MSA; 92 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-fps-wublast	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 60 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 61 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 63 if (! open(BLASTP, "$blastp $fafile $wrkdir/$msaname.query.fa $blastopts \|")) { 64 print "FAILED: $blastp $fafile $wrkdir/$msaname.query.fa $blastopts \|\n"; next QUERY; 93 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-iterate-psiblast	6 $wrkdir = shift; 39 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 48 …$cmd = "esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname"; … 49 …$cmd = "esl-seqstat --amino -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'"; … 53 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 55 …$cmd = "$blastpgp $blastopts1 -d $fafile.iter -i $wrkdir/$msaname.query -C $wrkdir/$msaname.asnt"… 56 …$cmd = "$blastpgp $blastopts2 -d $fafile -R $wrkdir/$msaname.asnt 2>/dev/null \|"; … 69 unlink "$wrkdir/$msaname.query"; 70 unlink "$wrkdir/$msaname.sto"; 71 unlink "$wrkdir/$msaname.asnt";
H A D	x-single-ncbiblast	6 $wrkdir = shift; 24 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 33 …$cmd = "esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname"; $output … 34 …$cmd = "esl-seqstat --amino -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'"; $output … 38 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 40 …$cmd = "$blastp -d $fafile -i $wrkdir/$msaname.query $blastopts \|"; if (! open(BLASTP… 53 unlink "$wrkdir/$msaname.sto"; 54 unlink "$wrkdir/$msaname.query";
H A D	x-single-ssearch	6 $wrkdir = shift; 25 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 34 …$cmd = "esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname"; $output … 35 …$cmd = "esl-seqstat --amino -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'"; $output … 39 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 42 …$cmd = "$fasta $opts $wrkdir/$msaname.query $fafile \|"; if (! open(FASTA… 53 unlink "$wrkdir/$msaname.query"; 54 unlink "$wrkdir/$msaname.sto";
/dports/biology/hmmer/hmmer-3.3/profmark/
H A D	x-psiblast	19 $wrkdir = shift; 37 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 51 $output = `esl-reformat -o $wrkdir/$msaname.pbl psiblast $wrkdir/$msaname.sto`; 60 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 66 $output = `esl-reformat -o $wrkdir/$msaname.fa fasta $wrkdir/$msaname.sto`; 69 $output = `seg $wrkdir/$msaname.fa -x > $wrkdir/$msaname.x.fa`; 75 …$output = `$blastpgp $blastopts1 -d $wrkdir/$msaname.x.fa -i $wrkdir/$msaname.query.fa -C $wrkdir/… 93 unlink "$wrkdir/$msaname.hmm"; 94 unlink "$wrkdir/$msaname.pbl"; 95 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-iterate-jackhmmer	10 $wrkdir = shift; 27 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 43 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 46 …$cmd = "$hmmbuild --hand $wrkdir/$msaname.hmm $wrkdir/$msaname.jck > /dev/null`"; … 47 …$cmd = "$hmmsearch $searchopts --tblout $wrkdir/$msaname.tmp $wrkdir/$msaname.hmm $fafile > /dev/n… 49 …if (! open(OUTPUT, "$wrkdir/$msaname.tmp")) { print "FAILED: to open $wrkdir/$msaname.tmp"; next M… 60 unlink "$wrkdir/$msaname.tmp"; 61 unlink "$wrkdir/$msaname.hmm"; 62 unlink "$wrkdir/$msaname.jck"; 63 unlink "$wrkdir/$msaname.query"; [all …]
H A D	x-sam	17 $wrkdir = shift; 31 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 40 $output = `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 44 $output = `esl-reformat -o $wrkdir/$msaname.a2m a2m $wrkdir/$msaname.sto`; 48 `(cd $wrkdir; $w05 $msaname.a2m $msaname.mod &> /dev/null)`; 73 unlink "$wrkdir/$msaname.1.$msaname.mod.dist"; 74 unlink "$wrkdir/$msaname.dist"; 75 unlink "$wrkdir/$msaname.mlib"; 76 unlink "$wrkdir/$msaname.mod"; 77 unlink "$wrkdir/$msaname.a2m"; [all …]
H A D	x-fps-phmmer	18 $wrkdir = shift; 31 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 44 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 48 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 57 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 58 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 60 …`$phmmer --cpu 1 --tblout $wrkdir/$msaname.tmp $opts $wrkdir/$msaname.query.fa $fafile > /dev/null… 63 …open(OUTPUT, "$wrkdir/$msaname.tmp") \|\| die "FAILED: to open $wrkdir/$msaname.tmp tabular output f… 87 unlink "$wrkdir/$msaname.tmp"; 88 unlink "$wrkdir/$msaname.query.fa"; [all …]
H A D	x-single-phmmer	5 $wrkdir = shift; 17 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 26 …$cmd = "esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname"; … 27 …$cmd = "esl-seqstat --amino -a $wrkdir/$msaname.sto \| grep '^=' \| awk '{print \$2}'"; … 31 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 33 …$cmd = "$phmmer $searchopts --tblout $wrkdir/$msaname.tmp $wrkdir/$msaname.query $fafile > /dev/nu… 35 …if (! open(OUTPUT, "$wrkdir/$msaname.tmp")) { print "FAILED: to open $wrkdir/$msaname.tmp"; next M… 47 unlink "$wrkdir/$msaname.tmp"; 48 unlink "$wrkdir/$msaname.query"; 49 unlink "$wrkdir/$msaname.sto";
H A D	x-fps-fasta	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 59 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 60 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 62 if (! open(FASTA, "$fasta $opts $wrkdir/$msaname.query.fa $fafile \|")) { 63 print "FAILED: $fasta $opts $wrkdir/$msaname.query.fa $fafile \|\n"; next MSA; 92 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-fps-ncbiblast	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 59 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 60 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 62 if (! open(BLASTP, "$blastp -d $fafile -i $wrkdir/$msaname.query.fa $blastopts \|")) { 63 print "FAILED: $blastp -d $fafile -i $wrkdir/$msaname.query.fa $blastopts \|\n"; next MSA; 92 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-fps-ncbiblast+	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 47 if ($?) { print "FAILED: esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname\n"; next MSA; } 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 59 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 60 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 62 if (! open(BLASTP, "$blastp -db $fafile -query $wrkdir/$msaname.query.fa $blastopts \|")) { 92 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-fps-ssearch	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 59 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 60 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 62 if (! open(FASTA, "$fasta $opts $wrkdir/$msaname.query.fa $fafile \|")) { 63 print "FAILED: $fasta $opts $wrkdir/$msaname.query.fa $fafile \|\n"; next MSA; 92 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-fps-wublast	18 $wrkdir = shift; 33 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 46 `esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname`; 50 $output = `esl-seqstat -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'`; 51 if ($?) { print "FAILED: esl-seqstat -a $wrkdir/$msaname.sto\n"; next MSA; } 60 $output = `esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname`; 61 …if ($?) { print "FAILED: esl-sfetch -o $wrkdir/$msaname.query.fa $wrkdir/$msaname.sto $qname\n"; n… 63 if (! open(BLASTP, "$blastp $fafile $wrkdir/$msaname.query.fa $blastopts \|")) { 64 print "FAILED: $blastp $fafile $wrkdir/$msaname.query.fa $blastopts \|\n"; next QUERY; 93 unlink "$wrkdir/$msaname.sto"; [all …]
H A D	x-iterate-psiblast	6 $wrkdir = shift; 39 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 48 …$cmd = "esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname"; … 49 …$cmd = "esl-seqstat --amino -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'"; … 53 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 55 …$cmd = "$blastpgp $blastopts1 -d $fafile.iter -i $wrkdir/$msaname.query -C $wrkdir/$msaname.asnt"… 56 …$cmd = "$blastpgp $blastopts2 -d $fafile -R $wrkdir/$msaname.asnt 2>/dev/null \|"; … 69 unlink "$wrkdir/$msaname.query"; 70 unlink "$wrkdir/$msaname.sto"; 71 unlink "$wrkdir/$msaname.asnt";
H A D	x-single-ncbiblast	6 $wrkdir = shift; 24 if (! -e $wrkdir) { die "$wrkdir doesn't exist"; } 33 …$cmd = "esl-afetch -o $wrkdir/$msaname.sto $msafile $msaname"; $output … 34 …$cmd = "esl-seqstat --amino -a $wrkdir/$msaname.sto \| grep "^=" \| awk '{print \$2}'"; $output … 38 …$cmd = "esl-sfetch -o $wrkdir/$msaname.query $wrkdir/$msaname.sto $qname > /dev/null"; … 40 …$cmd = "$blastp -d $fafile -i $wrkdir/$msaname.query $blastopts \|"; if (! open(BLASTP… 53 unlink "$wrkdir/$msaname.sto"; 54 unlink "$wrkdir/$msaname.query";

12 3 4