Updated 2019-06-24

hmmer.pbs

#PBS -N testHMMER
#PBS -l nodes=1:ppn=2
#PBS -l pmem=2gb
#PBS -l walltime=5:00
#PBS -q force-6
#PBS -j oe
#PBS -o testHMMER.out

cd $PBS_O_WORKDIR/tutorial
module load hmmer/3.1b1

# build a profile with hmmbuild
hmmbuild globins4.hmm globins4.sto

# search the sequence database with hmmsearch
hmmsearch globins4.hmm globins45.fa > globins4.out