cat << EOF > 17210.cya

info   := normal
name   := 17210
tf_type	:= 2
cut_tf := 1000.0
seed   := 35621
tolerance   := 0.03, 0.03, 0.4
cut_upl := 0.5

cyanalib
#dyanalib

read seq ./17210.seq

read prot ./17210_0204.prot

read upl 17210.upl

read aco 17210.aco

#distance select all
distance stat
distance check

nproc=4

structures select 1..20
structure violate structures=5 delete=100
calc_all 100 command=anneal steps=10000
overview name=17210_viol structures=20 range=8..110 full 
ramachandran file=rama_viol nobackground label

EOF

cyana << 17210.cya > 17210.out
17210
quit

grep Upper 17210.ovw | sort -n +10 -r > 17210.viol





