Starting phenix.real_space_refine on Tue Sep 24 00:00:44 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8a8j_15231/09_2024/8a8j_15231.cif Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8a8j_15231/09_2024/8a8j_15231.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.1 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8a8j_15231/09_2024/8a8j_15231.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8a8j_15231/09_2024/8a8j_15231.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8a8j_15231/09_2024/8a8j_15231.cif" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/8a8j_15231/09_2024/8a8j_15231.cif" } resolution = 3.1 write_initial_geo_file = False refinement { macro_cycles = 10 } qi { qm_restraints { package { program = *test } } } Starting job =============================================================================== ------------------------------------------------------------------------------- Citation: ********* Afonine PV, Poon BK, Read RJ, Sobolev OV, Terwilliger TC, Urzhumtsev A, Adams PD. (2018) Real-space refinement in PHENIX for cryo-EM and crystallography. Acta Cryst. D74:531-544. Validating inputs Origin is already at (0, 0, 0), no shifts will be applied ------------------------------------------------------------------------------- Processing inputs ***************** Set random seed Set to: 0 Set model cs if undefined Decide on map wrapping Map wrapping is set to: False Normalize map: mean=0, sd=1 Input map: mean= -0.000 sd= 0.077 Set stop_for_unknowns flag Set to: True Assert model is a single copy model Assert all atoms have isotropic ADPs Construct map_model_manager Extract box with map and model Check model and map are aligned Set scattering table Set to: electron Number of scattering types: 6 Type Number sf(0) Gaussians P 44 5.49 5 Mg 2 5.21 5 S 6 5.16 5 C 3748 2.51 5 N 1173 2.21 5 O 1186 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 4 residue(s): 0.03s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/chem_data/mon_lib" Total number of atoms: 6159 Number of models: 1 Model: "" Number of chains: 6 Chain: "A" Number of atoms: 2656 Number of conformers: 1 Conformer: "" Number of residues, atoms: 341, 2656 Classifications: {'peptide': 341} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 18, 'TRANS': 322} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 4 Chain: "B" Number of atoms: 2660 Number of conformers: 2 Conformer: "A" Number of residues, atoms: 341, 2652 Classifications: {'peptide': 341} Incomplete info: {'truncation_to_alanine': 2} Link IDs: {'PTRANS': 18, 'TRANS': 322} Unresolved non-hydrogen bonds: 8 Unresolved non-hydrogen angles: 9 Unresolved non-hydrogen dihedrals: 7 Planarities with less than four sites: {'GLU:plan': 1} Unresolved non-hydrogen planarities: 4 Conformer: "B" Number of residues, atoms: 341, 2652 Classifications: {'peptide': 341} Incomplete info: {'truncation_to_alanine': 2} Link IDs: {'PTRANS': 18, 'TRANS': 322} Unresolved non-hydrogen bonds: 8 Unresolved non-hydrogen angles: 9 Unresolved non-hydrogen dihedrals: 7 Planarities with less than four sites: {'GLU:plan': 1} Unresolved non-hydrogen planarities: 4 bond proxies already assigned to first conformer: 2692 Chain: "X" Number of atoms: 391 Number of conformers: 1 Conformer: "" Number of residues, atoms: 19, 391 Classifications: {'DNA': 19} Link IDs: {'rna3p': 18} Chain: "Y" Number of atoms: 388 Number of conformers: 1 Conformer: "" Number of residues, atoms: 19, 388 Classifications: {'DNA': 19} Link IDs: {'rna3p': 18} Chain: "A" Number of atoms: 32 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 32 Unusual residues: {' MG': 1, 'ANP': 1} Classifications: {'undetermined': 2} Link IDs: {None: 1} Chain: "B" Number of atoms: 32 Number of conformers: 1 Conformer: "" Number of residues, atoms: 2, 32 Unusual residues: {' MG': 1, 'ANP': 1} Classifications: {'undetermined': 2} Link IDs: {None: 1} Time building chain proxies: 6.13, per 1000 atoms: 1.00 Number of scatterers: 6159 At special positions: 0 Unit cell: (77.4, 89.44, 94.6, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 6 Type Number sf(0) S 6 16.00 P 44 15.00 Mg 2 11.99 O 1186 8.00 N 1173 7.00 C 3748 6.00 sf(0) = scattering factor at diffraction angle 0. Number of disulfides: simple=0, symmetry=0 Automatic linking Parameters for automatic linking Linking & cutoffs Metal : Auto - 3.00 Amino acid : True - 1.90 Carbohydrate : True - 1.99 Ligands : True - 1.99 Small molecules : False - 1.98 Amino acid - RNA/DNA : False Number of custom bonds: simple=0, symmetry=0 Time building additional restraints: 2.24 Conformation dependent library (CDL) restraints added in 1.4 seconds 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1238 Finding SS restraints... Secondary structure from input PDB file: 22 helices and 7 sheets defined 46.9% alpha, 15.2% beta 12 base pairs and 36 stacking pairs defined. Time for finding SS restraints: 1.98 Creating SS restraints... Processing helix chain 'A' and resid 35 through 47 Processing helix chain 'A' and resid 100 through 107 Processing helix chain 'A' and resid 114 through 118 Processing helix chain 'A' and resid 124 through 137 removed outlier: 3.636A pdb=" N ARG A 128 " --> pdb=" O PRO A 124 " (cutoff:3.500A) removed outlier: 3.692A pdb=" N ARG A 129 " --> pdb=" O LYS A 125 " (cutoff:3.500A) Processing helix chain 'A' and resid 141 through 162 Processing helix chain 'A' and resid 168 through 202 removed outlier: 3.832A pdb=" N ARG A 173 " --> pdb=" O SER A 169 " (cutoff:3.500A) removed outlier: 4.371A pdb=" N GLU A 174 " --> pdb=" O ALA A 170 " (cutoff:3.500A) removed outlier: 3.601A pdb=" N LEU A 175 " --> pdb=" O TRP A 171 " (cutoff:3.500A) removed outlier: 3.520A pdb=" N ARG A 177 " --> pdb=" O ARG A 173 " (cutoff:3.500A) removed outlier: 3.879A pdb=" N TYR A 178 " --> pdb=" O GLU A 174 " (cutoff:3.500A) removed outlier: 3.875A pdb=" N ARG A 192 " --> pdb=" O ARG A 188 " (cutoff:3.500A) removed outlier: 3.733A pdb=" N ALA A 194 " --> pdb=" O LEU A 190 " (cutoff:3.500A) Proline residue: A 195 - end of helix Processing helix chain 'A' and resid 220 through 228 Processing helix chain 'A' and resid 228 through 236 Processing helix chain 'A' and resid 241 through 245 removed outlier: 3.880A pdb=" N ASP A 245 " --> pdb=" O PRO A 242 " (cutoff:3.500A) Processing helix chain 'A' and resid 261 through 285 Processing helix chain 'A' and resid 300 through 312 Processing helix chain 'B' and resid 35 through 47 removed outlier: 3.698A pdb=" N LEU B 47 " --> pdb=" O ILE B 43 " (cutoff:3.500A) Processing helix chain 'B' and resid 101 through 107 removed outlier: 3.936A pdb=" N GLU B 106 " --> pdb=" O THR B 103 " (cutoff:3.500A) Processing helix chain 'B' and resid 114 through 120 removed outlier: 3.926A pdb=" N ALA B 120 " --> pdb=" O ASP B 117 " (cutoff:3.500A) Processing helix chain 'B' and resid 124 through 137 removed outlier: 3.661A pdb=" N LEU B 132 " --> pdb=" O ARG B 128 " (cutoff:3.500A) removed outlier: 3.530A pdb=" N ASP B 133 " --> pdb=" O ARG B 129 " (cutoff:3.500A) removed outlier: 3.570A pdb=" N ARG B 134 " --> pdb=" O ALA B 130 " (cutoff:3.500A) Processing helix chain 'B' and resid 140 through 163 removed outlier: 4.382A pdb=" N ALA B 144 " --> pdb=" O SER B 140 " (cutoff:3.500A) removed outlier: 3.851A pdb=" N ALA B 145 " --> pdb=" O ARG B 141 " (cutoff:3.500A) removed outlier: 3.734A pdb=" N GLU B 151 " --> pdb=" O LEU B 147 " (cutoff:3.500A) removed outlier: 3.831A pdb=" N LYS B 152 " --> pdb=" O SER B 148 " (cutoff:3.500A) removed outlier: 3.604A pdb=" N ALA B 153 " --> pdb=" O ALA B 149 " (cutoff:3.500A) removed outlier: 3.805A pdb=" N ALA B 163 " --> pdb=" O ALA B 159 " (cutoff:3.500A) Processing helix chain 'B' and resid 168 through 205 removed outlier: 3.610A pdb=" N ASP B 180 " --> pdb=" O ALA B 176 " (cutoff:3.500A) removed outlier: 4.165A pdb=" N GLU B 181 " --> pdb=" O ARG B 177 " (cutoff:3.500A) removed outlier: 3.629A pdb=" N ARG B 192 " --> pdb=" O ARG B 188 " (cutoff:3.500A) Proline residue: B 195 - end of helix removed outlier: 4.250A pdb=" N ALA B 203 " --> pdb=" O GLU B 199 " (cutoff:3.500A) removed outlier: 3.783A pdb=" N LEU B 204 " --> pdb=" O VAL B 200 " (cutoff:3.500A) Processing helix chain 'B' and resid 220 through 226 Processing helix chain 'B' and resid 228 through 236 Processing helix chain 'B' and resid 241 through 245 removed outlier: 3.642A pdb=" N ASP B 245 " --> pdb=" O PRO B 242 " (cutoff:3.500A) Processing helix chain 'B' and resid 261 through 285 removed outlier: 3.925A pdb=" N THR B 267 " --> pdb=" O GLY B 263 " (cutoff:3.500A) removed outlier: 3.676A pdb=" N LEU B 268 " --> pdb=" O GLU B 264 " (cutoff:3.500A) removed outlier: 3.547A pdb=" N GLU B 282 " --> pdb=" O ARG B 278 " (cutoff:3.500A) Processing helix chain 'B' and resid 300 through 312 removed outlier: 3.599A pdb=" N GLN B 312 " --> pdb=" O LEU B 308 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 3 through 9 removed outlier: 3.623A pdb=" N LEU A 5 " --> pdb=" O GLU A 71 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'A' and resid 24 through 26 removed outlier: 6.066A pdb=" N GLY A 24 " --> pdb=" O ALA A 317 " (cutoff:3.500A) removed outlier: 7.825A pdb=" N LEU A 319 " --> pdb=" O GLY A 24 " (cutoff:3.500A) removed outlier: 6.304A pdb=" N SER A 26 " --> pdb=" O LEU A 319 " (cutoff:3.500A) removed outlier: 3.807A pdb=" N LEU A 290 " --> pdb=" O GLN A 316 " (cutoff:3.500A) removed outlier: 6.164A pdb=" N VAL A 111 " --> pdb=" O LEU A 291 " (cutoff:3.500A) removed outlier: 7.726A pdb=" N ASP A 293 " --> pdb=" O VAL A 111 " (cutoff:3.500A) removed outlier: 6.796A pdb=" N VAL A 113 " --> pdb=" O ASP A 293 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'A' and resid 209 through 214 removed outlier: 3.511A pdb=" N ARG A 254 " --> pdb=" O LEU A 251 " (cutoff:3.500A) Processing sheet with id=AA4, first strand: chain 'A' and resid 332 through 334 Processing sheet with id=AA5, first strand: chain 'B' and resid 3 through 9 removed outlier: 5.747A pdb=" N LEU B 3 " --> pdb=" O GLU B 73 " (cutoff:3.500A) removed outlier: 6.435A pdb=" N GLU B 73 " --> pdb=" O LEU B 3 " (cutoff:3.500A) removed outlier: 3.543A pdb=" N LEU B 5 " --> pdb=" O GLU B 71 " (cutoff:3.500A) removed outlier: 3.548A pdb=" N ALA B 66 " --> pdb=" O LEU B 85 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'B' and resid 110 through 113 removed outlier: 6.045A pdb=" N VAL B 111 " --> pdb=" O LEU B 291 " (cutoff:3.500A) removed outlier: 7.507A pdb=" N ASP B 293 " --> pdb=" O VAL B 111 " (cutoff:3.500A) removed outlier: 6.664A pdb=" N VAL B 113 " --> pdb=" O ASP B 293 " (cutoff:3.500A) removed outlier: 3.934A pdb=" N LEU B 290 " --> pdb=" O GLN B 316 " (cutoff:3.500A) removed outlier: 3.718A pdb=" N ALA B 317 " --> pdb=" O GLY B 24 " (cutoff:3.500A) removed outlier: 3.741A pdb=" N LEU B 319 " --> pdb=" O SER B 26 " (cutoff:3.500A) removed outlier: 8.977A pdb=" N CYS B 331 " --> pdb=" O LEU B 25 " (cutoff:3.500A) removed outlier: 7.003A pdb=" N ALA B 27 " --> pdb=" O CYS B 331 " (cutoff:3.500A) removed outlier: 7.806A pdb=" N VAL B 333 " --> pdb=" O ALA B 27 " (cutoff:3.500A) removed outlier: 6.736A pdb=" N VAL B 29 " --> pdb=" O VAL B 333 " (cutoff:3.500A) removed outlier: 3.655A pdb=" N SER B 332 " --> pdb=" O LEU B 339 " (cutoff:3.500A) Processing sheet with id=AA7, first strand: chain 'B' and resid 209 through 214 264 hydrogen bonds defined for protein. 713 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 30 hydrogen bonds 60 hydrogen bond angles 0 basepair planarities 12 basepair parallelities 36 stacking parallelities Total time for adding SS restraints: 1.94 Time building geometry restraints manager: 1.89 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.22 - 1.33: 1694 1.33 - 1.45: 1156 1.45 - 1.57: 3409 1.57 - 1.69: 84 1.69 - 1.81: 10 Bond restraints: 6353 Sorted by residual: bond pdb=" O3A ANP A 901 " pdb=" PB ANP A 901 " ideal model delta sigma weight residual 1.700 1.556 0.144 2.00e-02 2.50e+03 5.17e+01 bond pdb=" O3A ANP B 901 " pdb=" PB ANP B 901 " ideal model delta sigma weight residual 1.700 1.557 0.143 2.00e-02 2.50e+03 5.10e+01 bond pdb=" N3B ANP B 901 " pdb=" PG ANP B 901 " ideal model delta sigma weight residual 1.795 1.693 0.102 2.00e-02 2.50e+03 2.58e+01 bond pdb=" N3B ANP A 901 " pdb=" PG ANP A 901 " ideal model delta sigma weight residual 1.795 1.694 0.101 2.00e-02 2.50e+03 2.56e+01 bond pdb=" O5' ANP B 901 " pdb=" PA ANP B 901 " ideal model delta sigma weight residual 1.655 1.573 0.082 2.00e-02 2.50e+03 1.69e+01 ... (remaining 6348 not shown) Histogram of bond angle deviations from ideal: 0.00 - 2.87: 8687 2.87 - 5.74: 69 5.74 - 8.61: 11 8.61 - 11.47: 4 11.47 - 14.34: 2 Bond angle restraints: 8773 Sorted by residual: angle pdb=" PB ANP B 901 " pdb=" N3B ANP B 901 " pdb=" PG ANP B 901 " ideal model delta sigma weight residual 126.95 112.61 14.34 3.00e+00 1.11e-01 2.29e+01 angle pdb=" PB ANP A 901 " pdb=" N3B ANP A 901 " pdb=" PG ANP A 901 " ideal model delta sigma weight residual 126.95 112.74 14.21 3.00e+00 1.11e-01 2.24e+01 angle pdb=" O1B ANP A 901 " pdb=" PB ANP A 901 " pdb=" O2B ANP A 901 " ideal model delta sigma weight residual 120.08 109.17 10.91 3.00e+00 1.11e-01 1.32e+01 angle pdb=" O1B ANP B 901 " pdb=" PB ANP B 901 " pdb=" O2B ANP B 901 " ideal model delta sigma weight residual 120.08 109.51 10.57 3.00e+00 1.11e-01 1.24e+01 angle pdb=" N3B ANP A 901 " pdb=" PG ANP A 901 " pdb=" O2G ANP A 901 " ideal model delta sigma weight residual 101.05 110.86 -9.81 3.00e+00 1.11e-01 1.07e+01 ... (remaining 8768 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 17.62: 3096 17.62 - 35.24: 413 35.24 - 52.86: 204 52.86 - 70.48: 37 70.48 - 88.10: 14 Dihedral angle restraints: 3764 sinusoidal: 1830 harmonic: 1934 Sorted by residual: dihedral pdb=" CG ARG A 192 " pdb=" CD ARG A 192 " pdb=" NE ARG A 192 " pdb=" CZ ARG A 192 " ideal model delta sinusoidal sigma weight residual -180.00 -135.54 -44.46 2 1.50e+01 4.44e-03 1.05e+01 dihedral pdb=" CA GLN A 8 " pdb=" C GLN A 8 " pdb=" N ARG A 9 " pdb=" CA ARG A 9 " ideal model delta harmonic sigma weight residual 180.00 163.97 16.03 0 5.00e+00 4.00e-02 1.03e+01 dihedral pdb=" N ARG B 198 " pdb=" CA BARG B 198 " pdb=" CB BARG B 198 " pdb=" CG BARG B 198 " ideal model delta sinusoidal sigma weight residual -60.00 -114.43 54.43 3 1.50e+01 4.44e-03 9.28e+00 ... (remaining 3761 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.056: 868 0.056 - 0.112: 90 0.112 - 0.169: 13 0.169 - 0.225: 0 0.225 - 0.281: 2 Chirality restraints: 973 Sorted by residual: chirality pdb=" C3' ANP A 901 " pdb=" C2' ANP A 901 " pdb=" C4' ANP A 901 " pdb=" O3' ANP A 901 " both_signs ideal model delta sigma weight residual False -2.36 -2.64 0.28 2.00e-01 2.50e+01 1.98e+00 chirality pdb=" C3' ANP B 901 " pdb=" C2' ANP B 901 " pdb=" C4' ANP B 901 " pdb=" O3' ANP B 901 " both_signs ideal model delta sigma weight residual False -2.36 -2.63 0.28 2.00e-01 2.50e+01 1.90e+00 chirality pdb=" CA GLU A 294 " pdb=" N GLU A 294 " pdb=" C GLU A 294 " pdb=" CB GLU A 294 " both_signs ideal model delta sigma weight residual False 2.51 2.34 0.17 2.00e-01 2.50e+01 6.92e-01 ... (remaining 970 not shown) Planarity restraints: 1011 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C GLY B 86 " 0.039 5.00e-02 4.00e+02 5.78e-02 5.34e+00 pdb=" N PRO B 87 " -0.100 5.00e-02 4.00e+02 pdb=" CA PRO B 87 " 0.029 5.00e-02 4.00e+02 pdb=" CD PRO B 87 " 0.032 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C GLY A 241 " -0.024 5.00e-02 4.00e+02 3.62e-02 2.10e+00 pdb=" N PRO A 242 " 0.063 5.00e-02 4.00e+02 pdb=" CA PRO A 242 " -0.019 5.00e-02 4.00e+02 pdb=" CD PRO A 242 " -0.020 5.00e-02 4.00e+02 delta sigma weight rms_deltas residual plane pdb=" C SER B 114 " -0.022 5.00e-02 4.00e+02 3.29e-02 1.73e+00 pdb=" N PRO B 115 " 0.057 5.00e-02 4.00e+02 pdb=" CA PRO B 115 " -0.017 5.00e-02 4.00e+02 pdb=" CD PRO B 115 " -0.018 5.00e-02 4.00e+02 ... (remaining 1008 not shown) Histogram of nonbonded interaction distances: 2.01 - 2.59: 25 2.59 - 3.17: 4322 3.17 - 3.74: 8502 3.74 - 4.32: 11898 4.32 - 4.90: 20269 Nonbonded interactions: 45016 Sorted by model distance: nonbonded pdb=" O2G ANP A 901 " pdb="MG MG A 902 " model vdw 2.009 2.170 nonbonded pdb=" OG1 THR B 37 " pdb="MG MG B 902 " model vdw 2.032 2.170 nonbonded pdb=" O2G ANP B 901 " pdb="MG MG B 902 " model vdw 2.048 2.170 nonbonded pdb=" O1B ANP A 901 " pdb=" O2G ANP A 901 " model vdw 2.150 3.040 nonbonded pdb=" O1B ANP B 901 " pdb=" O2G ANP B 901 " model vdw 2.155 3.040 ... (remaining 45011 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model Time spend for trying shortcut: 0.00 Found NCS groups: ncs_group { reference = (chain 'A' and (resid 1 through 74 or (resid 75 and (name N or name CA or name C \ or name O or name CB )) or resid 76 through 197 or resid 199 through 341 or res \ id 901 through 902)) selection = (chain 'B' and (resid 1 through 197 or resid 199 through 341 or resid 901 throug \ h 902)) } Set up NCS constraints No NCS constraints will be used in refinement. Set refine NCS operators Adjust number of macro_cycles Number of macro_cycles: 10 Reset NCS operators Extract rigid body selections Check and reset occupancies Occupancies: min=0.50 max=1.00 mean=1.00 Load rotamer database and sin/cos tables Set ADP refinement strategy ADPs will be refined as individual isotropic Make a string to write initial .geo file Internal consistency checks Time: Set random seed: 0.000 Set model cs if undefined: 0.000 Decide on map wrapping: 0.000 Normalize map: mean=0, sd=1: 2.650 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.000 Construct map_model_manager: 0.010 Extract box with map and model: 0.250 Check model and map are aligned: 0.050 Set scattering table: 0.070 Process input model: 22.250 Find NCS groups from input model: 0.220 Set up NCS constraints: 0.030 Set refine NCS operators: 0.000 Adjust number of macro_cycles: 0.000 Reset NCS operators: 0.000 Extract rigid body selections: 0.000 Check and reset occupancies: 0.000 Load rotamer database and sin/cos tables:2.920 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 28.450 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8039 moved from start: 0.0000 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.144 6353 Z= 0.283 Angle : 0.706 14.342 8773 Z= 0.316 Chirality : 0.039 0.281 973 Planarity : 0.005 0.058 1011 Dihedral : 19.597 88.099 2526 Min Nonbonded Distance : 2.009 Molprobity Statistics. All-atom Clashscore : 4.14 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.42 % Favored : 95.58 % Rotamer: Outliers : 0.19 % Allowed : 0.76 % Favored : 99.05 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.13 (0.33), residues: 679 helix: 1.46 (0.32), residues: 275 sheet: 0.34 (0.45), residues: 138 loop : -1.55 (0.36), residues: 266 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.006 0.001 TRP A 105 HIS 0.005 0.001 HIS A 277 PHE 0.005 0.001 PHE A 259 TYR 0.009 0.001 TYR B 150 ARG 0.009 0.000 ARG B 223 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 53 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 53 time to evaluate : 0.668 Fit side-chains outliers start: 0 outliers final: 0 residues processed: 53 average time/residue: 0.3871 time to fit residues: 24.2085 Evaluate side-chains 46 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 46 time to evaluate : 0.596 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 70 random chunks: chunk 59 optimal weight: 0.8980 chunk 53 optimal weight: 0.3980 chunk 29 optimal weight: 0.3980 chunk 18 optimal weight: 1.9990 chunk 35 optimal weight: 0.7980 chunk 28 optimal weight: 3.9990 chunk 54 optimal weight: 0.8980 chunk 21 optimal weight: 0.0970 chunk 33 optimal weight: 0.7980 chunk 40 optimal weight: 0.6980 chunk 63 optimal weight: 2.9990 overall best weight: 0.4778 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 83 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** B 10 ASN Total number of N/Q/H flips: 1 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7991 moved from start: 0.0980 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.033 6353 Z= 0.158 Angle : 0.527 5.759 8773 Z= 0.272 Chirality : 0.036 0.130 973 Planarity : 0.004 0.032 1011 Dihedral : 20.202 89.475 1243 Min Nonbonded Distance : 2.108 Molprobity Statistics. All-atom Clashscore : 3.98 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.10 % Favored : 96.90 % Rotamer: Outliers : 0.57 % Allowed : 4.39 % Favored : 95.04 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.53 (0.33), residues: 679 helix: 1.78 (0.32), residues: 277 sheet: 0.97 (0.47), residues: 128 loop : -1.52 (0.36), residues: 274 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.005 0.001 TRP B 171 HIS 0.004 0.001 HIS A 69 PHE 0.006 0.001 PHE A 139 TYR 0.008 0.001 TYR B 150 ARG 0.006 0.000 ARG B 2 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 51 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 49 time to evaluate : 0.683 Fit side-chains outliers start: 2 outliers final: 2 residues processed: 51 average time/residue: 0.3409 time to fit residues: 21.1141 Evaluate side-chains 47 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 45 time to evaluate : 0.645 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 308 LEU Chi-restraints excluded: chain B residue 116 LEU Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 70 random chunks: chunk 35 optimal weight: 2.9990 chunk 19 optimal weight: 0.8980 chunk 52 optimal weight: 1.9990 chunk 43 optimal weight: 0.6980 chunk 17 optimal weight: 2.9990 chunk 63 optimal weight: 0.8980 chunk 68 optimal weight: 5.9990 chunk 56 optimal weight: 4.9990 chunk 21 optimal weight: 0.9990 chunk 51 optimal weight: 0.7980 chunk 62 optimal weight: 0.0870 overall best weight: 0.6758 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 83 HIS B 10 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8004 moved from start: 0.1162 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.026 6353 Z= 0.173 Angle : 0.512 5.477 8773 Z= 0.262 Chirality : 0.037 0.148 973 Planarity : 0.004 0.033 1011 Dihedral : 19.739 86.634 1243 Min Nonbonded Distance : 2.013 Molprobity Statistics. All-atom Clashscore : 3.56 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.39 % Favored : 96.61 % Rotamer: Outliers : 0.95 % Allowed : 7.63 % Favored : 91.41 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.65 (0.33), residues: 679 helix: 1.86 (0.32), residues: 279 sheet: 1.12 (0.46), residues: 128 loop : -1.51 (0.36), residues: 272 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP B 171 HIS 0.005 0.001 HIS A 277 PHE 0.004 0.001 PHE A 139 TYR 0.008 0.001 TYR B 150 ARG 0.007 0.000 ARG B 223 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 50 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 46 time to evaluate : 0.659 Fit side-chains REVERT: B 229 ARG cc_start: 0.7577 (mtm110) cc_final: 0.7108 (mtt180) outliers start: 4 outliers final: 4 residues processed: 49 average time/residue: 0.3451 time to fit residues: 20.3371 Evaluate side-chains 47 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 43 time to evaluate : 0.604 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 308 LEU Chi-restraints excluded: chain B residue 101 LEU Chi-restraints excluded: chain B residue 116 LEU Chi-restraints excluded: chain B residue 216 THR Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 70 random chunks: chunk 47 optimal weight: 1.9990 chunk 33 optimal weight: 0.0470 chunk 7 optimal weight: 5.9990 chunk 30 optimal weight: 0.9980 chunk 42 optimal weight: 0.6980 chunk 63 optimal weight: 3.9990 chunk 67 optimal weight: 2.9990 chunk 60 optimal weight: 0.8980 chunk 18 optimal weight: 1.9990 chunk 56 optimal weight: 6.9990 chunk 38 optimal weight: 0.6980 overall best weight: 0.6678 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... No N/Q/H corrections needed this macrocycle ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8000 moved from start: 0.1260 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.003 0.022 6353 Z= 0.169 Angle : 0.496 5.435 8773 Z= 0.254 Chirality : 0.036 0.137 973 Planarity : 0.004 0.033 1011 Dihedral : 19.609 87.876 1243 Min Nonbonded Distance : 2.025 Molprobity Statistics. All-atom Clashscore : 3.48 Ramachandran Plot: Outliers : 0.00 % Allowed : 3.69 % Favored : 96.31 % Rotamer: Outliers : 1.15 % Allowed : 9.16 % Favored : 89.69 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.75 (0.33), residues: 679 helix: 1.93 (0.32), residues: 279 sheet: 1.11 (0.45), residues: 134 loop : -1.46 (0.36), residues: 266 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP B 171 HIS 0.005 0.001 HIS A 277 PHE 0.004 0.001 PHE A 139 TYR 0.008 0.001 TYR B 150 ARG 0.007 0.000 ARG B 223 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 50 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 45 time to evaluate : 0.676 Fit side-chains REVERT: B 339 LEU cc_start: 0.7662 (pt) cc_final: 0.7221 (tp) outliers start: 5 outliers final: 4 residues processed: 49 average time/residue: 0.3259 time to fit residues: 19.3761 Evaluate side-chains 45 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 41 time to evaluate : 0.661 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 98 ARG Chi-restraints excluded: chain A residue 308 LEU Chi-restraints excluded: chain B residue 116 LEU Chi-restraints excluded: chain B residue 216 THR Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 70 random chunks: chunk 0 optimal weight: 3.9990 chunk 50 optimal weight: 2.9990 chunk 27 optimal weight: 0.8980 chunk 57 optimal weight: 5.9990 chunk 46 optimal weight: 3.9990 chunk 34 optimal weight: 3.9990 chunk 60 optimal weight: 1.9990 chunk 17 optimal weight: 4.9990 chunk 22 optimal weight: 4.9990 chunk 13 optimal weight: 0.8980 chunk 39 optimal weight: 0.0010 overall best weight: 1.3590 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: A 83 HIS B 10 ASN Total number of N/Q/H flips: 2 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8063 moved from start: 0.1208 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.034 6353 Z= 0.276 Angle : 0.582 5.984 8773 Z= 0.293 Chirality : 0.039 0.138 973 Planarity : 0.005 0.036 1011 Dihedral : 19.711 89.689 1243 Min Nonbonded Distance : 1.989 Molprobity Statistics. All-atom Clashscore : 3.31 Ramachandran Plot: Outliers : 0.00 % Allowed : 4.72 % Favored : 95.28 % Rotamer: Outliers : 1.15 % Allowed : 10.31 % Favored : 88.55 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 0.00 % Cis-general : 0.00 % Twisted Proline : 0.00 % Twisted General : 0.00 % Rama-Z (Ramachandran plot Z-score): Interpretation: bad |Rama-Z| > 3; suspicious 2 < |Rama-Z| < 3; good |Rama-Z| < 2. Scores for whole/helix/sheet/loop are scaled independently; therefore, the values are not related in a simple manner. whole: 0.50 (0.32), residues: 679 helix: 1.72 (0.32), residues: 277 sheet: 1.29 (0.46), residues: 124 loop : -1.66 (0.34), residues: 278 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.004 0.001 TRP A 105 HIS 0.008 0.002 HIS A 277 PHE 0.007 0.001 PHE B 193 TYR 0.011 0.001 TYR B 150 ARG 0.008 0.000 ARG B 223 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1358 Ramachandran restraints generated. 679 Oldfield, 0 Emsley, 679 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 48 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 43 time to evaluate : 0.554 Fit side-chains REVERT: B 122 LEU cc_start: 0.8847 (OUTLIER) cc_final: 0.8425 (pp) outliers start: 5 outliers final: 4 residues processed: 47 average time/residue: 0.3429 time to fit residues: 19.3197 Evaluate side-chains 47 residues out of total 525 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 42 time to evaluate : 0.591 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 98 ARG Chi-restraints excluded: chain A residue 308 LEU Chi-restraints excluded: chain B residue 116 LEU Chi-restraints excluded: chain B residue 122 LEU Chi-restraints excluded: chain B residue 216 THR Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in run_program(real_space_refine.Program) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/iotbx/cli_parser.py", line 946, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/phenix/phenix/programs/real_space_refine.py", line 200, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 306, in __init__ self.caller(self.refine_xyz) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 702, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 766, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 51, in __init__ gradients_method = gradients_method) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/mmtbx/refinement/real_space/weight.py", line 87, in __init__ rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 368, in refine rms_angles_limit = rms_angles_limit) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 112, in refine states_collector = self.states_accumulator) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 174, in __init__ exception_handling_params=lbfgs_exception_handling_params) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 131, in run_c_plus_plus f, g = target_evaluator.compute_functional_and_gradients() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/cctbx/maptbx/real_space_refinement_simple.py", line 245, in compute_functional_and_gradients compute_gradients=True) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1479, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1381, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5461/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1262, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 50.0418 > 50: distance: 33 - 175: 20.818 distance: 44 - 188: 25.521 distance: 75 - 177: 24.596 distance: 86 - 180: 5.627 distance: 102 - 117: 10.690 distance: 106 - 129: 18.686 distance: 114 - 117: 8.910 distance: 117 - 118: 8.811 distance: 118 - 119: 8.156 distance: 118 - 121: 9.949 distance: 119 - 120: 8.880 distance: 119 - 129: 13.620 distance: 121 - 122: 5.400 distance: 122 - 124: 4.460 distance: 126 - 127: 3.943 distance: 129 - 130: 31.377 distance: 130 - 131: 16.765 distance: 130 - 133: 30.406 distance: 131 - 132: 19.169 distance: 131 - 134: 30.609 distance: 134 - 135: 28.245 distance: 135 - 136: 21.377 distance: 135 - 138: 28.720 distance: 136 - 137: 25.185 distance: 136 - 143: 20.559 distance: 138 - 139: 13.363 distance: 139 - 140: 10.372 distance: 140 - 141: 7.252 distance: 140 - 142: 11.755 distance: 143 - 144: 26.082 distance: 144 - 145: 33.604 distance: 144 - 147: 50.042 distance: 145 - 146: 36.891 distance: 145 - 148: 30.645 distance: 148 - 149: 22.799 distance: 149 - 150: 27.347 distance: 149 - 152: 25.498 distance: 150 - 151: 10.054 distance: 150 - 156: 27.958 distance: 152 - 153: 18.040 distance: 153 - 154: 19.471 distance: 153 - 155: 24.067 distance: 156 - 157: 24.673 distance: 156 - 162: 17.551 distance: 157 - 158: 20.782 distance: 157 - 160: 32.989 distance: 158 - 159: 6.170 distance: 158 - 163: 8.598 distance: 160 - 161: 13.500 distance: 161 - 162: 16.142 distance: 164 - 165: 9.017 distance: 165 - 166: 6.770 distance: 165 - 172: 6.943 distance: 167 - 168: 8.641 distance: 168 - 169: 10.778 distance: 169 - 171: 4.974 distance: 172 - 173: 12.218 distance: 173 - 174: 3.918 distance: 173 - 176: 22.427 distance: 174 - 175: 22.653 distance: 174 - 177: 26.840 distance: 177 - 178: 24.194 distance: 178 - 179: 20.247 distance: 178 - 181: 17.350 distance: 179 - 180: 16.991 distance: 179 - 185: 18.154 distance: 181 - 182: 28.323 distance: 181 - 183: 10.213 distance: 182 - 184: 29.824 distance: 185 - 186: 21.914 distance: 186 - 187: 21.233 distance: 186 - 189: 29.742 distance: 187 - 188: 15.574 distance: 187 - 193: 14.413 distance: 189 - 190: 15.391 distance: 190 - 191: 17.529 distance: 190 - 192: 19.320