#!/bin/sh
#pid=13265
pid=$$
blastout=mqblast.out
listfile=mdpin.list
resultfile=mdp.result
mdpout=mdp.mergeout
restfile=mdp.rest
outfile=mdp.out

seqdbfile=newprotseq
qseqfile=

PATH=${PATH}:/dbb/project/MBGD/package/mqblast

mqblast -TMPDIR=/scratch/uchiyama/blasttmp \
	-OUTDIR2=/tmp/uchiyama/mqblast.$pid/odir \
	-OUTFILE="$blastout" \
	-mode=qsub -CopyAll -retainTmp \
	-PROCNUM=24 \
	$seqdbfile $qseqfile

create_list $blastout > $listfile
execmdp.pl $listfile $seqdbfile $outfile
update_list.pl -listfile=$listfile -outfiles=$outfile \
	-resultout=$resultfile -nonuniq > $restfile
