Updated 2021-02-10

Hmmer pbs

#PBS -N testHMMER
#PBS -A [Account]
#PBS -l nodes=1:ppn=2
#PBS -l pmem=2gb
#PBS -l walltime=5:00
#PBS -q inferno
#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