Starting phenix.real_space_refine on Sat Jan 18 03:56:10 2025 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/cci-nas-00/data/ceres_data/7xd1_33132/01_2025/7xd1_33132.cif Found real_map, /net/cci-nas-00/data/ceres_data/7xd1_33132/01_2025/7xd1_33132.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=3.2 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/cci-nas-00/data/ceres_data/7xd1_33132/01_2025/7xd1_33132.map" default_real_map = "/net/cci-nas-00/data/ceres_data/7xd1_33132/01_2025/7xd1_33132.map" model { file = "/net/cci-nas-00/data/ceres_data/7xd1_33132/01_2025/7xd1_33132.cif" } default_model = "/net/cci-nas-00/data/ceres_data/7xd1_33132/01_2025/7xd1_33132.cif" } resolution = 3.2 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.003 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: 5 Type Number sf(0) Gaussians P 294 5.49 5 S 14 5.16 5 C 6648 2.51 5 N 2268 2.21 5 O 2824 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped. Time to flip 5 residue(s): 0.02s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/chem_data/mon_lib" Total number of atoms: 12048 Number of models: 1 Model: "" Number of chains: 10 Chain: "A" Number of atoms: 810 Number of conformers: 1 Conformer: "" Number of residues, atoms: 98, 810 Classifications: {'peptide': 98} Link IDs: {'PTRANS': 4, 'TRANS': 93} Chain: "B" Number of atoms: 641 Number of conformers: 1 Conformer: "" Number of residues, atoms: 80, 641 Classifications: {'peptide': 80} Link IDs: {'PTRANS': 1, 'TRANS': 78} Chain: "C" Number of atoms: 836 Number of conformers: 1 Conformer: "" Number of residues, atoms: 109, 836 Classifications: {'peptide': 109} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 5, 'TRANS': 103} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 4 Chain: "E" Number of atoms: 801 Number of conformers: 1 Conformer: "" Number of residues, atoms: 97, 801 Classifications: {'peptide': 97} Link IDs: {'PTRANS': 3, 'TRANS': 93} Chain: "F" Number of atoms: 627 Number of conformers: 1 Conformer: "" Number of residues, atoms: 79, 627 Classifications: {'peptide': 79} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 1, 'TRANS': 77} Unresolved non-hydrogen bonds: 6 Unresolved non-hydrogen angles: 7 Unresolved non-hydrogen dihedrals: 5 Planarities with less than four sites: {'ARG:plan': 1} Unresolved non-hydrogen planarities: 5 Chain: "G" Number of atoms: 836 Number of conformers: 1 Conformer: "" Number of residues, atoms: 109, 836 Classifications: {'peptide': 109} Incomplete info: {'truncation_to_alanine': 1} Link IDs: {'PTRANS': 5, 'TRANS': 103} Unresolved non-hydrogen bonds: 4 Unresolved non-hydrogen angles: 4 Unresolved non-hydrogen dihedrals: 4 Chain: "I" Number of atoms: 3031 Number of conformers: 1 Conformer: "" Number of residues, atoms: 147, 3031 Classifications: {'DNA': 147} Link IDs: {'rna3p': 146} Chain: "J" Number of atoms: 2996 Number of conformers: 1 Conformer: "" Number of residues, atoms: 147, 2996 Classifications: {'DNA': 147} Link IDs: {'rna3p': 146} Chain: "D" Number of atoms: 735 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 735 Classifications: {'peptide': 94} Link IDs: {'PTRANS': 2, 'TRANS': 91} Chain: "H" Number of atoms: 735 Number of conformers: 1 Conformer: "" Number of residues, atoms: 94, 735 Classifications: {'peptide': 94} Link IDs: {'PTRANS': 2, 'TRANS': 91} Time building chain proxies: 6.94, per 1000 atoms: 0.58 Number of scatterers: 12048 At special positions: 0 Unit cell: (119.214, 127.806, 91.29, 90, 90, 90) Space group: P 1 (No. 1) Number of sites at special positions: 0 Number of scattering types: 5 Type Number sf(0) S 14 16.00 P 294 15.00 O 2824 8.00 N 2268 7.00 C 6648 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.96 Conformation dependent library (CDL) restraints added in 979.8 milliseconds 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 1420 Finding SS restraints... Secondary structure from input PDB file: 35 helices and 10 sheets defined 68.8% alpha, 3.2% beta 129 base pairs and 276 stacking pairs defined. Time for finding SS restraints: 6.93 Creating SS restraints... Processing helix chain 'A' and resid 44 through 55 removed outlier: 3.927A pdb=" N TYR A 54 " --> pdb=" O GLU A 50 " (cutoff:3.500A) Processing helix chain 'A' and resid 63 through 76 Processing helix chain 'A' and resid 85 through 114 removed outlier: 3.525A pdb=" N ALA A 98 " --> pdb=" O GLU A 94 " (cutoff:3.500A) removed outlier: 3.531A pdb=" N LEU A 109 " --> pdb=" O GLU A 105 " (cutoff:3.500A) Processing helix chain 'A' and resid 120 through 131 removed outlier: 3.847A pdb=" N ILE A 124 " --> pdb=" O MET A 120 " (cutoff:3.500A) removed outlier: 3.755A pdb=" N GLN A 125 " --> pdb=" O PRO A 121 " (cutoff:3.500A) removed outlier: 3.538A pdb=" N ARG A 131 " --> pdb=" O ALA A 127 " (cutoff:3.500A) Processing helix chain 'B' and resid 25 through 29 Processing helix chain 'B' and resid 30 through 41 removed outlier: 3.581A pdb=" N GLY B 41 " --> pdb=" O LEU B 37 " (cutoff:3.500A) Processing helix chain 'B' and resid 49 through 76 removed outlier: 4.132A pdb=" N GLU B 53 " --> pdb=" O LEU B 49 " (cutoff:3.500A) removed outlier: 4.137A pdb=" N ARG B 67 " --> pdb=" O GLU B 63 " (cutoff:3.500A) removed outlier: 3.563A pdb=" N ASP B 68 " --> pdb=" O ASN B 64 " (cutoff:3.500A) removed outlier: 3.646A pdb=" N THR B 71 " --> pdb=" O ARG B 67 " (cutoff:3.500A) removed outlier: 3.677A pdb=" N ALA B 76 " --> pdb=" O TYR B 72 " (cutoff:3.500A) Processing helix chain 'B' and resid 82 through 92 removed outlier: 3.637A pdb=" N VAL B 86 " --> pdb=" O THR B 82 " (cutoff:3.500A) Processing helix chain 'C' and resid 16 through 21 removed outlier: 3.641A pdb=" N ALA C 21 " --> pdb=" O ARG C 17 " (cutoff:3.500A) Processing helix chain 'C' and resid 26 through 36 removed outlier: 3.737A pdb=" N VAL C 30 " --> pdb=" O PRO C 26 " (cutoff:3.500A) Processing helix chain 'C' and resid 46 through 72 removed outlier: 3.535A pdb=" N ALA C 53 " --> pdb=" O VAL C 49 " (cutoff:3.500A) Processing helix chain 'C' and resid 79 through 90 removed outlier: 3.646A pdb=" N LEU C 83 " --> pdb=" O ILE C 79 " (cutoff:3.500A) Processing helix chain 'C' and resid 90 through 97 removed outlier: 3.639A pdb=" N LYS C 95 " --> pdb=" O GLU C 91 " (cutoff:3.500A) Processing helix chain 'E' and resid 45 through 55 Processing helix chain 'E' and resid 63 through 75 removed outlier: 3.588A pdb=" N PHE E 67 " --> pdb=" O ARG E 63 " (cutoff:3.500A) removed outlier: 3.610A pdb=" N LEU E 70 " --> pdb=" O PRO E 66 " (cutoff:3.500A) removed outlier: 3.526A pdb=" N GLU E 73 " --> pdb=" O ARG E 69 " (cutoff:3.500A) Processing helix chain 'E' and resid 76 through 78 No H-bonds generated for 'chain 'E' and resid 76 through 78' Processing helix chain 'E' and resid 85 through 114 removed outlier: 3.565A pdb=" N ALA E 111 " --> pdb=" O THR E 107 " (cutoff:3.500A) Processing helix chain 'E' and resid 120 through 131 removed outlier: 3.845A pdb=" N ILE E 124 " --> pdb=" O MET E 120 " (cutoff:3.500A) Processing helix chain 'F' and resid 25 through 29 Processing helix chain 'F' and resid 30 through 42 Processing helix chain 'F' and resid 50 through 76 removed outlier: 3.999A pdb=" N ARG F 55 " --> pdb=" O TYR F 51 " (cutoff:3.500A) removed outlier: 3.960A pdb=" N GLY F 56 " --> pdb=" O GLU F 52 " (cutoff:3.500A) removed outlier: 3.531A pdb=" N LYS F 59 " --> pdb=" O ARG F 55 " (cutoff:3.500A) removed outlier: 3.848A pdb=" N ARG F 67 " --> pdb=" O GLU F 63 " (cutoff:3.500A) removed outlier: 3.756A pdb=" N ASP F 68 " --> pdb=" O ASN F 64 " (cutoff:3.500A) removed outlier: 3.789A pdb=" N THR F 71 " --> pdb=" O ARG F 67 " (cutoff:3.500A) removed outlier: 3.551A pdb=" N HIS F 75 " --> pdb=" O THR F 71 " (cutoff:3.500A) Processing helix chain 'F' and resid 82 through 92 removed outlier: 3.749A pdb=" N VAL F 86 " --> pdb=" O THR F 82 " (cutoff:3.500A) removed outlier: 3.557A pdb=" N ARG F 92 " --> pdb=" O TYR F 88 " (cutoff:3.500A) Processing helix chain 'G' and resid 16 through 21 removed outlier: 3.641A pdb=" N ALA G 21 " --> pdb=" O ARG G 17 " (cutoff:3.500A) Processing helix chain 'G' and resid 26 through 36 removed outlier: 3.737A pdb=" N VAL G 30 " --> pdb=" O PRO G 26 " (cutoff:3.500A) Processing helix chain 'G' and resid 46 through 72 removed outlier: 3.535A pdb=" N ALA G 53 " --> pdb=" O VAL G 49 " (cutoff:3.500A) Processing helix chain 'G' and resid 79 through 90 removed outlier: 3.646A pdb=" N LEU G 83 " --> pdb=" O ILE G 79 " (cutoff:3.500A) Processing helix chain 'G' and resid 90 through 97 removed outlier: 3.639A pdb=" N LYS G 95 " --> pdb=" O GLU G 91 " (cutoff:3.500A) Processing helix chain 'D' and resid 37 through 49 removed outlier: 3.572A pdb=" N LYS D 46 " --> pdb=" O TYR D 42 " (cutoff:3.500A) Processing helix chain 'D' and resid 58 through 84 removed outlier: 3.640A pdb=" N ARG D 79 " --> pdb=" O GLY D 75 " (cutoff:3.500A) Processing helix chain 'D' and resid 91 through 102 Processing helix chain 'D' and resid 104 through 124 removed outlier: 3.580A pdb=" N LYS D 108 " --> pdb=" O GLY D 104 " (cutoff:3.500A) removed outlier: 3.599A pdb=" N ALA D 110 " --> pdb=" O LEU D 106 " (cutoff:3.500A) removed outlier: 3.541A pdb=" N SER D 112 " --> pdb=" O LYS D 108 " (cutoff:3.500A) removed outlier: 3.663A pdb=" N ALA D 117 " --> pdb=" O GLU D 113 " (cutoff:3.500A) Processing helix chain 'H' and resid 37 through 49 removed outlier: 3.572A pdb=" N LYS H 46 " --> pdb=" O TYR H 42 " (cutoff:3.500A) Processing helix chain 'H' and resid 58 through 84 removed outlier: 3.640A pdb=" N ARG H 79 " --> pdb=" O GLY H 75 " (cutoff:3.500A) Processing helix chain 'H' and resid 91 through 102 Processing helix chain 'H' and resid 104 through 124 removed outlier: 3.580A pdb=" N LYS H 108 " --> pdb=" O GLY H 104 " (cutoff:3.500A) removed outlier: 3.599A pdb=" N ALA H 110 " --> pdb=" O LEU H 106 " (cutoff:3.500A) removed outlier: 3.541A pdb=" N SER H 112 " --> pdb=" O LYS H 108 " (cutoff:3.500A) removed outlier: 3.663A pdb=" N ALA H 117 " --> pdb=" O GLU H 113 " (cutoff:3.500A) Processing sheet with id=AA1, first strand: chain 'A' and resid 83 through 84 removed outlier: 6.907A pdb=" N ARG A 83 " --> pdb=" O VAL B 81 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA1 Processing sheet with id=AA2, first strand: chain 'A' and resid 118 through 119 Processing sheet with id=AA3, first strand: chain 'B' and resid 96 through 98 removed outlier: 6.597A pdb=" N THR B 96 " --> pdb=" O THR G 101 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA3 Processing sheet with id=AA4, first strand: chain 'C' and resid 42 through 43 removed outlier: 7.291A pdb=" N ARG C 42 " --> pdb=" O ILE D 89 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA4 Processing sheet with id=AA5, first strand: chain 'C' and resid 77 through 78 Processing sheet with id=AA6, first strand: chain 'C' and resid 100 through 102 Processing sheet with id=AA7, first strand: chain 'E' and resid 83 through 84 removed outlier: 7.013A pdb=" N ARG E 83 " --> pdb=" O VAL F 81 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA7 Processing sheet with id=AA8, first strand: chain 'E' and resid 118 through 119 Processing sheet with id=AA9, first strand: chain 'G' and resid 42 through 43 removed outlier: 7.220A pdb=" N ARG G 42 " --> pdb=" O ILE H 89 " (cutoff:3.500A) No H-bonds generated for sheet with id=AA9 Processing sheet with id=AB1, first strand: chain 'G' and resid 77 through 78 343 hydrogen bonds defined for protein. 1017 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 317 hydrogen bonds 630 hydrogen bond angles 0 basepair planarities 129 basepair parallelities 276 stacking parallelities Total time for adding SS restraints: 4.65 Time building geometry restraints manager: 3.34 seconds NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Histogram of bond lengths: 1.21 - 1.33: 2152 1.33 - 1.45: 4325 1.45 - 1.57: 5772 1.57 - 1.69: 588 1.69 - 1.81: 24 Bond restraints: 12861 Sorted by residual: bond pdb=" CB GLU D 113 " pdb=" CG GLU D 113 " ideal model delta sigma weight residual 1.520 1.598 -0.078 3.00e-02 1.11e+03 6.78e+00 bond pdb=" CB GLU H 113 " pdb=" CG GLU H 113 " ideal model delta sigma weight residual 1.520 1.597 -0.077 3.00e-02 1.11e+03 6.66e+00 bond pdb=" C3' DT I -59 " pdb=" O3' DT I -59 " ideal model delta sigma weight residual 1.422 1.492 -0.070 3.00e-02 1.11e+03 5.42e+00 bond pdb=" C3' DG I 20 " pdb=" O3' DG I 20 " ideal model delta sigma weight residual 1.422 1.488 -0.066 3.00e-02 1.11e+03 4.90e+00 bond pdb=" CG1 ILE B 50 " pdb=" CD1 ILE B 50 " ideal model delta sigma weight residual 1.513 1.435 0.078 3.90e-02 6.57e+02 4.04e+00 ... (remaining 12856 not shown) Histogram of bond angle deviations from ideal: 0.00 - 2.23: 18104 2.23 - 4.46: 478 4.46 - 6.70: 28 6.70 - 8.93: 12 8.93 - 11.16: 7 Bond angle restraints: 18629 Sorted by residual: angle pdb=" CB GLU D 113 " pdb=" CG GLU D 113 " pdb=" CD GLU D 113 " ideal model delta sigma weight residual 112.60 123.17 -10.57 1.70e+00 3.46e-01 3.87e+01 angle pdb=" CB GLU H 113 " pdb=" CG GLU H 113 " pdb=" CD GLU H 113 " ideal model delta sigma weight residual 112.60 123.16 -10.56 1.70e+00 3.46e-01 3.86e+01 angle pdb=" CA GLU H 113 " pdb=" CB GLU H 113 " pdb=" CG GLU H 113 " ideal model delta sigma weight residual 114.10 125.26 -11.16 2.00e+00 2.50e-01 3.11e+01 angle pdb=" CA GLU D 113 " pdb=" CB GLU D 113 " pdb=" CG GLU D 113 " ideal model delta sigma weight residual 114.10 125.19 -11.09 2.00e+00 2.50e-01 3.08e+01 angle pdb=" C THR E 80 " pdb=" N ASP E 81 " pdb=" CA ASP E 81 " ideal model delta sigma weight residual 121.54 131.03 -9.49 1.91e+00 2.74e-01 2.47e+01 ... (remaining 18624 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.23: 5612 35.23 - 70.45: 1353 70.45 - 105.68: 8 105.68 - 140.90: 2 140.90 - 176.13: 2 Dihedral angle restraints: 6977 sinusoidal: 4775 harmonic: 2202 Sorted by residual: dihedral pdb=" CA GLU A 50 " pdb=" C GLU A 50 " pdb=" N ILE A 51 " pdb=" CA ILE A 51 " ideal model delta harmonic sigma weight residual 180.00 153.63 26.37 0 5.00e+00 4.00e-02 2.78e+01 dihedral pdb=" CA PRO E 121 " pdb=" C PRO E 121 " pdb=" N LYS E 122 " pdb=" CA LYS E 122 " ideal model delta harmonic sigma weight residual 180.00 158.63 21.37 0 5.00e+00 4.00e-02 1.83e+01 dihedral pdb=" CA ILE F 29 " pdb=" C ILE F 29 " pdb=" N THR F 30 " pdb=" CA THR F 30 " ideal model delta harmonic sigma weight residual 180.00 -159.33 -20.67 0 5.00e+00 4.00e-02 1.71e+01 ... (remaining 6974 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.051: 1523 0.051 - 0.102: 513 0.102 - 0.153: 71 0.153 - 0.204: 9 0.204 - 0.255: 3 Chirality restraints: 2119 Sorted by residual: chirality pdb=" CA ASP E 81 " pdb=" N ASP E 81 " pdb=" C ASP E 81 " pdb=" CB ASP E 81 " both_signs ideal model delta sigma weight residual False 2.51 2.25 0.26 2.00e-01 2.50e+01 1.63e+00 chirality pdb=" CA ASP A 81 " pdb=" N ASP A 81 " pdb=" C ASP A 81 " pdb=" CB ASP A 81 " both_signs ideal model delta sigma weight residual False 2.51 2.27 0.24 2.00e-01 2.50e+01 1.49e+00 chirality pdb=" P DC I -58 " pdb=" OP1 DC I -58 " pdb=" OP2 DC I -58 " pdb=" O5' DC I -58 " both_signs ideal model delta sigma weight residual True 2.35 -2.56 -0.22 2.00e-01 2.50e+01 1.16e+00 ... (remaining 2116 not shown) Planarity restraints: 1340 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" CB TYR F 88 " 0.017 2.00e-02 2.50e+03 1.83e-02 6.68e+00 pdb=" CG TYR F 88 " -0.045 2.00e-02 2.50e+03 pdb=" CD1 TYR F 88 " 0.017 2.00e-02 2.50e+03 pdb=" CD2 TYR F 88 " 0.008 2.00e-02 2.50e+03 pdb=" CE1 TYR F 88 " -0.004 2.00e-02 2.50e+03 pdb=" CE2 TYR F 88 " 0.004 2.00e-02 2.50e+03 pdb=" CZ TYR F 88 " 0.002 2.00e-02 2.50e+03 pdb=" OH TYR F 88 " 0.000 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB TYR B 51 " -0.016 2.00e-02 2.50e+03 1.54e-02 4.72e+00 pdb=" CG TYR B 51 " 0.038 2.00e-02 2.50e+03 pdb=" CD1 TYR B 51 " -0.008 2.00e-02 2.50e+03 pdb=" CD2 TYR B 51 " -0.004 2.00e-02 2.50e+03 pdb=" CE1 TYR B 51 " -0.003 2.00e-02 2.50e+03 pdb=" CE2 TYR B 51 " -0.007 2.00e-02 2.50e+03 pdb=" CZ TYR B 51 " -0.002 2.00e-02 2.50e+03 pdb=" OH TYR B 51 " 0.003 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" CB PHE C 25 " -0.012 2.00e-02 2.50e+03 1.46e-02 3.71e+00 pdb=" CG PHE C 25 " 0.033 2.00e-02 2.50e+03 pdb=" CD1 PHE C 25 " -0.010 2.00e-02 2.50e+03 pdb=" CD2 PHE C 25 " -0.011 2.00e-02 2.50e+03 pdb=" CE1 PHE C 25 " 0.001 2.00e-02 2.50e+03 pdb=" CE2 PHE C 25 " 0.001 2.00e-02 2.50e+03 pdb=" CZ PHE C 25 " -0.003 2.00e-02 2.50e+03 ... (remaining 1337 not shown) Histogram of nonbonded interaction distances: 2.34 - 2.85: 3453 2.85 - 3.37: 9714 3.37 - 3.88: 22898 3.88 - 4.39: 28712 4.39 - 4.90: 39523 Nonbonded interactions: 104300 Sorted by model distance: nonbonded pdb=" OG1 THR D 90 " pdb=" OE1 GLU D 93 " model vdw 2.343 3.040 nonbonded pdb=" OG1 THR H 90 " pdb=" OE1 GLU H 93 " model vdw 2.344 3.040 nonbonded pdb=" NE2 GLN A 93 " pdb=" OE2 GLU A 97 " model vdw 2.345 3.120 nonbonded pdb=" O6 DG I -21 " pdb=" N4 DC J 21 " model vdw 2.381 3.120 nonbonded pdb=" O6 DG I -10 " pdb=" N4 DC J 10 " model vdw 2.391 3.120 ... (remaining 104295 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 38 through 134) selection = chain 'E' } ncs_group { reference = (chain 'B' and ((resid 23 and (name N or name CA or name C or name O or name CB \ )) or resid 24 through 101)) selection = chain 'F' } ncs_group { reference = chain 'C' selection = chain 'G' } ncs_group { reference = chain 'D' selection = chain 'H' } 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=1.00 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: 3.480 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.000 Extract box with map and model: 0.000 Check model and map are aligned: 0.000 Set scattering table: 1.260 Process input model: 31.930 Find NCS groups from input model: 0.260 Set up NCS constraints: 0.050 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.010 Load rotamer database and sin/cos tables:2.830 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 39.820 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.7997 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.009 0.080 12861 Z= 0.524 Angle : 0.871 11.158 18629 Z= 0.501 Chirality : 0.049 0.255 2119 Planarity : 0.005 0.044 1340 Dihedral : 26.766 176.131 5557 Min Nonbonded Distance : 2.343 Molprobity Statistics. All-atom Clashscore : 2.77 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.34 % Favored : 98.66 % Rotamer: Outliers : 1.59 % Allowed : 3.34 % Favored : 95.07 % 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: -5.83 (0.19), residues: 744 helix: -4.27 (0.11), residues: 524 sheet: None (None), residues: 0 loop : -2.04 (0.34), residues: 220 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.012 0.003 HIS C 82 PHE 0.033 0.004 PHE C 25 TYR 0.045 0.004 TYR F 88 ARG 0.014 0.002 ARG G 81 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 195 residues out of total 632 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 10 poor density : 185 time to evaluate : 0.987 Fit side-chains REVERT: A 56 LYS cc_start: 0.7914 (ttpt) cc_final: 0.7632 (ttpp) REVERT: A 107 THR cc_start: 0.9122 (t) cc_final: 0.8878 (p) REVERT: A 120 MET cc_start: 0.8199 (mtm) cc_final: 0.7969 (mtp) REVERT: B 93 GLN cc_start: 0.8220 (mt0) cc_final: 0.8005 (mt0) REVERT: C 50 TYR cc_start: 0.9089 (t80) cc_final: 0.8850 (t80) REVERT: C 76 THR cc_start: 0.8500 (p) cc_final: 0.8023 (t) REVERT: E 97 GLU cc_start: 0.7922 (mt-10) cc_final: 0.7645 (mt-10) REVERT: G 75 LYS cc_start: 0.7932 (mttm) cc_final: 0.7661 (mtmt) REVERT: D 40 TYR cc_start: 0.8764 (m-80) cc_final: 0.7998 (m-80) REVERT: D 116 LYS cc_start: 0.8499 (ttpt) cc_final: 0.8288 (mtpt) REVERT: H 37 TYR cc_start: 0.8612 (m-80) cc_final: 0.8385 (m-80) REVERT: H 120 LYS cc_start: 0.7980 (tttp) cc_final: 0.7540 (ttmm) outliers start: 10 outliers final: 0 residues processed: 189 average time/residue: 0.4181 time to fit residues: 98.9882 Evaluate side-chains 134 residues out of total 632 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 0 poor density : 134 time to evaluate : 1.066 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 83 optimal weight: 0.9990 chunk 75 optimal weight: 2.9990 chunk 41 optimal weight: 0.6980 chunk 25 optimal weight: 1.9990 chunk 50 optimal weight: 0.6980 chunk 40 optimal weight: 0.7980 chunk 77 optimal weight: 3.9990 chunk 30 optimal weight: 0.9980 chunk 47 optimal weight: 0.6980 chunk 57 optimal weight: 3.9990 chunk 89 optimal weight: 0.9990 overall best weight: 0.7780 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** A 68 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** B 25 ASN ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 84 GLN C 104 GLN E 39 HIS E 108 ASN G 24 GLN G 82 HIS D 47 GLN H 95 GLN Total number of N/Q/H flips: 9 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3702 r_free = 0.3702 target = 0.147930 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 40)----------------| | r_work = 0.3221 r_free = 0.3221 target = 0.110556 restraints weight = 15829.281| |-----------------------------------------------------------------------------| r_work (start): 0.3175 rms_B_bonded: 0.86 r_work: 0.3047 rms_B_bonded: 2.34 restraints_weight: 0.5000 r_work: 0.2941 rms_B_bonded: 3.79 restraints_weight: 0.2500 r_work (final): 0.2941 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8202 moved from start: 0.2454 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.042 12861 Z= 0.214 Angle : 0.615 5.685 18629 Z= 0.370 Chirality : 0.037 0.155 2119 Planarity : 0.005 0.035 1340 Dihedral : 30.712 170.708 4059 Min Nonbonded Distance : 2.141 Molprobity Statistics. All-atom Clashscore : 6.97 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.21 % Favored : 98.79 % Rotamer: Outliers : 1.27 % Allowed : 7.15 % Favored : 91.57 % 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: -2.50 (0.26), residues: 744 helix: -1.53 (0.19), residues: 538 sheet: None (None), residues: 0 loop : -1.49 (0.37), residues: 206 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.004 0.001 HIS C 82 PHE 0.010 0.002 PHE A 67 TYR 0.023 0.002 TYR B 51 ARG 0.003 0.000 ARG F 40 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 161 residues out of total 632 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 8 poor density : 153 time to evaluate : 0.811 Fit side-chains REVERT: A 56 LYS cc_start: 0.8198 (ttpt) cc_final: 0.7963 (ttpp) REVERT: A 76 GLN cc_start: 0.7912 (mt0) cc_final: 0.7682 (mt0) REVERT: A 107 THR cc_start: 0.9238 (t) cc_final: 0.9015 (p) REVERT: A 120 MET cc_start: 0.8654 (mtm) cc_final: 0.8402 (mtp) REVERT: B 79 LYS cc_start: 0.8780 (mptt) cc_final: 0.8064 (mptt) REVERT: C 76 THR cc_start: 0.8706 (p) cc_final: 0.8430 (t) REVERT: E 90 MET cc_start: 0.7867 (mmp) cc_final: 0.7426 (mmp) REVERT: G 75 LYS cc_start: 0.7894 (mttm) cc_final: 0.7662 (mtmt) REVERT: D 36 SER cc_start: 0.9114 (p) cc_final: 0.8778 (p) REVERT: D 43 LYS cc_start: 0.8756 (mtmm) cc_final: 0.8372 (mttt) REVERT: H 120 LYS cc_start: 0.8017 (tttp) cc_final: 0.7661 (ttmm) outliers start: 8 outliers final: 4 residues processed: 159 average time/residue: 0.3826 time to fit residues: 77.3217 Evaluate side-chains 138 residues out of total 632 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 4 poor density : 134 time to evaluate : 0.884 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 46 VAL Chi-restraints excluded: chain F residue 96 THR Chi-restraints excluded: chain D residue 112 SER Chi-restraints excluded: chain H residue 47 GLN Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 99 random chunks: chunk 77 optimal weight: 3.9990 chunk 86 optimal weight: 2.9990 chunk 52 optimal weight: 0.9980 chunk 22 optimal weight: 0.7980 chunk 80 optimal weight: 4.9990 chunk 58 optimal weight: 3.9990 chunk 55 optimal weight: 4.9990 chunk 15 optimal weight: 1.9990 chunk 61 optimal weight: 3.9990 chunk 59 optimal weight: 4.9990 chunk 35 optimal weight: 1.9990 overall best weight: 1.7586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** C 82 HIS ** both conformations clash, **PLEASE CHECK MANUALLY** C 104 GLN D 109 HIS H 84 ASN H 95 GLN Total number of N/Q/H flips: 4 ------------------------------------------------------------------------------- ADP refinement ************** |-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------| | r_work = 0.3591 r_free = 0.3591 target = 0.139457 restraints weight = None | |-----------------------------------------------------------------------------| |-group b-factor refinement (macro cycle = 1; iterations = 45)----------------| | r_work = 0.3091 r_free = 0.3091 target = 0.101519 restraints weight = 16177.827| |-----------------------------------------------------------------------------| r_work (start): 0.3097 rms_B_bonded: 0.95 r_work: 0.2963 rms_B_bonded: 2.43 restraints_weight: 0.5000 r_work: 0.2851 rms_B_bonded: 3.91 restraints_weight: 0.2500 r_work (final): 0.2851 ------------------------------------------------------------------------------- Occupancy refinement ******************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8263 moved from start: 0.2780 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.006 0.053 12861 Z= 0.327 Angle : 0.650 5.872 18629 Z= 0.382 Chirality : 0.040 0.171 2119 Planarity : 0.006 0.045 1340 Dihedral : 31.036 171.024 4059 Min Nonbonded Distance : 2.131 Molprobity Statistics. All-atom Clashscore : 5.49 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.88 % Favored : 98.12 % Rotamer: Outliers : 1.43 % Allowed : 10.02 % 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.83 (0.29), residues: 744 helix: -0.18 (0.21), residues: 531 sheet: None (None), residues: 0 loop : -1.20 (0.38), residues: 213 Max deviation from planes: Type MaxDev MeanDev LineInFile HIS 0.004 0.001 HIS G 31 PHE 0.009 0.002 PHE F 100 TYR 0.020 0.002 TYR B 51 ARG 0.005 0.001 ARG B 45 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 1488 Ramachandran restraints generated. 744 Oldfield, 0 Emsley, 744 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 146 residues out of total 632 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 9 poor density : 137 time to evaluate : 0.822 Fit side-chains revert: symmetry clash REVERT: A 56 LYS cc_start: 0.8301 (ttpt) cc_final: 0.8050 (ttpp) REVERT: A 107 THR cc_start: 0.9274 (t) cc_final: 0.9031 (p) REVERT: A 120 MET cc_start: 0.8637 (mtm) cc_final: 0.8406 (mtp) REVERT: B 79 LYS cc_start: 0.8839 (mptt) cc_final: 0.7995 (mptt) REVERT: C 76 THR cc_start: 0.8654 (p) cc_final: 0.8278 (t) REVERT: G 75 LYS cc_start: 0.7963 (mttm) cc_final: 0.7737 (mtmt) REVERT: D 36 SER cc_start: 0.9140 (p) cc_final: 0.8857 (p) REVERT: D 43 LYS cc_start: 0.8708 (mtmm) cc_final: 0.8369 (mttt) outliers start: 9 outliers final: 8 residues processed: 143 average time/residue: 0.3757 time to fit residues: 68.7803 Evaluate side-chains 138 residues out of total 632 (non-[ALA,GLY,PRO]) need to be fit. rotamer outliers: 8 poor density : 130 time to evaluate : 0.907 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain A residue 46 VAL Chi-restraints excluded: chain A residue 97 GLU Chi-restraints excluded: chain B residue 30 THR Chi-restraints excluded: chain F residue 49 LEU Chi-restraints excluded: chain F residue 96 THR Chi-restraints excluded: chain G residue 101 THR Chi-restraints excluded: chain D residue 123 SER Chi-restraints excluded: chain H residue 47 GLN Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/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-5571/modules/cctbx_project/iotbx/cli_parser.py", line 946, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/phenix/phenix/programs/real_space_refine.py", line 210, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/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-5571/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/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-5571/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 767, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/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-5571/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-5571/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-5571/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-5571/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-5571/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-5571/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/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-5571/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-5571/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1525, in energies_sites flags=flags, sites_cart=sites_cart, site_labels=site_labels) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1427, in pair_proxies check_bonded_distance_cutoff(sites_frac=sites_frac) File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5571/modules/cctbx_project/cctbx/geometry_restraints/manager.py", line 1308, in check_bonded_distance_cutoff raise RuntimeError(msg) RuntimeError: Bond distance > max_reasonable_bond_distance: 70.447 > 50: distance: 140 - 142: 19.551 distance: 142 - 143: 4.365 distance: 143 - 144: 21.565 distance: 143 - 146: 22.359 distance: 144 - 145: 13.093 distance: 144 - 153: 30.152 distance: 146 - 147: 22.675 distance: 147 - 148: 18.522 distance: 148 - 149: 10.219 distance: 149 - 150: 5.568 distance: 150 - 151: 3.845 distance: 153 - 154: 11.970 distance: 154 - 155: 38.738 distance: 154 - 157: 48.470 distance: 155 - 156: 18.818 distance: 155 - 160: 28.378 distance: 157 - 158: 35.446 distance: 157 - 159: 31.536 distance: 160 - 161: 21.046 distance: 161 - 162: 26.883 distance: 161 - 164: 14.539 distance: 162 - 163: 22.814 distance: 162 - 167: 30.497 distance: 164 - 165: 13.234 distance: 164 - 166: 4.388 distance: 167 - 168: 25.358 distance: 168 - 169: 26.449 distance: 168 - 171: 30.377 distance: 169 - 170: 16.434 distance: 169 - 175: 32.382 distance: 171 - 172: 18.302 distance: 171 - 173: 10.398 distance: 172 - 174: 10.825 distance: 175 - 176: 20.958 distance: 176 - 177: 33.097 distance: 176 - 179: 10.801 distance: 177 - 178: 17.098 distance: 177 - 180: 29.409 distance: 180 - 181: 17.329 distance: 181 - 182: 29.477 distance: 181 - 184: 17.976 distance: 182 - 183: 16.432 distance: 182 - 189: 22.345 distance: 184 - 185: 9.951 distance: 185 - 186: 14.218 distance: 186 - 187: 10.451 distance: 186 - 188: 8.501 distance: 189 - 190: 22.628 distance: 190 - 191: 13.198 distance: 191 - 192: 9.056 distance: 191 - 193: 13.092 distance: 193 - 194: 9.765 distance: 194 - 195: 19.035 distance: 195 - 196: 16.323 distance: 195 - 197: 23.906 distance: 197 - 198: 25.207 distance: 198 - 199: 35.689 distance: 198 - 201: 22.654 distance: 199 - 200: 21.791 distance: 199 - 204: 23.396 distance: 201 - 202: 12.202 distance: 201 - 203: 10.400 distance: 204 - 205: 8.106 distance: 205 - 206: 29.161 distance: 205 - 208: 44.569 distance: 206 - 207: 15.910 distance: 206 - 212: 4.601 distance: 208 - 209: 8.672 distance: 209 - 210: 33.267 distance: 209 - 211: 15.777