ó
‹™xTc           @   sr   d  d l  m Z d  d l m Z d  d l m Z d  d l m Z m Z d „  Z	 d „  Z
 e d k rn e
 ƒ  n  d S(	   iÿÿÿÿ(   t   division(   t   pdb(   t   flex(   t   Exception_expectedt	   show_diffc          C   sK  t  j d d  d t j d ƒ ƒ }  |  j ƒ  } | j d |  j ƒ  ƒ } | j | j	 ƒ  k s` t
 ‚ | j } | d ƒ j ƒ  d k s‡ t
 ‚ | d ƒ j ƒ  | j k s¨ t
 ‚ | d ƒ j ƒ  d k sÆ t
 ‚ | d	 d
 t ƒj ƒ  d k sê t
 ‚ | d	 d
 t ƒd  k st
 ‚ | d ƒ j ƒ  d k s&t
 ‚ | d ƒ j ƒ  | j k sGt
 ‚ t | d ƒ ƒ d d d d d g k stt
 ‚ t | d ƒ ƒ g  k s’t
 ‚ t | d ƒ ƒ d d d d d g k s¿t
 ‚ t | d ƒ ƒ d d d d d g k sìt
 ‚ t | d ƒ ƒ d d g k st
 ‚ x´ d d g D]¦ } t | d | ƒ ƒ d g k sHt
 ‚ t | d | ƒ ƒ d g k smt
 ‚ t | d | ƒ ƒ d d d  g k s˜t
 ‚ t | d! | ƒ ƒ d d d  g k st
 ‚ qWt | d" ƒ ƒ d# g k sèt
 ‚ t | d$ ƒ ƒ d% d& d' d( d) g k st
 ‚ t | d* ƒ ƒ d% d& d' d( d) g k sBt
 ‚ t | d+ ƒ ƒ d% d& d' d( d) g k sot
 ‚ t | d, ƒ ƒ d% d& d' d( d) g k sœt
 ‚ t | d- ƒ ƒ d. g k s½t
 ‚ t | d/ ƒ ƒ d0 d. g k sát
 ‚ t | d1 ƒ ƒ d0 d. g k st
 ‚ t | d2 ƒ ƒ d0 d. g k s)t
 ‚ t | d3 ƒ ƒ d0 d. g k sMt
 ‚ t | d4 ƒ ƒ d5 d6 d7 d8 d9 g k szt
 ‚ t | d: ƒ ƒ d; d< d= d> d? d@ dA dB g k s°t
 ‚ t | dC ƒ ƒ dB d' d( d) g k sÚt
 ‚ t | dD ƒ ƒ d% d& d' d( d) g k st
 ‚ t | dD ƒ ƒ d% d& d' d( d) g k s4t
 ‚ t | dE ƒ ƒ d d5 dF dG dH dI dJ dK dL dM dN dO g k svt
 ‚ t | dP ƒ ƒ t d d  ƒ dQ dR dS dT dU dV g k s³t
 ‚ t | dW ƒ ƒ dX dY g k s×t
 ‚ t | dZ ƒ ƒ t d[ ƒ t dJ d\ ƒ k st
 ‚ t | d] ƒ ƒ t d d ƒ dQ dR dS dT dU dV g k sEt
 ‚ t | d^ ƒ ƒ d% d& d' d( d) g k srt
 ‚ t | d_ ƒ ƒ d# g k s“t
 ‚ t | d` ƒ ƒ d% d& d' d( d) d# g k sÃt
 ‚ t | da ƒ ƒ d\ db g k sçt
 ‚ t | dc ƒ ƒ db dd g k st
 ‚ t | de ƒ ƒ t df ƒ k s/t
 ‚ t | dg ƒ ƒ d d dh di dj g k s\t
 ‚ t | dk ƒ ƒ dl dm dn do g k s†t
 ‚ t | dp ƒ ƒ d d dh di dj g k s³t
 ‚ t | dq ƒ ƒ dl dm dn do g k sÝt
 ‚ t | dr ƒ ƒ dH dl dm dn do g k s
t
 ‚ t | ds ƒ ƒ d; d& g k s.t
 ‚ t | dt ƒ ƒ d0 dm d. du g k sXt
 ‚ t | dv ƒ ƒ d% d& d' d( d) g k s…t
 ‚ t | dw ƒ ƒ d  d0 dx dy dm d. du g k s¸t
 ‚ t | dz ƒ ƒ dN dO d{ g k sßt
 ‚ t | d| ƒ ƒ dN dO g k s	t
 ‚ t | d} ƒ ƒ d{ g k s$	t
 ‚ t | d~ ƒ ƒ g  k sB	t
 ‚ t | d ƒ ƒ t d  ƒ dJ dK dQ dR dL dM dS dN dT dU dO dV g k sŽ	t
 ‚ t d d  ƒ dQ dR dS dN dT dU dO dV g } t | d€ ƒ ƒ | k s×	t
 ‚ t | d� ƒ ƒ | k sõ	t
 ‚ t | d‚ ƒ ƒ | k s
t
 ‚ t d d[ ƒ dQ dR dX dS dN dT dU dO dV dY d{ g } t | dƒ ƒ ƒ | k se
t
 ‚ d6 d0 d8 d9 g } t | d„ ƒ ƒ | k s•
t
 ‚ t | d… ƒ ƒ | k s³
t
 ‚ t | d† ƒ ƒ | k sÑ
t
 ‚ t | d‡ ƒ ƒ | k sï
t
 ‚ d5 d6 d0 g } t | dˆ ƒ ƒ | k st
 ‚ t | d‰ ƒ ƒ | k s:t
 ‚ t | dŠ ƒ ƒ | k sXt
 ‚ t | d‹ ƒ ƒ | k svt
 ‚ d0 d8 d9 g } t | dŒ ƒ ƒ | k s£t
 ‚ t | d� ƒ ƒ | k sÁt
 ‚ t | dŽ ƒ ƒ | k sßt
 ‚ t | d� ƒ ƒ | k sýt
 ‚ t | d� ƒ ƒ d0 g k st
 ‚ d5 d6 d0 d8 d9 g } t | d‘ ƒ ƒ | k sQt
 ‚ t | d’ ƒ ƒ | k sot
 ‚ t | d“ ƒ ƒ d# g k s�t
 ‚ t | d” ƒ ƒ d& d' d( d) g k sºt
 ‚ t | d• ƒ ƒ d% g k sÛt
 ‚ t | d– ƒ ƒ d5 dG g k sÿt
 ‚ t | d— ƒ ƒ t d[ ƒ k s#t
 ‚ t | d˜ ƒ ƒ d# dX dS dN dO dV dY d\ db dd d{ g k sbt
 ‚ t | d™ ƒ ƒ t | dš ƒ ƒ k oëd[ d; d< d= d> d d d? d d@ d dA d d› dœ d d� dž dB d dŸ d  d¡ d¢ d% d& d' d( d) d# g k n söt
 ‚ t | d£ ƒ ƒ d6 d¤ d¥ d  d0 dx dy dm d. du d d8 d d9 dR dS dU g k sGt
 ‚ y | d¦ ƒ Wn: t  j j k
 r‘} t | ƒ j d§ ƒ d k s˜t
 ‚ n Xt ‚ y | d¨ ƒ Wn: t  j j k
 râ} t | ƒ j d© ƒ d k sét
 ‚ n Xt ‚ y | dª ƒ Wn: t  j j k
 r3} t | ƒ j d« ƒ d k s:t
 ‚ n Xt ‚ xŠ d¬ d­ d® d¯ d° d± d² d³ d´ dµ g
 D]d } y | d¶ | ƒ WnD t  j j k
 r¼} t | ƒ j d· | ƒ sÃt
 t | ƒ ƒ ‚ q_Xt ‚ q_W| j d¸ d¹ ƒ } t | ƒ d5 k sñt
 ‚ t | d ƒ d5 d6 d0 d8 d9 g k st
 ‚ | j ƒ  }	 x+ | d D] }
 |	 |
 j d¹ k s3t
 ‚ q3W| j d¦ dº ƒ } t | ƒ d5 k s€t
 ‚ t | d ƒ dX dY g k s¢t
 ‚ |	 j ƒ  | j d» ƒ } | j t ƒ d5 k sÖt
 ‚ t  j j d¼ | d½ |	 ƒ } | j t ƒ d k st
 ‚ g  d¾ j ƒ  D] } t  j  j! d¿ | ƒ ^ q} d% g d& g g d% g d' g g g  g  g g } xQ t" | | ƒ D]@ \ } } g  | j# | ƒ D] } t | ƒ ^ q“| k swt
 ‚ qwWt  j d d  d t j dÀ ƒ ƒ j ƒ  } | j ƒ  } | j } t | dÁ ƒ ƒ d d5 dF g k st
 ‚ t | dÂ ƒ ƒ dG dH dI g k sEt
 ‚ t | dÃ ƒ ƒ t d ƒ k sit
 ‚ t | dÄ ƒ ƒ d5 dF dI g k s�t
 ‚ t | dÅ ƒ ƒ d5 dF dI g k s·t
 ‚ t | dÆ ƒ ƒ d5 dF dI g k sÞt
 ‚ t  j d d  d t j dÇ ƒ ƒ j ƒ  } | j ƒ  } | j } t | dÄ ƒ ƒ d5 dF dI g k sAt
 ‚ t | dÅ ƒ ƒ d5 dF dI g k sht
 ‚ t | dÆ ƒ ƒ d5 dF dI g k s�t
 ‚ t | dÁ ƒ ƒ d d5 dF g k s¶t
 ‚ t | dÂ ƒ ƒ g  k sÔt
 ‚ t | dÈ ƒ ƒ g  k sòt
 ‚ t | dÉ ƒ ƒ dG dH dI g k st
 ‚ t  j d d  d t j dÊ ƒ ƒ j ƒ  } | j ƒ  } | j } t | d— ƒ ƒ d d5 dI d d6 d¤ d¥ dl d7 g	 k sŽt
 ‚ xj dË dÌ g D]\ } y | | ƒ Wn? t  j j k
 rð} t$ t | ƒ dÍ | | f ƒ s÷t
 ‚ q›Xt ‚ q›Wy | dÎ ƒ Wn5 t  j j k
 r@} t$ t | ƒ dÏ ƒ sGt
 ‚ n Xt ‚ t  j d d  d t j dÐ ƒ ƒ j ƒ  } | j ƒ  } | j d˜ ƒ j ƒ  d k s›t
 ‚ t  j d d  d t j dÑ ƒ ƒ j ƒ  } | j ƒ  } | j dÒ ƒ } | j ƒ  dH k sõt
 ‚ | j dÓ ƒ } | j ƒ  d6 k st
 ‚ y | j dÔ ƒ } Wn t  j j k
 rJ} n Xt ‚ t  j d d  d t j dÕ ƒ ƒ j ƒ  } | j ƒ  } | j dÖ ƒ } t | ƒ dF d d6 d¤ d¥ g k sºt
 ‚ | j d× ƒ } t | ƒ d d6 d¤ d¥ g k sít
 ‚ t  j d d  d t j dØ ƒ ƒ j ƒ  } | j ƒ  } | j } t | dÙ ƒ ƒ d5 dF d d6 g k sSt
 ‚ t | dÚ ƒ ƒ dI d d6 g k szt
 ‚ t | dÛ ƒ ƒ dG dH g k sžt
 ‚ t | dÜ ƒ ƒ d d5 dF dG dH dI g k sÎt
 ‚ t | dÝ ƒ ƒ dG d d6 g k sõt
 ‚ t | dÞ ƒ ƒ d d6 g k st
 ‚ y | dß ƒ Wn t  j j k
 r@n Xt ‚ d  S(à   Nt   source_infot   linessŠ  CRYST1   50.800   50.800  155.300  90.00  90.00  90.00 P 43 21 2     8
MODEL        1
ATOM      4  N   SER     1       8.753  29.755  61.685  1.00 49.13
ATOM      5  CA  SER     1       9.242  30.200  62.974  1.00 46.62
ANISOU    5  CA  SER     1    343    490   2719    -45   -169    617
ATOM      6  C   SER     1      10.453  29.500  63.579  1.00 41.99
ATOM      7  O   SER     1      10.593  29.607  64.814  1.00 43.24
ANISOU    7  O   SER     1    343    490   2719    -45   -169    617
ATOM      8  CB  SER     1       8.052  30.189  63.974  1.00 53.00
ATOM      9  OG  SER     1       7.294  31.409  63.930  1.00 57.79
ATOM     10  N   ARG     2      11.360  28.819  62.827  1.00 36.48
ATOM     11  CA  ARG     2      12.548  28.316  63.532  1.00 30.20
ATOM     12  C   ARG     2      13.502  29.501  63.500  1.00 25.54
ATOM     13  O   ARG     2      13.730  30.037  62.407  1.00 23.86
ATOM     14  CB  ARG     2      13.241  27.119  62.861  1.00 27.44
ATOM     15  CG  ARG     2      12.412  25.849  62.964  1.00 23.66
ATOM     16  CD  ARG     2      13.267  24.651  63.266  1.00 23.98
ATOM     17  NE  ARG     2      13.948  24.115  62.135  1.00 22.71
ATOM     18  CZ  ARG     2      15.114  23.487  62.201  1.00 21.38
ATOM     19  NH1 ARG     2      15.845  23.331  63.301  1.00 19.34
ATOM     20  NH2 ARG     2      15.575  23.030  61.051  1.00 26.66
ATOM     21  N   PRO     3J     13.947  29.997  64.680  1.00 22.94
ATOM     22  CA  PRO     3J     14.902  31.100  64.827  1.00 20.19
ATOM     23  C   PRO     3J     16.195  30.718  64.086  1.00 18.44
ATOM     24  O   PRO     3J     16.545  29.521  64.086  1.00 19.76
ATOM     25  CB  PRO     3J     15.133  31.218  66.313  1.00 19.17
ATOM     26  CG  PRO     3J     14.065  30.364  66.951  1.00 15.12
ATOM     27  CD  PRO     3J     13.816  29.289  65.966  1.00 19.56
ATOM     28  N  AILE     4      16.953  31.648  63.512  1.00 15.29
ATOM     29  CA AILE     4      18.243  31.372  62.859  1.00 14.32
ATOM     30  C  AILE     4      19.233  32.112  63.743  1.00 13.54
ATOM     31  O  AILE     4      19.105  33.315  64.009  1.00 11.84
ATOM     32  CB AILE     4      18.298  31.951  61.406  1.00 13.62
ATOM     33  CG1AILE     4      17.157  31.300  60.620  1.00 18.39
ATOM     34  CG2AILE     4      19.661  31.747  60.743  1.00 13.64
ATOM     35  CD1AILE     4      16.879  32.102  59.355  1.00 16.69
ATOM     28  N  BILE     4      16.953  31.648  63.512  1.00 15.29
ATOM     29  CA BILE     4      18.243  31.372  62.859  1.00 14.32
ATOM     30  C  BILE     4      19.233  32.112  63.743  1.00 13.54
ATOM     31  O  BILE     4      19.105  33.315  64.009  1.00 11.84
ATOM     32  CB BILE     4      18.298  31.951  61.406  1.00 13.62
ATOM     33  CG1BILE     4      17.157  31.300  60.620  1.00 18.39
ATOM     34  CG2BILE     4      19.661  31.747  60.743  1.00 13.64
ATOM1200035  CD1BILE     4      16.879  32.102  59.355  1.00 16.69
TER      36      ILE     4
ENDMDL
MODEL        2
HETATM 1451  PA  5GP H 187      29.875  44.488  69.823  1.00 19.62
HETATM 1452  O1A 5GP H 187      28.526  44.888  69.143  1.00 19.86
HETATM 1453  O2A 5GP H 187      30.764  44.617  68.702  1.00 23.42
HETATM 1454  O3A 5GP H 187      30.319  45.004  71.073  1.00 20.20
HETATM 1455  O5* 5GP H 187      29.683  43.016  70.027  1.00 20.32
HETATM 1456  C5* 5GP H 187      30.740  42.297  70.837  1.00 21.47
HETATM 1457  C4* 5GP H 187      30.677  40.747  70.770  1.00 21.56
HETATM 1458  O4* 5GP H 187      29.608  40.160  71.599  1.00 20.50
HETATM 1459  C3* 5GP H 187      30.547  40.121  69.352  1.00 20.18
HETATM 1460  O3* 5GP H 187      31.228  38.864  69.416  1.00 23.65
HETATM 1461  C2* 5GP H 187      29.031  39.871  69.248  1.00 18.78
HETATM 1462  O2* 5GP H 187      28.685  38.690  68.496  1.00 20.45
HETATM 1463  C1* 5GP H 187      28.634  39.641  70.688  1.00 17.09
HETATM 1464  N9  5GP H 187      27.238  39.525  71.076  1.00 15.35
HETATM 1465  C8  5GP H 187      26.330  40.535  70.852  1.00 12.57
HETATM 1466  N7' 5GP H 187      25.175  40.314  71.417  1.00 12.88
HETATM 1467  C5  5GP H 187      25.278  39.082  72.070  1.00 10.75
HETATM 1468  C6  5GP H 187      24.326  38.354  72.827  1.00  9.77
HETATM 1469  O6  5GP H 187      23.169  38.678  73.029  1.00  8.66
HETATM 1470  N1' 5GP H 187      24.836  37.190  73.270  1.00  9.67
HETATM 1471  C2  5GP H 187      26.075  36.701  73.001  1.00  9.84
HETATM 1472  N2  5GP H 187      26.361  35.490  73.520  1.00  9.77
HETATM 1473  N3  5GP H 187      27.005  37.353  72.310  1.00 10.31
HETATM 1474  C4  5GP H 187      26.583  38.559  71.844  1.00 12.50
ENDMDL
MODEL        3
HETATM 1475  S   SO4 S 188      31.424  42.923  60.396  1.00 55.69           S4+
HETATM 1476  O1  SO4 S 188      31.631  41.513  60.336  1.00 59.84           o1-
HETATM 1477  O2  SO4 S 188      32.533  43.699  59.932  1.00 49.98           O1-
HETATM 1478  O3  SO4 S 188      31.128  43.217  61.738  1.00 59.44           O1-
HETATM 1479  O4  SO4 S 188      30.353  43.201  59.539  1.00 60.54           O1-
HETATM 1480  O   HOH W 200      29.478  23.354  61.364  1.00  8.67      WATE
ATOM   2000  A1  AAA X   1       8.753  29.755  61.685  1.00 49.13
ATOM   2001  A2  AAA X   1       9.242  30.200  62.974  1.00 46.62
ATOM   2002  A1  BBB X   2      11.360  28.819  62.827  1.00 36.48
ATOM   2003  A2  BBB X   2      12.548  28.316  63.532  1.00 30.20
ATOM   2004  A1  AAA Y   1       8.753  29.755  61.685  1.00 49.13
ATOM   2005  A2  AAA Y   1       9.242  30.200  62.974  1.00 46.62
ATOM   2006  A1  CCC Y   5       9.242  30.200  62.974  1.00 46.62
ATOM   2007  A2  BBB Y   2      12.548  28.316  63.532  1.00 30.20
ATOM   2008  A1  AAA Z   1K      8.753  29.755  61.685  1.00 49.13
ATOM   2009  A1  BBB Z   2      11.360  28.819  62.827  1.00 36.48
ATOM   2010  A2  BBB Z   2      12.548  28.316  63.532  1.00 30.20
ATOM   2011  A1  AAAZZ   1K      8.753  29.755  61.685  1.00 49.13
ATOM   2012  A1  BBBZZ   2      11.360  28.819  62.827  1.00 36.48
ATOM   2013  A1  CCCZZ   5       9.242  30.200  62.974  1.00 46.62
ATOM   2014  A1  CCCZZA001       9.242  30.200  62.974  1.00 46.62
ATOM   2015  A1  CCCZZA002       9.242  30.200  62.974  1.00 46.62
ATOM   2016  A1  CCCZZA003       9.242  30.200  62.974  1.00 46.62
ATOM   2017  A1  AAAUU  1K       8.753  29.755  61.685  1.00 49.13
ENDMDL
END
t   special_position_settingst    i    t   allt   nones   optional nonet   optionals   not alls   not nones	   name c?\*i-   i.   i0   i2   i4   s   name 'C?\*'s   name ' C?\*'s   name ' c?\*'s   name n?'i7   i;   s   and s   altloc a %sname ni   s   altloc b %sname ni    s   altloc ' ' %sname ni   i   s   altid ' ' %sname ns   resname hohiE   s   resname SO4i@   iA   iB   iC   iD   s   resname so4s   resname So4s   resname S?4s   resname pro and name cgi   s$   resname pro and (name cg or name ca)i   s$   resname pro AND (name cg or name ca)s$   resname pro and (name cg OR name ca)s$   resname pro AND (name cg OR name ca)s(   not resname pro and (name cg or name ca)i   i   i   i   i!   s   chain h and name o*i)   i*   i+   i,   i/   i1   i3   i:   s%   (chain h or chain s) and name o[2-46]s
   resseq 188s
   resseq 1:1i   i   i   i   iF   iG   iJ   iK   iN   iQ   s
   resseq 2:2iH   iI   iM   iO   iP   iR   s
   resseq 5:5iL   iS   s
   resseq 1:5i(   iT   s
   resseq 2:3s   resseq 188:188s   resseq 200:200s   resseq 188:200s   resseq 9999:A002iU   s   resseq A002:A003iV   s   resseq :iX   s   resseq :2 and name n*i   i   i   s   resseq 2: and name cbi
   i   i   i$   s   resseq 1:2 and name n*s   resseq 2:4 and name cbs   model 1 and name cbs   model 2:3 and name o1*s   icode j and name c?i   s	   resid 188s   resid 3Ji   i   s   resid 1KiW   s   resid '   1K's   resid '  1K 's   resi '  1K's	   resid 1:2s
   resid 1K:2s   resid '   1K:2's   resid '  1K:2's   resi '  1K:  1K 's   resseq 2:4 and name cas   resseq 2 : 4 and name cas   resseq 2: 4 and name cas   resseq 2 :4 and name cas   resseq :3 and name cas   resseq : 3 and name cas   name ca and resseq :3s   name ca and resseq : 3s   resseq 3: and name cas   resseq 3 : and name cas   name ca and resseq 3:s   name ca and resseq 3 :s   name ca and resseq 3s   resseq : and name cas   name ca and resseq :s
   segid wates	   element os	   charge 4+t   anisout   pepnamest   single_atom_residuet   heterot   hetatmi5   i6   i8   i9   i<   i=   i>   i?   s   within(3, resname PRO)i   i	   t   resseqs   Missing argument for resseq.s
   resseq 3:2s/   range with first index > last index: resseq 3:2s   resid ' 1K :2's1   range with first index > last index: resid  1K :2s   altloc a and and name ns   altloc a and or name ns   altloc a or and name ns   altloc a or or name ns
   and name ns	   or name nt   nots   not nots   altloc a optionals   optional optional altloc at   strings,   Atom selection string leading to error:
  %st   names    CA s      5s   resname 5GP and name C6t	   selectiont	   pdb_atomssÂ   LINK         S   SO4 S 188                 O1  SO4 S 188
LINK         S   SO4 S 188                 O2  SO4 S 188
LINK         NZ  LYS A 680        1.260    C4A PLP D   1                LYS-PLP
t   pdb_strs%  CRYST1   21.937    4.866   23.477  90.00 107.08  90.00 P 1 21 1      2
ATOM      2  CA  GLY A   1      -9.052   4.207   4.651  1.00 16.57           C
ATOM      6  CA  ASN A   2      -6.522   2.038   2.831  1.00 14.10           C
ATOM     14  CA  Asn A   3      -3.193   1.904   4.589  1.00 11.74           C
ATOM     22  CA  GLN a   4       0.384   1.888   3.199  1.00 10.53           C
ATOM     31  CA  GLN a   5       3.270   2.361   5.640  1.00 11.39           C
ATOM     40  CA  ASN a   6       6.831   2.310   4.318  1.00 12.30           C
END
s   chain As   chain as   name cas   resname asns   resname ASNs   resname Asns%  CRYST1   21.937    4.866   23.477  90.00 107.08  90.00 P 1 21 1      2
ATOM      2  CA  GLY A   1      -9.052   4.207   4.651  1.00 16.57           C
ATOM      6  CA  ASN A   2      -6.522   2.038   2.831  1.00 14.10           C
ATOM     14  CA  ASN A   3      -3.193   1.904   4.589  1.00 11.74           C
ATOM     22  CA  GLN b   4       0.384   1.888   3.199  1.00 10.53           C
ATOM     31  CA  GLN b   5       3.270   2.361   5.640  1.00 11.39           C
ATOM     40  CA  ASN b   6       6.831   2.310   4.318  1.00 12.30           C
END
s   chain Bs   chain bsg  ATOM      5  CA  SER     1       9.242  30.200  62.974  1.00 46.62
ATOM     11  CA  ARG     2      12.548  28.316  63.532  1.00 30.20
ATOM     21  N   NON     3J     13.947  29.997  64.680  1.00 22.94
ATOM     22  CA  NON     3J     14.902  31.100  64.827  1.00 20.19
ATOM     24  O   NON     3J     16.545  29.521  64.086  1.00 19.76
ATOM     28  N  AILE     4      16.953  31.648  63.512  1.00 15.29
ATOM     29  CA AILE     4      18.243  31.372  62.859  1.00 14.32
ATOM     30  C  AILE     4      19.233  32.112  63.743  1.00 13.54
ATOM     31  O  AILE     4      19.105  33.315  64.009  1.00 11.84
ATOM     41  N   PRO     5      13.947  29.997  64.680  1.00 22.94
ATOM     42  CA  PRO     5      14.902  31.100  64.827  1.00 20.19
ATOM     44  O   PRO     5      16.545  29.521  64.086  1.00 19.76
ATOM     45  CA  CA      6      16.545  29.521  64.086  1.00 19.76
t   peptidet   proteins…   Sorry: "%s" atom selection keyword not available:
  Please try using "pepnames" instead.
Atom selection string leading to error:
  %ss    chain A or (peptyde and name ca)s‡   RuntimeError: Atom selection syntax error at word "peptyde".
Atom selection string leading to error:
  chain A or (peptyde and name ca)sQ   ATOM    459  CA  SER A  58
ATOM    463  OG ASER A  58
ATOM    464  OG BSER A  58
s  CRYST1   21.937    4.866   23.477  90.00 107.08  90.00 P 1 21 1      2
ATOM      2  CA  GLY A 666      -9.052   4.207   4.651  1.00 16.57           C
ATOM      6  CA  ASN A 777      -6.522   2.038   2.831  1.00 14.10           C
ATOM     14  CA  ASN A   1      -3.193   1.904   4.589  1.00 11.74           C
ATOM     22  CA  GLN A   2       0.384   1.888   3.199  1.00 10.53           C
ATOM     31  CA  GLN A   3       3.270   2.361   5.640  1.00 11.39           C
ATOM     40  CA  ASN A   4       6.831   2.310   4.318  1.00 12.30           C
TER
ATOM      6  CA  ASN B 777      -6.522   2.038   2.831  1.00 14.10           C
ATOM     14  CA  ASN B   1      -3.193   1.904   4.589  1.00 11.74           C
ATOM     22  CA  GLN B   2       0.384   1.888   3.199  1.00 10.53           C
END
s   resid 777 through 3 and chain As   resid 777 through 3s   resid 777 through chain Asi  CRYST1   21.937    4.866   23.477  90.00 107.08  90.00 P 1 21 1      2
ATOM      6  CA  ASN B 777      -6.522   2.038   2.831  1.00 14.10           C
ATOM     14  CA  ASN B   1      -3.193   1.904   4.589  1.00 11.74           C
ATOM     22  CA  GLN B   2       0.384   1.888   3.199  1.00 10.53           C
TER
ATOM      2  CA  GLY A 666      -9.052   4.207   4.651  1.00 16.57           C
ATOM      6  CA  ASN A 777      -6.522   2.038   2.831  1.00 14.10           C
ATOM     14  CA  ASN A   1      -3.193   1.904   4.589  1.00 11.74           C
ATOM     22  CA  GLN A   2       0.384   1.888   3.199  1.00 10.53           C
ATOM     31  CA  GLN A   3       3.270   2.361   5.640  1.00 11.39           C
ATOM     40  CA AASN A   4       6.831   2.310   4.318  1.00 12.30           C
ATOM     41  CA BASN A   4       6.831   2.310   4.318  1.00 12.30           C
TER
END
s   resid 2 through 4s   chain A and resid 2 through 4s$  ATOM      5  CA  SER     1       9.242  30.200  62.974  1.00 46.62
ATOM     11  O   HOH     2      12.548  28.316  63.532  1.00 30.20
ATOM     21  O   WAT C   3J     13.947  29.997  64.680  1.00 22.94           O
ATOM     22  CA  NON     3J     14.902  31.100  64.827  0.50 20.19
ATOM     42  CA  PRO     5      14.902  31.100  64.827  1.00 20.19
ATOM     44  O   PRO     5      16.545  29.521  64.086  1.00 19.76
HETATM   45  O  ZHOH     6      16.545  29.521  64.086  0.00 19.76
HETATM   4   O  KHOH     6      16.545  29.521  64.086  0.00 19.76
t   waters   bfactor < 20s   name CA and bfactor = 20.19s   occupancy > 0s   occupancy < 1s   occupancy = 0s   within(5, resname PRO)(%   R   t   inputt   NoneR   t   split_linest   construct_hierarchyt   atom_selection_cacheR   t   n_seqt
   atoms_sizet   AssertionErrort
   iselectiont   sizet   Truet   Falset   listt   ranget   atom_selectiont   AtomSelectionErrort   strt   findR   t   endswitht
   get_labelst   lent   atomsR   t   reset_i_seqR   t   countt&   expand_selection_to_entire_atom_groupst
   splitlinest   recordst   linkt   zipt   link_iselectionsR   (   t   pdb_inpt	   hierarchyt	   sel_cachet   iselt   conjt   expectedt   et   st   selR0   t   i_seqt   sel_completeR   t   link_recordst   expected_resultst   link_recordt   sele(    (    sJ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/pdb/tst_atom_selection.pyt   exercise_selection   sÜ   d		!$!---$%%+/!----!$$$$6*--B=$1=-!0$$$-*-*-$*-3'$!=+4!!*!$$0&95'*    #
+"
	+%		''$'''		'''''	9
!'$		*'$0'$c           C   s   t  ƒ  d GHd  S(   Nt   OK(   RH   (    (    (    sJ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/pdb/tst_atom_selection.pyt   run±  s    t   __main__N(   t
   __future__R    t   iotbxR   t   cctbx.array_familyR   t   libtbx.test_utilsR   R   RH   RJ   t   __name__(    (    (    sJ   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/iotbx/pdb/tst_atom_selection.pyt   <module>   s   	ÿ ¬	