#!/bin/sh

# First run of SAPI to find list of peaks

sapi \
HKLIN $CEXAM/rnase/rnase25 SAPIPKS $CCP4_SCR/sapi_run1.peaks <<eof
TITLE Test of sapi locating Hg sites in rnase 25 data
LABIN DF=FHG2DEL SIGDF=SDFHG2DEL
END
eof

# The output is a set of peaks in the SAPIPKS file
# Run SAPI again to refine these sites and obtain scores
# using the KAR keyword