#!/bin/sh -l
#PBS -q bioinf
#PBS -l nodes=1:ppn=1
#PBS -l walltime=4:00:00
#PBS -N run_step2_kobas
#PBS -o run_step2_kobas.out
#PBS -e run_step2_kobas.err
#PBS -m bea
cd $PBS_O_WORKDIR
module load bioinfo python
#/group/bioinfo/apps/apps/kobas-2.1.0/scripts/annotate.py
module load bioinfo KOBAS
module load blast

annotate.py -i diffPeak_DiffBind_H3K9me2_all_genes.decre.fil5k.fa -t fasta:pro -s pvu -o diffPeak_DiffBind_H3K9me2_all_genes.decre.fil5k.anno1
