Starting phenix.real_space_refine on Thu Jan 18 05:46:01 2024 by dcliebschner =============================================================================== Processing files: ------------------------------------------------------------------------------- Found model, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7utn_26782/01_2024/7utn_26782.pdb Found real_map, /net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7utn_26782/01_2024/7utn_26782.map Processing PHIL parameters: ------------------------------------------------------------------------------- Adding command-line PHIL: ------------------------- refinement.macro_cycles=10 scattering_table=electron resolution=2.74 write_initial_geo_file=False Final processed PHIL parameters: ------------------------------------------------------------------------------- data_manager { real_map_files = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7utn_26782/01_2024/7utn_26782.map" default_real_map = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7utn_26782/01_2024/7utn_26782.map" model { file = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7utn_26782/01_2024/7utn_26782.pdb" } default_model = "/net/marbles/raid1/dorothee/rerefine/cryoem/data_no_H/7utn_26782/01_2024/7utn_26782.pdb" } resolution = 2.74 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.001 sd= 0.026 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 228 5.49 5 S 21 5.16 5 C 4183 2.51 5 N 1470 2.21 5 O 2101 1.98 5 sf(0) = scattering factor at diffraction angle 0. Process input model Symmetric amino acids flipped Residue "D ARG 19": "NH1" <-> "NH2" Residue "D GLU 29": "OE1" <-> "OE2" Residue "D GLU 47": "OE1" <-> "OE2" Residue "D GLU 50": "OE1" <-> "OE2" Residue "D GLU 51": "OE1" <-> "OE2" Residue "D GLU 125": "OE1" <-> "OE2" Residue "D GLU 192": "OE1" <-> "OE2" Residue "D TYR 340": "CD1" <-> "CD2" "CE1" <-> "CE2" Residue "D GLU 410": "OE1" <-> "OE2" Residue "D GLU 411": "OE1" <-> "OE2" Residue "D GLU 459": "OE1" <-> "OE2" Residue "D GLU 473": "OE1" <-> "OE2" Time to flip residues: 0.01s Monomer Library directory: "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/modules/chem_data/mon_lib" Total number of atoms: 8003 Number of models: 1 Model: "" Number of chains: 4 Chain: "D" Number of atoms: 3174 Number of conformers: 1 Conformer: "" Number of residues, atoms: 397, 3174 Classifications: {'peptide': 397} Link IDs: {'PCIS': 1, 'PTRANS': 17, 'TRANS': 378} Chain breaks: 1 Chain: "C" Number of atoms: 3937 Number of conformers: 1 Conformer: "" Number of residues, atoms: 184, 3937 Classifications: {'RNA': 184} Modifications used: {'rna2p_pur': 14, 'rna2p_pyr': 12, 'rna3p_pur': 90, 'rna3p_pyr': 68} Link IDs: {'rna2p': 25, 'rna3p': 158} Chain breaks: 2 Chain: "A" Number of atoms: 579 Number of conformers: 1 Conformer: "" Number of residues, atoms: 29, 579 Classifications: {'DNA': 29} Link IDs: {'rna3p': 28} Chain: "B" Number of atoms: 313 Number of conformers: 1 Conformer: "" Number of residues, atoms: 15, 313 Classifications: {'DNA': 15} Link IDs: {'rna3p': 14} Time building chain proxies: 5.18, per 1000 atoms: 0.65 Number of scatterers: 8003 At special positions: 0 Unit cell: (90.95, 105.93, 131.61, 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 21 16.00 P 228 15.00 O 2101 8.00 N 1470 7.00 C 4183 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.50 Amino acid : False - 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: 4.65 Conformation dependent library (CDL) restraints added in 821.6 milliseconds 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. Adding C-beta torsion restraints... Number of C-beta restraints generated: 744 Finding SS restraints... Secondary structure from input PDB file: 9 helices and 6 sheets defined 30.0% alpha, 20.4% beta 76 base pairs and 137 stacking pairs defined. Time for finding SS restraints: 3.71 Creating SS restraints... Processing helix chain 'D' and resid 21 through 29 Processing helix chain 'D' and resid 88 through 109 removed outlier: 4.313A pdb=" N LYS D 92 " --> pdb=" O LYS D 88 " (cutoff:3.500A) Processing helix chain 'D' and resid 109 through 121 Processing helix chain 'D' and resid 163 through 184 removed outlier: 3.662A pdb=" N ILE D 183 " --> pdb=" O LYS D 179 " (cutoff:3.500A) Processing helix chain 'D' and resid 298 through 314 Proline residue: D 304 - end of helix Processing helix chain 'D' and resid 323 through 334 Processing helix chain 'D' and resid 340 through 348 removed outlier: 4.588A pdb=" N TYR D 344 " --> pdb=" O TYR D 340 " (cutoff:3.500A) Processing helix chain 'D' and resid 408 through 416 Processing helix chain 'D' and resid 417 through 422 Processing sheet with id=AA1, first strand: chain 'D' and resid 5 through 7 removed outlier: 6.264A pdb=" N TYR D 5 " --> pdb=" O ILE D 42 " (cutoff:3.500A) removed outlier: 7.458A pdb=" N LEU D 44 " --> pdb=" O TYR D 5 " (cutoff:3.500A) removed outlier: 7.028A pdb=" N ILE D 7 " --> pdb=" O LEU D 44 " (cutoff:3.500A) removed outlier: 6.402A pdb=" N THR D 41 " --> pdb=" O VAL D 35 " (cutoff:3.500A) Processing sheet with id=AA2, first strand: chain 'D' and resid 318 through 322 removed outlier: 6.231A pdb=" N VAL D 189 " --> pdb=" O CYS D 319 " (cutoff:3.500A) removed outlier: 7.629A pdb=" N THR D 321 " --> pdb=" O VAL D 189 " (cutoff:3.500A) removed outlier: 6.312A pdb=" N LEU D 191 " --> pdb=" O THR D 321 " (cutoff:3.500A) removed outlier: 6.246A pdb=" N LEU D 57 " --> pdb=" O VAL D 190 " (cutoff:3.500A) removed outlier: 7.845A pdb=" N GLU D 192 " --> pdb=" O LEU D 57 " (cutoff:3.500A) removed outlier: 7.460A pdb=" N ILE D 59 " --> pdb=" O GLU D 192 " (cutoff:3.500A) removed outlier: 3.607A pdb=" N PHE D 79 " --> pdb=" O VAL D 70 " (cutoff:3.500A) removed outlier: 3.611A pdb=" N GLN D 369 " --> pdb=" O GLU D 84 " (cutoff:3.500A) removed outlier: 4.864A pdb=" N MET D 366 " --> pdb=" O TYR D 494 " (cutoff:3.500A) Processing sheet with id=AA3, first strand: chain 'D' and resid 128 through 131 Processing sheet with id=AA4, first strand: chain 'D' and resid 378 through 380 removed outlier: 6.559A pdb=" N GLN D 432 " --> pdb=" O HIS D 379 " (cutoff:3.500A) Processing sheet with id=AA5, first strand: chain 'D' and resid 392 through 395 removed outlier: 6.610A pdb=" N TYR D 387 " --> pdb=" O VAL D 394 " (cutoff:3.500A) Processing sheet with id=AA6, first strand: chain 'D' and resid 476 through 478 removed outlier: 4.836A pdb=" N VAL D 466 " --> pdb=" O GLY D 460 " (cutoff:3.500A) removed outlier: 7.246A pdb=" N GLY D 460 " --> pdb=" O VAL D 466 " (cutoff:3.500A) removed outlier: 6.659A pdb=" N TYR D 468 " --> pdb=" O SER D 458 " (cutoff:3.500A) removed outlier: 4.662A pdb=" N SER D 458 " --> pdb=" O TYR D 468 " (cutoff:3.500A) removed outlier: 6.750A pdb=" N VAL D 470 " --> pdb=" O SER D 456 " (cutoff:3.500A) removed outlier: 6.694A pdb=" N ILE D 444 " --> pdb=" O LEU D 485 " (cutoff:3.500A) 118 hydrogen bonds defined for protein. 330 hydrogen bond angles defined for protein. Restraints generated for nucleic acids: 196 hydrogen bonds 364 hydrogen bond angles 0 basepair planarities 76 basepair parallelities 137 stacking parallelities Total time for adding SS restraints: 3.41 Time building geometry restraints manager: 4.83 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: 1144 1.33 - 1.45: 3372 1.45 - 1.56: 3635 1.56 - 1.68: 451 1.68 - 1.80: 33 Bond restraints: 8635 Sorted by residual: bond pdb=" O3' U C 205 " pdb=" P U C 206 " ideal model delta sigma weight residual 1.607 1.572 0.035 1.50e-02 4.44e+03 5.50e+00 bond pdb=" N ASN D 65 " pdb=" CA ASN D 65 " ideal model delta sigma weight residual 1.466 1.443 0.023 1.36e-02 5.41e+03 2.93e+00 bond pdb=" CA ARG D 151 " pdb=" CB ARG D 151 " ideal model delta sigma weight residual 1.527 1.497 0.030 1.75e-02 3.27e+03 2.93e+00 bond pdb=" SD MET D 68 " pdb=" CE MET D 68 " ideal model delta sigma weight residual 1.791 1.748 0.043 2.50e-02 1.60e+03 2.92e+00 bond pdb=" N ARG D 330 " pdb=" CA ARG D 330 " ideal model delta sigma weight residual 1.459 1.441 0.018 1.17e-02 7.31e+03 2.43e+00 ... (remaining 8630 not shown) Histogram of bond angle deviations from ideal: 96.49 - 104.00: 837 104.00 - 111.50: 4587 111.50 - 119.01: 2929 119.01 - 126.51: 3708 126.51 - 134.02: 675 Bond angle restraints: 12736 Sorted by residual: angle pdb=" C4' U C 206 " pdb=" C3' U C 206 " pdb=" O3' U C 206 " ideal model delta sigma weight residual 109.40 117.98 -8.58 1.50e+00 4.44e-01 3.27e+01 angle pdb=" C3' A C 128 " pdb=" C2' A C 128 " pdb=" C1' A C 128 " ideal model delta sigma weight residual 101.50 105.10 -3.60 1.00e+00 1.00e+00 1.30e+01 angle pdb=" C4' A C 128 " pdb=" C3' A C 128 " pdb=" O3' A C 128 " ideal model delta sigma weight residual 109.40 114.76 -5.36 1.50e+00 4.44e-01 1.27e+01 angle pdb=" C3' A C 78 " pdb=" C2' A C 78 " pdb=" C1' A C 78 " ideal model delta sigma weight residual 101.30 104.71 -3.41 1.00e+00 1.00e+00 1.16e+01 angle pdb=" C2' A C 78 " pdb=" C1' A C 78 " pdb=" N9 A C 78 " ideal model delta sigma weight residual 112.00 106.96 5.04 1.50e+00 4.44e-01 1.13e+01 ... (remaining 12731 not shown) Histogram of dihedral angle deviations from ideal: 0.00 - 35.57: 4595 35.57 - 71.15: 477 71.15 - 106.72: 52 106.72 - 142.29: 4 142.29 - 177.87: 10 Dihedral angle restraints: 5138 sinusoidal: 3983 harmonic: 1155 Sorted by residual: dihedral pdb=" O4' C C 60 " pdb=" C1' C C 60 " pdb=" N1 C C 60 " pdb=" C2 C C 60 " ideal model delta sinusoidal sigma weight residual 200.00 36.86 163.14 1 1.50e+01 4.44e-03 8.35e+01 dihedral pdb=" O4' U C 204 " pdb=" C1' U C 204 " pdb=" N1 U C 204 " pdb=" C2 U C 204 " ideal model delta sinusoidal sigma weight residual 200.00 46.30 153.70 1 1.50e+01 4.44e-03 8.09e+01 dihedral pdb=" O4' U C 180 " pdb=" C1' U C 180 " pdb=" N1 U C 180 " pdb=" C2 U C 180 " ideal model delta sinusoidal sigma weight residual -160.00 -28.49 -131.51 1 1.50e+01 4.44e-03 7.09e+01 ... (remaining 5135 not shown) Histogram of chiral volume deviations from ideal: 0.000 - 0.072: 1455 0.072 - 0.144: 101 0.144 - 0.216: 7 0.216 - 0.288: 3 0.288 - 0.361: 1 Chirality restraints: 1567 Sorted by residual: chirality pdb=" C3' A C 128 " pdb=" C4' A C 128 " pdb=" O3' A C 128 " pdb=" C2' A C 128 " both_signs ideal model delta sigma weight residual False -2.74 -2.38 -0.36 2.00e-01 2.50e+01 3.25e+00 chirality pdb=" C2' A C 78 " pdb=" C3' A C 78 " pdb=" O2' A C 78 " pdb=" C1' A C 78 " both_signs ideal model delta sigma weight residual False -2.75 -2.48 -0.28 2.00e-01 2.50e+01 1.92e+00 chirality pdb=" C1' G C 121 " pdb=" O4' G C 121 " pdb=" C2' G C 121 " pdb=" N9 G C 121 " both_signs ideal model delta sigma weight residual False 2.46 2.22 0.23 2.00e-01 2.50e+01 1.38e+00 ... (remaining 1564 not shown) Planarity restraints: 787 Sorted by residual: delta sigma weight rms_deltas residual plane pdb=" C1' A C 78 " -0.042 2.00e-02 2.50e+03 2.12e-02 1.24e+01 pdb=" N9 A C 78 " 0.055 2.00e-02 2.50e+03 pdb=" C8 A C 78 " -0.002 2.00e-02 2.50e+03 pdb=" N7 A C 78 " -0.005 2.00e-02 2.50e+03 pdb=" C5 A C 78 " -0.004 2.00e-02 2.50e+03 pdb=" C6 A C 78 " -0.005 2.00e-02 2.50e+03 pdb=" N6 A C 78 " -0.007 2.00e-02 2.50e+03 pdb=" N1 A C 78 " -0.002 2.00e-02 2.50e+03 pdb=" C2 A C 78 " 0.002 2.00e-02 2.50e+03 pdb=" N3 A C 78 " 0.004 2.00e-02 2.50e+03 pdb=" C4 A C 78 " 0.006 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C1' G C 91 " -0.034 2.00e-02 2.50e+03 1.45e-02 6.27e+00 pdb=" N9 G C 91 " 0.034 2.00e-02 2.50e+03 pdb=" C8 G C 91 " 0.004 2.00e-02 2.50e+03 pdb=" N7 G C 91 " -0.001 2.00e-02 2.50e+03 pdb=" C5 G C 91 " 0.001 2.00e-02 2.50e+03 pdb=" C6 G C 91 " -0.003 2.00e-02 2.50e+03 pdb=" O6 G C 91 " -0.010 2.00e-02 2.50e+03 pdb=" N1 G C 91 " -0.003 2.00e-02 2.50e+03 pdb=" C2 G C 91 " 0.001 2.00e-02 2.50e+03 pdb=" N2 G C 91 " 0.001 2.00e-02 2.50e+03 pdb=" N3 G C 91 " 0.006 2.00e-02 2.50e+03 pdb=" C4 G C 91 " 0.006 2.00e-02 2.50e+03 delta sigma weight rms_deltas residual plane pdb=" C THR D 163 " -0.037 5.00e-02 4.00e+02 5.61e-02 5.04e+00 pdb=" N PRO D 164 " 0.097 5.00e-02 4.00e+02 pdb=" CA PRO D 164 " -0.029 5.00e-02 4.00e+02 pdb=" CD PRO D 164 " -0.031 5.00e-02 4.00e+02 ... (remaining 784 not shown) Histogram of nonbonded interaction distances: 2.08 - 2.64: 226 2.64 - 3.21: 6504 3.21 - 3.77: 15156 3.77 - 4.34: 20706 4.34 - 4.90: 28300 Nonbonded interactions: 70892 Sorted by model distance: nonbonded pdb=" O2' A C 78 " pdb=" O4' G C 79 " model vdw 2.080 2.440 nonbonded pdb=" O4 U C 30 " pdb=" N6 A C 46 " model vdw 2.215 2.520 nonbonded pdb=" NH1 ARG D 115 " pdb=" OP2 C C 58 " model vdw 2.223 2.520 nonbonded pdb=" OH TYR D 46 " pdb=" OP1 U C 206 " model vdw 2.258 2.440 nonbonded pdb=" OD1 ASN D 487 " pdb=" N ASN D 488 " model vdw 2.321 2.520 ... (remaining 70887 not shown) NOTE: a complete listing of the restraints can be obtained by requesting output of .geo file. Find NCS groups from input model ========== WARNING! ============ No NCS relation were found !!! ================================ Found NCS groups: found none. 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.010 Normalize map: mean=0, sd=1: 1.450 Set stop_for_unknowns flag: 0.000 Assert model is a single copy model: 0.000 Assert all atoms have isotropic ADPs: 0.010 Construct map_model_manager: 0.000 Extract box with map and model: 5.080 Check model and map are aligned: 0.370 Set scattering table: 0.070 Process input model: 29.720 Find NCS groups from input model: 0.160 Set up NCS constraints: 0.020 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:12.190 Set ADP refinement strategy: 0.000 Make a string to write initial .geo file:0.000 Internal consistency checks: 0.000 Total: 49.080 ------------------------------------------------------------------------------- Set refinement monitor ********************** ------------------------------------------------------------------------------- Setup refinement engine *********************** ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8643 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.007 0.043 8635 Z= 0.407 Angle : 0.679 8.584 12736 Z= 0.394 Chirality : 0.041 0.361 1567 Planarity : 0.004 0.056 787 Dihedral : 22.970 177.868 4394 Min Nonbonded Distance : 2.080 Molprobity Statistics. All-atom Clashscore : 8.58 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.78 % Favored : 98.22 % Rotamer: Outliers : 0.00 % Allowed : 0.00 % Favored : 100.00 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.56 % 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.01 (0.40), residues: 393 helix: 1.09 (0.47), residues: 117 sheet: -0.44 (0.62), residues: 63 loop : -0.40 (0.39), residues: 213 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.002 TRP D 161 HIS 0.003 0.001 HIS D 104 PHE 0.010 0.002 PHE D 453 TYR 0.010 0.001 TYR D 5 ARG 0.007 0.001 ARG D 486 *********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 30 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 30 time to evaluate : 0.539 Fit side-chains REVERT: D 102 MET cc_start: 0.8634 (mmt) cc_final: 0.8320 (mmt) outliers start: 0 outliers final: 0 residues processed: 30 average time/residue: 1.6898 time to fit residues: 53.3704 Evaluate side-chains 22 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 0 poor density : 22 time to evaluate : 0.536 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Rotamers are restrained with sigma=5.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 60 random chunks: chunk 50 optimal weight: 2.9990 chunk 45 optimal weight: 4.9990 chunk 25 optimal weight: 4.9990 chunk 15 optimal weight: 3.9990 chunk 30 optimal weight: 0.8980 chunk 24 optimal weight: 1.9990 chunk 47 optimal weight: 2.9990 chunk 18 optimal weight: 0.3980 chunk 28 optimal weight: 2.9990 chunk 35 optimal weight: 1.9990 chunk 54 optimal weight: 3.9990 overall best weight: 1.6586 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** D 87 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 492 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8633 moved from start: 0.0614 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.005 0.057 8635 Z= 0.295 Angle : 0.605 6.506 12736 Z= 0.330 Chirality : 0.038 0.325 1567 Planarity : 0.005 0.056 787 Dihedral : 23.852 177.019 3581 Min Nonbonded Distance : 2.141 Molprobity Statistics. All-atom Clashscore : 8.00 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.78 % Favored : 98.22 % Rotamer: Outliers : 0.58 % Allowed : 4.65 % Favored : 94.77 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.56 % 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.39 (0.40), residues: 393 helix: 1.67 (0.47), residues: 118 sheet: -0.19 (0.60), residues: 74 loop : -0.37 (0.40), residues: 201 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.009 0.001 TRP D 161 HIS 0.003 0.001 HIS D 104 PHE 0.008 0.001 PHE D 453 TYR 0.013 0.001 TYR D 327 ARG 0.004 0.001 ARG D 363 *********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 27 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 25 time to evaluate : 0.528 Fit side-chains REVERT: D 102 MET cc_start: 0.8663 (mmt) cc_final: 0.8366 (mmt) outliers start: 2 outliers final: 2 residues processed: 25 average time/residue: 1.5905 time to fit residues: 42.1259 Evaluate side-chains 23 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 21 time to evaluate : 0.568 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 68 MET Chi-restraints excluded: chain D residue 361 LYS Rotamers are restrained with sigma=4.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 60 random chunks: chunk 30 optimal weight: 3.9990 chunk 16 optimal weight: 3.9990 chunk 45 optimal weight: 5.9990 chunk 37 optimal weight: 3.9990 chunk 15 optimal weight: 3.9990 chunk 54 optimal weight: 4.9990 chunk 58 optimal weight: 2.9990 chunk 48 optimal weight: 3.9990 chunk 18 optimal weight: 0.9990 chunk 43 optimal weight: 3.9990 chunk 53 optimal weight: 2.9990 overall best weight: 2.9990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** D 87 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8679 moved from start: 0.0770 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.008 0.078 8635 Z= 0.496 Angle : 0.687 7.229 12736 Z= 0.371 Chirality : 0.043 0.391 1567 Planarity : 0.005 0.061 787 Dihedral : 23.904 177.184 3581 Min Nonbonded Distance : 2.139 Molprobity Statistics. All-atom Clashscore : 8.50 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.53 % Favored : 98.47 % Rotamer: Outliers : 0.87 % Allowed : 5.81 % Favored : 93.31 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.56 % 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.36 (0.40), residues: 393 helix: 1.57 (0.46), residues: 119 sheet: -0.22 (0.60), residues: 75 loop : -0.33 (0.41), residues: 199 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.012 0.002 TRP D 161 HIS 0.004 0.001 HIS D 428 PHE 0.009 0.002 PHE D 453 TYR 0.023 0.002 TYR D 327 ARG 0.004 0.001 ARG D 363 *********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 26 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 23 time to evaluate : 0.550 Fit side-chains REVERT: D 102 MET cc_start: 0.8717 (mmt) cc_final: 0.8386 (mmt) outliers start: 3 outliers final: 3 residues processed: 24 average time/residue: 2.0628 time to fit residues: 51.9876 Evaluate side-chains 24 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 21 time to evaluate : 0.643 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 36 GLU Chi-restraints excluded: chain D residue 68 MET Chi-restraints excluded: chain D residue 361 LYS Rotamers are restrained with sigma=4.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 60 random chunks: chunk 41 optimal weight: 3.9990 chunk 28 optimal weight: 0.8980 chunk 6 optimal weight: 0.8980 chunk 26 optimal weight: 0.7980 chunk 36 optimal weight: 0.9980 chunk 54 optimal weight: 3.9990 chunk 57 optimal weight: 3.9990 chunk 51 optimal weight: 3.9990 chunk 15 optimal weight: 2.9990 chunk 48 optimal weight: 2.9990 chunk 32 optimal weight: 1.9990 overall best weight: 1.1182 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** D 87 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 492 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8609 moved from start: 0.1062 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.043 8635 Z= 0.207 Angle : 0.562 7.300 12736 Z= 0.305 Chirality : 0.035 0.267 1567 Planarity : 0.004 0.058 787 Dihedral : 23.706 176.928 3581 Min Nonbonded Distance : 2.164 Molprobity Statistics. All-atom Clashscore : 8.36 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.27 % Favored : 98.73 % Rotamer: Outliers : 0.87 % Allowed : 6.98 % Favored : 92.15 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.56 % 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.70 (0.41), residues: 393 helix: 1.96 (0.46), residues: 119 sheet: 0.07 (0.63), residues: 70 loop : -0.24 (0.40), residues: 204 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP D 161 HIS 0.002 0.001 HIS D 104 PHE 0.007 0.001 PHE D 453 TYR 0.009 0.001 TYR D 484 ARG 0.005 0.000 ARG D 363 *********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 28 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 25 time to evaluate : 0.503 Fit side-chains outliers start: 3 outliers final: 2 residues processed: 26 average time/residue: 1.4840 time to fit residues: 41.1010 Evaluate side-chains 25 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 23 time to evaluate : 0.712 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 68 MET Chi-restraints excluded: chain D residue 361 LYS Rotamers are restrained with sigma=3.50 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 60 random chunks: chunk 0 optimal weight: 4.9990 chunk 43 optimal weight: 3.9990 chunk 23 optimal weight: 5.9990 chunk 49 optimal weight: 3.9990 chunk 40 optimal weight: 3.9990 chunk 29 optimal weight: 2.9990 chunk 52 optimal weight: 5.9990 chunk 14 optimal weight: 4.9990 chunk 19 optimal weight: 3.9990 chunk 11 optimal weight: 0.9990 chunk 34 optimal weight: 0.9990 overall best weight: 2.5990 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** D 87 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 492 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8664 moved from start: 0.0951 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.007 0.071 8635 Z= 0.431 Angle : 0.643 7.421 12736 Z= 0.348 Chirality : 0.041 0.362 1567 Planarity : 0.005 0.059 787 Dihedral : 23.813 176.945 3581 Min Nonbonded Distance : 2.132 Molprobity Statistics. All-atom Clashscore : 8.72 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.78 % Favored : 98.22 % Rotamer: Outliers : 1.16 % Allowed : 8.72 % Favored : 90.12 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.56 % 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.60 (0.41), residues: 393 helix: 1.84 (0.46), residues: 119 sheet: -0.15 (0.61), residues: 75 loop : -0.20 (0.41), residues: 199 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.011 0.002 TRP D 161 HIS 0.003 0.001 HIS D 104 PHE 0.009 0.001 PHE D 453 TYR 0.010 0.001 TYR D 5 ARG 0.005 0.001 ARG D 363 *********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 28 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 4 poor density : 24 time to evaluate : 0.577 Fit side-chains outliers start: 4 outliers final: 3 residues processed: 26 average time/residue: 1.7018 time to fit residues: 46.6995 Evaluate side-chains 24 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 3 poor density : 21 time to evaluate : 0.546 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 36 GLU Chi-restraints excluded: chain D residue 68 MET Chi-restraints excluded: chain D residue 361 LYS Rotamers are restrained with sigma=3.00 ------------------------------------------------------------------------------- XYZ refinement ************** Weight determination summary: number of chunks: 60 random chunks: chunk 14 optimal weight: 0.7980 chunk 58 optimal weight: 2.9990 chunk 48 optimal weight: 3.9990 chunk 26 optimal weight: 1.9990 chunk 4 optimal weight: 0.9980 chunk 19 optimal weight: 4.9990 chunk 30 optimal weight: 0.9990 chunk 55 optimal weight: 10.0000 chunk 6 optimal weight: 1.9990 chunk 33 optimal weight: 0.9980 chunk 42 optimal weight: 3.9990 overall best weight: 1.1584 ------------------------------------------------------------------------------- NQH flips ********* Analyzing N/Q/H residues for possible flip corrections... Flipped N/Q/H residues before XYZ refinement: ** D 87 ASN ** both conformations clash, **PLEASE CHECK MANUALLY** ** D 492 GLN ** both conformations clash, **PLEASE CHECK MANUALLY** Total number of N/Q/H flips: 0 ------------------------------------------------------------------------------- Overall statistics ****************** model-to-map fit, CC_mask: 0.8613 moved from start: 0.1165 Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2 Deviations from Ideal Values - rmsd, rmsZ for bonds and angles. Bond : 0.004 0.044 8635 Z= 0.220 Angle : 0.567 7.459 12736 Z= 0.308 Chirality : 0.035 0.274 1567 Planarity : 0.004 0.057 787 Dihedral : 23.698 176.928 3581 Min Nonbonded Distance : 2.155 Molprobity Statistics. All-atom Clashscore : 8.65 Ramachandran Plot: Outliers : 0.00 % Allowed : 1.27 % Favored : 98.73 % Rotamer: Outliers : 1.45 % Allowed : 8.14 % Favored : 90.41 % Cbeta Deviations : 0.00 % Peptide Plane: Cis-proline : 5.56 % 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.79 (0.41), residues: 393 helix: 2.06 (0.47), residues: 119 sheet: 0.15 (0.65), residues: 70 loop : -0.22 (0.40), residues: 204 Max deviation from planes: Type MaxDev MeanDev LineInFile TRP 0.008 0.001 TRP D 161 HIS 0.002 0.001 HIS D 104 PHE 0.007 0.001 PHE D 453 TYR 0.009 0.001 TYR D 484 ARG 0.005 0.000 ARG D 363 *********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************ ------------------------------------------------------------------------------- Update Rama plot phi/psi targets (oldfield only) ************************************************ 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Re-set Ramachandran plot restraints *********************************** favored: oldfield allowed: oldfield outlier: oldfield 786 Ramachandran restraints generated. 393 Oldfield, 0 Emsley, 393 emsley8k and 0 Phi/Psi/2. ------------------------------------------------------------------------------- Optimize residue side-chains **************************** Evaluate side-chains 30 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 5 poor density : 25 time to evaluate : 0.577 Fit side-chains outliers start: 5 outliers final: 2 residues processed: 27 average time/residue: 1.5308 time to fit residues: 43.9157 Evaluate side-chains 22 residues out of total 344 non-(ALA, GLY, PRO) need fitting. rotamer outliers: 2 poor density : 20 time to evaluate : 0.573 ------------------------------------------------------------------------------- Set rotamer restraints ********************** Chi-restraints excluded: chain D residue 68 MET Chi-restraints excluded: chain D residue 361 LYS Rotamers are restrained with sigma=2.50 ------------------------------------------------------------------------------- XYZ refinement ************** Traceback (most recent call last): File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/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-1.21-5211/modules/cctbx_project/iotbx/cli_parser.py", line 931, in run_program task.run() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/modules/phenix/phenix/programs/real_space_refine.py", line 191, in run log = self.logger) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__ log = log) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/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-1.21-5211/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller func() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 700, in refine_xyz self.minimization_no_ncs() File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 764, in minimization_no_ncs log = self.log) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/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-1.21-5211/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-1.21-5211/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-1.21-5211/modules/cctbx_project/mmtbx/refinement/real_space/individual_sites.py", line 235, in __init__ weight = weight) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/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-1.21-5211/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-1.21-5211/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run line_search) File "/net/cci-filer2/raid1/xp/phenix/phenix-1.21-5211/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-1.21-5211/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-1.21-5211/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-1.21-5211/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-1.21-5211/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: 63.0284 > 50: distance: 31 - 35: 10.624 distance: 35 - 36: 7.958 distance: 36 - 39: 8.326 distance: 37 - 38: 3.718 distance: 37 - 41: 3.659 distance: 39 - 40: 11.425 distance: 42 - 43: 7.199 distance: 43 - 44: 6.310 distance: 43 - 52: 5.066 distance: 45 - 46: 10.136 distance: 46 - 47: 3.942 distance: 48 - 49: 3.315 distance: 49 - 51: 3.719 distance: 52 - 53: 8.281 distance: 53 - 54: 19.660 distance: 53 - 56: 3.815 distance: 54 - 55: 10.945 distance: 54 - 60: 18.500 distance: 56 - 57: 17.400 distance: 56 - 58: 15.757 distance: 57 - 59: 8.581 distance: 60 - 61: 4.822 distance: 61 - 62: 22.407 distance: 61 - 64: 23.317 distance: 62 - 63: 22.893 distance: 62 - 68: 5.772 distance: 64 - 65: 18.427 distance: 65 - 66: 27.703 distance: 66 - 67: 5.002 distance: 68 - 69: 17.344 distance: 68 - 74: 15.520 distance: 69 - 70: 12.904 distance: 69 - 72: 15.751 distance: 70 - 71: 16.283 distance: 70 - 75: 20.453 distance: 72 - 73: 22.722 distance: 73 - 74: 13.677 distance: 75 - 76: 29.663 distance: 76 - 77: 13.136 distance: 77 - 78: 28.938 distance: 77 - 79: 6.840 distance: 79 - 80: 8.654 distance: 80 - 81: 4.855 distance: 80 - 83: 9.296 distance: 81 - 82: 10.537 distance: 81 - 85: 9.484 distance: 83 - 84: 23.095 distance: 85 - 86: 12.875 distance: 86 - 87: 6.682 distance: 86 - 89: 8.914 distance: 87 - 88: 6.675 distance: 87 - 93: 22.837 distance: 88 - 186: 8.969 distance: 89 - 90: 12.867 distance: 89 - 91: 17.295 distance: 90 - 92: 9.616 distance: 93 - 94: 19.317 distance: 93 - 132: 24.375 distance: 94 - 95: 21.063 distance: 94 - 97: 9.184 distance: 95 - 96: 18.186 distance: 95 - 101: 16.563 distance: 96 - 129: 18.763 distance: 97 - 98: 6.527 distance: 98 - 99: 20.353 distance: 98 - 100: 22.096 distance: 101 - 102: 10.087 distance: 101 - 176: 27.693 distance: 102 - 103: 8.820 distance: 102 - 105: 12.349 distance: 103 - 104: 15.815 distance: 103 - 108: 24.532 distance: 104 - 173: 26.322 distance: 105 - 106: 18.154 distance: 105 - 107: 20.357