#!/bin/sh

# Second run of SAPI to refine/score possible sites
# The output is a set of peaks in the SAPIPKS file
# The list of sites after the KAR keyword was obtained from
# sapi_run1.exam

sapi \
HKLIN $CEXAM/rnase/rnase25 SAPIPKS $CCP4_SCR/sapi_run2.peaks <<eof
TITLE Test of sapi refining Hg sites in rnase 25 data
LABIN DF=FHG2DEL SIGDF=SDFHG2DEL
KAR
    C    0.94457    0.21078    0.88274    1
    C    0.37739    0.16940    0.88159    2
    C    0.59956    0.24052    0.88245    3
    C    0.37793    0.19901    0.12397    4
    C    0.44887    0.19629    0.62509    5
    C    0.80046    0.09821    0.81837    6
    C    0.19947    0.17463    0.87489    7
    C    0.57133    0.08057    0.95709    8
    C    0.88168    0.23685    0.88885    9
    C    0.34479    0.13627    0.01350   10
    C    0.87973    0.22975    0.90961   11
    C    0.97571    0.04062    0.61159   12
    C    0.84639    0.05830    0.01431   13
    C    0.55053    0.04797    0.67305   14
    C    0.40993    0.02790    0.38379   15
    C    0.09616    0.23974    0.83284   16
    C    0.52910    0.11167    0.59637   17
    C    0.66983    0.17143    0.63231   18
    C    0.36417    0.06948    0.39542   19
    C    0.90900    0.16416    0.88284   20
    C    0.62325    0.12709    0.08309   21
    C    0.54520    0.05597    0.24936   22
    C    0.90430    0.13145    0.12407   23
    C    0.44817    0.15435    0.85332   24
    C    0.80456    0.04025    0.86747   25
    C    0.11532    0.19172    0.10521   26
    C    0.19393    0.07721    0.78510   27
    C    0.58834    0.12542    0.45870   28
    C    0.55652    0.06065    0.10324   29
    C    0.79940    0.05451    0.31214   30
END
eof