cat << EOF > at3g01050_viol.cya

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

cyanalib
#dyanalib

read seq ./at3g01050.seq

read prot ./at3g01050.prot

read upl at3g01050.upl

read aco at3g01050.aco

#distance select all
distance stat
distance check

nproc=4

structures select 1..20
structure violate structures=10 delete=100
calc_all 100 command=anneal steps=10000
overview name=at3g01050_viol structures=20 range=10..108 full 
ramachandran file=rama_viol nobackground label

EOF

cyana << at3g01050_viol.cya > at3g01050_viol.out
at3g01050_viol
quit







