Starting phenix.real_space_refine
on Wed Jun 11 13:06:18 2025 by dcliebschner
===============================================================================

Processing files:
-------------------------------------------------------------------------------

  Found model, /net/cci-nas-00/data/ceres_data/8vdo_43152/06_2025/8vdo_43152.cif
  Found real_map, /net/cci-nas-00/data/ceres_data/8vdo_43152/06_2025/8vdo_43152.map

Processing PHIL parameters:
-------------------------------------------------------------------------------

  Adding command-line PHIL:
  -------------------------
    refinement.macro_cycles=10
    scattering_table=electron
    resolution=2.7
    write_initial_geo_file=False

Final processed PHIL parameters:
-------------------------------------------------------------------------------
  data_manager {
    real_map_files = "/net/cci-nas-00/data/ceres_data/8vdo_43152/06_2025/8vdo_43152.map"
    default_real_map = "/net/cci-nas-00/data/ceres_data/8vdo_43152/06_2025/8vdo_43152.map"
    model {
      file = "/net/cci-nas-00/data/ceres_data/8vdo_43152/06_2025/8vdo_43152.cif"
    }
    default_model = "/net/cci-nas-00/data/ceres_data/8vdo_43152/06_2025/8vdo_43152.cif"
  }
  resolution = 2.7
  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.040
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: 4
    Type Number    sf(0)   Gaussians
     S      62      5.16       5
     C    7794      2.51       5
     N    2267      2.21       5
     O    2568      1.98       5
    sf(0) = scattering factor at diffraction angle 0.
Process input model

  Symmetric amino acids flipped. Time to flip 61 residue(s): 0.06s

  Monomer Library directory:
    "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/chem_data/mon_lib"
  Total number of atoms: 12691
  Number of models: 1
  Model: ""
    Number of chains: 1
    Chain: "A"
      Number of atoms: 12691
      Number of conformers: 2
      Conformer: "A"
        Number of residues, atoms: 1740, 12658
          Classifications: {'peptide': 1740}
          Incomplete info: {'truncation_to_alanine': 2}
          Link IDs: {'PTRANS': 65, 'TRANS': 1674}
          Chain breaks: 1
          Unresolved non-hydrogen bonds: 6
          Unresolved non-hydrogen angles: 8
          Unresolved non-hydrogen dihedrals: 4
          Unresolved non-hydrogen chiralities: 2
      Conformer: "B"
        Number of residues, atoms: 1740, 12658
          Classifications: {'peptide': 1740}
          Incomplete info: {'truncation_to_alanine': 2}
          Link IDs: {'PTRANS': 65, 'TRANS': 1674}
          Chain breaks: 1
          Unresolved non-hydrogen bonds: 6
          Unresolved non-hydrogen angles: 8
          Unresolved non-hydrogen dihedrals: 4
          Unresolved non-hydrogen chiralities: 2
          bond proxies already assigned to first conformer: 12754
  Residues with excluded nonbonded symmetry interactions: 3
    residue:
      pdb=" N  AARG A 589 " occ=0.50
      ... (20 atoms not shown)
      pdb=" NH2BARG A 589 " occ=0.50
    residue:
      pdb=" N  AARG A 891 " occ=0.50
      ... (20 atoms not shown)
      pdb=" NH2BARG A 891 " occ=0.50
    residue:
      pdb=" N  AARG A1136 " occ=0.50
      ... (20 atoms not shown)
      pdb=" NH2BARG A1136 " occ=0.50
  Time building chain proxies: 14.90, per 1000 atoms: 1.17
  Number of scatterers: 12691
  At special positions: 0
  Unit cell: (111.744, 119.808, 148.608, 90, 90, 90)
  Space group: P 1 (No. 1)
  Number of sites at special positions: 0
  Number of scattering types: 4
    Type Number    sf(0)
     S      62     16.00
     O    2568      8.00
     N    2267      7.00
     C    7794      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: 3.96
  Conformation dependent library (CDL) restraints added in 3.6 seconds
  

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
  Adding C-beta torsion restraints...
  Number of C-beta restraints generated:  3266

  Finding SS restraints...
    Secondary structure from input PDB file:
      57 helices and 2 sheets defined
      84.9% alpha, 1.8% beta
      0 base pairs and 0 stacking pairs defined.
    Time for finding SS restraints: 1.34
  Creating SS restraints...
    Processing helix  chain 'A' and resid 384 through 400
    Processing helix  chain 'A' and resid 489 through 515
      removed outlier: 3.835A  pdb=" N   ASP A 514 " --> pdb=" O   GLN A 510 " (cutoff:3.500A)
      removed outlier: 3.786A  pdb=" N   ASP A 515 " --> pdb=" O   ALA A 511 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 525 through 561
    Processing helix  chain 'A' and resid 569 through 600
    Processing helix  chain 'A' and resid 601 through 603
      No H-bonds generated for 'chain 'A' and resid 601 through 603'
    Processing helix  chain 'A' and resid 605 through 626
    Processing helix  chain 'A' and resid 633 through 655
      removed outlier: 3.928A  pdb=" N   LEU A 637 " --> pdb=" O   PRO A 633 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 661 through 693
    Processing helix  chain 'A' and resid 695 through 724
    Processing helix  chain 'A' and resid 729 through 756
    Processing helix  chain 'A' and resid 759 through 789
      removed outlier: 3.747A  pdb=" N   ALA A 789 " --> pdb=" O   VAL A 785 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 799 through 815
    Processing helix  chain 'A' and resid 818 through 847
    Processing helix  chain 'A' and resid 849 through 880
    Processing helix  chain 'A' and resid 883 through 905
      removed outlier: 3.576A  pdb=" N   GLN A 887 " --> pdb=" O   SER A 883 " (cutoff:3.500A)
      removed outlier: 3.512A  pdb=" N   ALA A 903 " --> pdb=" O   MET A 899 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 907 through 941
    Processing helix  chain 'A' and resid 948 through 976
      Proline residue:  A 964  - end of helix
    Processing helix  chain 'A' and resid 979 through 1008
      removed outlier: 4.451A  pdb=" N   GLN A 995 " --> pdb=" O   GLN A 991 " (cutoff:3.500A)
      Proline residue:  A 996  - end of helix
    Processing helix  chain 'A' and resid 1009 through 1011
      No H-bonds generated for 'chain 'A' and resid 1009 through 1011'
    Processing helix  chain 'A' and resid 1013 through 1046
    Processing helix  chain 'A' and resid 1048 through 1073
      removed outlier: 3.771A  pdb=" N   LEU A1054 " --> pdb=" O   MET A1050 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1083 through 1112
    Processing helix  chain 'A' and resid 1113 through 1142
    Processing helix  chain 'A' and resid 1144 through 1175
      removed outlier: 3.553A  pdb=" N   GLN A1148 " --> pdb=" O   ASP A1144 " (cutoff:3.500A)
      removed outlier: 3.538A  pdb=" N   HIS A1175 " --> pdb=" O   LYS A1171 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1178 through 1202
    Processing helix  chain 'A' and resid 1203 through 1205
      No H-bonds generated for 'chain 'A' and resid 1203 through 1205'
    Processing helix  chain 'A' and resid 1206 through 1225
    Processing helix  chain 'A' and resid 1234 through 1260
    Processing helix  chain 'A' and resid 1263 through 1291
    Processing helix  chain 'A' and resid 1294 through 1325
      removed outlier: 3.889A  pdb=" N   THR A1324 " --> pdb=" O   LYS A1320 " (cutoff:3.500A)
      removed outlier: 3.633A  pdb=" N   ASP A1325 " --> pdb=" O   ALA A1321 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1328 through 1353
    Processing helix  chain 'A' and resid 1359 through 1374
      removed outlier: 3.582A  pdb=" N   CYS A1363 " --> pdb=" O   GLY A1359 " (cutoff:3.500A)
      removed outlier: 3.930A  pdb=" N   THR A1372 " --> pdb=" O   ARG A1368 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1375 through 1378
      removed outlier: 3.629A  pdb=" N   GLU A1378 " --> pdb=" O   GLU A1375 " (cutoff:3.500A)
      No H-bonds generated for 'chain 'A' and resid 1375 through 1378'
    Processing helix  chain 'A' and resid 1388 through 1417
    Processing helix  chain 'A' and resid 1418 through 1450
    Processing helix  chain 'A' and resid 1463 through 1483
    Processing helix  chain 'A' and resid 1487 through 1516
    Processing helix  chain 'A' and resid 1518 through 1547
    Processing helix  chain 'A' and resid 1551 through 1577
      removed outlier: 3.884A  pdb=" N   ALA A1563 " --> pdb=" O   ARG A1559 " (cutoff:3.500A)
      Proline residue:  A1564  - end of helix
    Processing helix  chain 'A' and resid 1589 through 1623
      Proline residue:  A1598  - end of helix
    Processing helix  chain 'A' and resid 1626 through 1655
      removed outlier: 3.720A  pdb=" N   SER A1631 " --> pdb=" O   PRO A1627 " (cutoff:3.500A)
      removed outlier: 3.584A  pdb=" N   ALA A1655 " --> pdb=" O   MET A1651 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1657 through 1684
    Processing helix  chain 'A' and resid 1694 through 1722
      removed outlier: 3.555A  pdb=" N   LEU A1698 " --> pdb=" O   SER A1694 " (cutoff:3.500A)
      removed outlier: 4.596A  pdb=" N   GLU A1714 " --> pdb=" O   SER A1710 " (cutoff:3.500A)
      Proline residue:  A1715  - end of helix
    Processing helix  chain 'A' and resid 1723 through 1736
      removed outlier: 3.507A  pdb=" N   LEU A1727 " --> pdb=" O   GLU A1723 " (cutoff:3.500A)
      removed outlier: 4.356A  pdb=" N   GLY A1728 " --> pdb=" O   ALA A1724 " (cutoff:3.500A)
      removed outlier: 4.343A  pdb=" N   HIS A1729 " --> pdb=" O   SER A1725 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1737 through 1752
    Processing helix  chain 'A' and resid 1754 through 1783
    Processing helix  chain 'A' and resid 1789 through 1821
      removed outlier: 3.721A  pdb=" N   ALA A1821 " --> pdb=" O   ALA A1817 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1826 through 1838
    Processing helix  chain 'A' and resid 1849 through 1878
      removed outlier: 3.924A  pdb=" N   ASN A1876 " --> pdb=" O   VAL A1872 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1881 through 1907
      removed outlier: 3.985A  pdb=" N   LYS A1901 " --> pdb=" O   ALA A1897 " (cutoff:3.500A)
      Proline residue:  A1902  - end of helix
    Processing helix  chain 'A' and resid 1909 through 1940
    Processing helix  chain 'A' and resid 1943 through 1972
      removed outlier: 3.520A  pdb=" N   LYS A1947 " --> pdb=" O   ASP A1943 " (cutoff:3.500A)
      removed outlier: 3.747A  pdb=" N   GLY A1971 " --> pdb=" O   ALA A1967 " (cutoff:3.500A)
      removed outlier: 4.127A  pdb=" N   ASN A1972 " --> pdb=" O   LEU A1968 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 1972 through 2001
      removed outlier: 4.143A  pdb=" N   ALA A2001 " --> pdb=" O   MET A1997 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 2014 through 2037
      removed outlier: 4.268A  pdb=" N   GLY A2018 " --> pdb=" O   ASP A2014 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 2040 through 2070
      removed outlier: 3.562A  pdb=" N   GLY A2070 " --> pdb=" O   ALA A2066 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 2073 through 2102
      removed outlier: 3.547A  pdb=" N   GLN A2077 " --> pdb=" O   ASP A2073 " (cutoff:3.500A)
      removed outlier: 3.790A  pdb=" N   LEU A2080 " --> pdb=" O   THR A2076 " (cutoff:3.500A)
    Processing helix  chain 'A' and resid 2108 through 2135
      removed outlier: 4.123A  pdb=" N   TRP A2112 " --> pdb=" O   ASP A2108 " (cutoff:3.500A)
    Processing sheet with id=AA1, first strand: chain 'A' and resid 311 through 318
      removed outlier: 3.670A  pdb=" N   SER A 311 " --> pdb=" O   ILE A 332 " (cutoff:3.500A)
      removed outlier: 7.146A  pdb=" N   VAL A 337 " --> pdb=" O   GLU A 350 " (cutoff:3.500A)
      removed outlier: 5.083A  pdb=" N   GLU A 350 " --> pdb=" O   VAL A 337 " (cutoff:3.500A)
      removed outlier: 6.802A  pdb=" N   ARG A 339 " --> pdb=" O   ILE A 348 " (cutoff:3.500A)
    Processing sheet with id=AA2, first strand: chain 'A' and resid 358 through 361

    1226 hydrogen bonds defined for protein.
    3672 hydrogen bond angles defined for protein.
    Restraints generated for nucleic acids:
      0 hydrogen bonds
      0 hydrogen bond angles
      0 basepair planarities
      0 basepair parallelities
      0 stacking parallelities
  Total time for adding SS restraints: 5.82

  Time building geometry restraints manager: 4.63 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: 2651
        1.33 -     1.45: 2458
        1.45 -     1.57: 7620
        1.57 -     1.69: 0
        1.69 -     1.81: 97
  Bond restraints: 12826
  Sorted by residual:
  bond pdb=" C   GLU A1463 "
       pdb=" N   PRO A1464 "
    ideal  model  delta    sigma   weight residual
    1.336  1.371 -0.035 1.20e-02 6.94e+03 8.75e+00
  bond pdb=" C   PRO A1464 "
       pdb=" O   PRO A1464 "
    ideal  model  delta    sigma   weight residual
    1.237  1.206  0.031 1.20e-02 6.94e+03 6.68e+00
  bond pdb=" N   THR A1465 "
       pdb=" CA  THR A1465 "
    ideal  model  delta    sigma   weight residual
    1.459  1.488 -0.029 1.21e-02 6.83e+03 5.57e+00
  bond pdb=" N   ARG A1469 "
       pdb=" CA  ARG A1469 "
    ideal  model  delta    sigma   weight residual
    1.459  1.487 -0.028 1.21e-02 6.83e+03 5.46e+00
  bond pdb=" N   GLN A1466 "
       pdb=" CA  GLN A1466 "
    ideal  model  delta    sigma   weight residual
    1.459  1.487 -0.028 1.21e-02 6.83e+03 5.22e+00
  ... (remaining 12821 not shown)

  Histogram of bond angle deviations from ideal:
        0.00 -     2.02: 16968
        2.02 -     4.04: 373
        4.04 -     6.06: 49
        6.06 -     8.08: 8
        8.08 -    10.09: 2
  Bond angle restraints: 17400
  Sorted by residual:
  angle pdb=" C   THR A 772 "
        pdb=" CA  THR A 772 "
        pdb=" CB  THR A 772 "
      ideal   model   delta    sigma   weight residual
     110.79  117.82   -7.03 1.66e+00 3.63e-01 1.79e+01
  angle pdb=" CA  LEU A1020 "
        pdb=" C   LEU A1020 "
        pdb=" O   LEU A1020 "
      ideal   model   delta    sigma   weight residual
     120.55  116.77    3.78 1.06e+00 8.90e-01 1.27e+01
  angle pdb=" C   GLN A 783 "
        pdb=" N   HIS A 784 "
        pdb=" CA  HIS A 784 "
      ideal   model   delta    sigma   weight residual
     122.06  115.49    6.57 1.86e+00 2.89e-01 1.25e+01
  angle pdb=" CA  SER A1021 "
        pdb=" C   SER A1021 "
        pdb=" O   SER A1021 "
      ideal   model   delta    sigma   weight residual
     120.55  116.89    3.66 1.06e+00 8.90e-01 1.19e+01
  angle pdb=" CB  MET A1759 "
        pdb=" CG  MET A1759 "
        pdb=" SD  MET A1759 "
      ideal   model   delta    sigma   weight residual
     112.70  102.61   10.09 3.00e+00 1.11e-01 1.13e+01
  ... (remaining 17395 not shown)

  Histogram of dihedral angle deviations from ideal:
        0.00 -    17.96: 7287
       17.96 -    35.92: 532
       35.92 -    53.88: 102
       53.88 -    71.85: 59
       71.85 -    89.81: 21
  Dihedral angle restraints: 8001
    sinusoidal: 2977
      harmonic: 5024
  Sorted by residual:
  dihedral pdb=" CA  GLU A 599 "
           pdb=" C   GLU A 599 "
           pdb=" N   ASP A 600 "
           pdb=" CA  ASP A 600 "
      ideal   model   delta  harmonic     sigma   weight residual
     180.00  159.61   20.39     0      5.00e+00 4.00e-02 1.66e+01
  dihedral pdb=" CA  VAL A 380 "
           pdb=" C   VAL A 380 "
           pdb=" N   GLN A 381 "
           pdb=" CA  GLN A 381 "
      ideal   model   delta  harmonic     sigma   weight residual
     180.00 -161.09  -18.91     0      5.00e+00 4.00e-02 1.43e+01
  dihedral pdb=" C   THR A 772 "
           pdb=" N   THR A 772 "
           pdb=" CA  THR A 772 "
           pdb=" CB  THR A 772 "
      ideal   model   delta  harmonic     sigma   weight residual
    -122.00 -131.32    9.32     0      2.50e+00 1.60e-01 1.39e+01
  ... (remaining 7998 not shown)

  Histogram of chiral volume deviations from ideal:
       0.000 -    0.043: 1591
       0.043 -    0.086: 406
       0.086 -    0.129: 84
       0.129 -    0.172: 15
       0.172 -    0.215: 6
  Chirality restraints: 2102
  Sorted by residual:
  chirality pdb=" CG  LEU A 609 "
            pdb=" CB  LEU A 609 "
            pdb=" CD1 LEU A 609 "
            pdb=" CD2 LEU A 609 "
    both_signs  ideal   model   delta    sigma   weight residual
      False     -2.59   -2.38   -0.21 2.00e-01 2.50e+01 1.15e+00
  chirality pdb=" CB  VAL A 380 "
            pdb=" CA  VAL A 380 "
            pdb=" CG1 VAL A 380 "
            pdb=" CG2 VAL A 380 "
    both_signs  ideal   model   delta    sigma   weight residual
      False     -2.63   -2.42   -0.21 2.00e-01 2.50e+01 1.08e+00
  chirality pdb=" CA  ARG A1035 "
            pdb=" N   ARG A1035 "
            pdb=" C   ARG A1035 "
            pdb=" CB  ARG A1035 "
    both_signs  ideal   model   delta    sigma   weight residual
      False      2.51    2.71   -0.20 2.00e-01 2.50e+01 1.02e+00
  ... (remaining 2099 not shown)

  Planarity restraints: 2313
  Sorted by residual:
                                 delta    sigma   weight rms_deltas residual
  plane pdb=" CA  THR A 772 "    0.013 2.00e-02 2.50e+03   2.67e-02 7.14e+00
        pdb=" C   THR A 772 "   -0.046 2.00e-02 2.50e+03
        pdb=" O   THR A 772 "    0.017 2.00e-02 2.50e+03
        pdb=" N   ALA A 773 "    0.016 2.00e-02 2.50e+03
                                 delta    sigma   weight rms_deltas residual
  plane pdb=" C   ASP A 564 "   -0.039 5.00e-02 4.00e+02   5.85e-02 5.48e+00
        pdb=" N   PRO A 565 "    0.101 5.00e-02 4.00e+02
        pdb=" CA  PRO A 565 "   -0.030 5.00e-02 4.00e+02
        pdb=" CD  PRO A 565 "   -0.033 5.00e-02 4.00e+02
                                 delta    sigma   weight rms_deltas residual
  plane pdb=" CA  MET A1018 "   -0.011 2.00e-02 2.50e+03   2.16e-02 4.65e+00
        pdb=" C   MET A1018 "    0.037 2.00e-02 2.50e+03
        pdb=" O   MET A1018 "   -0.014 2.00e-02 2.50e+03
        pdb=" N   GLN A1019 "   -0.013 2.00e-02 2.50e+03
  ... (remaining 2310 not shown)

  Histogram of nonbonded interaction distances:
        2.20 -     2.74: 1296
        2.74 -     3.28: 13686
        3.28 -     3.82: 23089
        3.82 -     4.36: 25619
        4.36 -     4.90: 44352
  Nonbonded interactions: 108042
  Sorted by model distance:
  nonbonded pdb=" ND2 ASN A1250 "
            pdb=" OE1 GLU A1371 "
     model   vdw
     2.198 3.120
  nonbonded pdb=" OE2 GLU A1168 "
            pdb=" NH1 ARG A1184 "
     model   vdw
     2.200 3.120
  nonbonded pdb=" OD1 ASN A1211 "
            pdb=" NH2 ARG A1214 "
     model   vdw
     2.203 3.120
  nonbonded pdb=" O   VAL A1640 "
            pdb=" OG  SER A1643 "
     model   vdw
     2.205 3.040
  nonbonded pdb=" O   LYS A2031 "
            pdb=" NE2 GLN A2035 "
     model   vdw
     2.231 3.120
  ... (remaining 108037 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=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.230
  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.550
  Check model and map are aligned:         0.090
  Set scattering table:                    0.140
  Process input model:                     45.660
  Find NCS groups from input model:        0.210
  Set up NCS constraints:                  0.080
  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:15.210
  Set ADP refinement strategy:             0.000
  Make a string to write initial .geo file:0.000
  Internal consistency checks:             0.000
  Total:   64.190
-------------------------------------------------------------------------------
Set refinement monitor
**********************
-------------------------------------------------------------------------------
Setup refinement engine
***********************
-------------------------------------------------------------------------------
Overall statistics
******************
model-to-map fit, CC_mask: 0.7035
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.075  12826  Z= 0.228
  Angle     :  0.718  10.095  17400  Z= 0.425
  Chirality :  0.041   0.215   2102
  Planarity :  0.005   0.059   2313
  Dihedral  : 15.243  89.807   4735
  Min Nonbonded Distance : 2.198

Molprobity Statistics.
  All-atom Clashscore : 9.48
  Ramachandran Plot:
    Outliers :  0.00 %
    Allowed  :  0.92 %
    Favored  : 99.08 %
  Rotamer:
    Outliers :  0.45 %
    Allowed  :  0.30 %
    Favored  : 99.25 %
  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:  1.94 (0.19), residues: 1745
  helix:  1.81 (0.13), residues: 1428
  sheet: -0.22 (0.70), residues: 53
  loop : -1.23 (0.34), residues: 264

Max deviation from planes:
   Type  MaxDev  MeanDev LineInFile
 TRP   0.011   0.002   TRP A 351 
 HIS   0.004   0.001   HIS A 784 
 PHE   0.022   0.002   PHE A1390 
 TYR   0.023   0.002   TYR A1893 
 ARG   0.008   0.001   ARG A 889 

Details of bonding type rmsd
  hydrogen bonds       : bond        0.12064 ( 1226)
  hydrogen bonds       : angle       4.78087 ( 3672)
  covalent geometry    : bond        0.00481 (12826)
  covalent geometry    : angle       0.71807 (17400)

*********************** REFINEMENT MACRO_CYCLE 1 OF 10 ************************

-------------------------------------------------------------------------------
Update Rama plot phi/psi targets (oldfield only)
************************************************

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Re-set Ramachandran plot restraints
***********************************
  favored: oldfield
  allowed: oldfield
  outlier: oldfield

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Optimize residue side-chains
****************************
Evaluate side-chains
  155 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 2
    poor density    : 153
  time to evaluate  : 1.446 
Fit side-chains
REVERT: A  772 THR cc_start: 0.7184 (OUTLIER) cc_final: 0.6740 (t)
REVERT: A 1184 ARG cc_start: 0.6547 (mtm110) cc_final: 0.5996 (mtp85)
REVERT: A 1304 ASN cc_start: 0.7159 (m-40) cc_final: 0.6952 (m-40)
REVERT: A 1531 GLU cc_start: 0.7420 (tp30) cc_final: 0.6990 (tp30)
  outliers start: 2
  outliers final: 0
  residues processed: 154
  average time/residue: 1.2787
  time to fit residues: 215.3046
Evaluate side-chains
  113 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 1
    poor density    : 112
  time to evaluate  : 1.491 
-------------------------------------------------------------------------------
Set rotamer restraints
**********************
Chi-restraints excluded: chain A residue  772 THR
Rotamers are restrained with sigma=5.00
-------------------------------------------------------------------------------
XYZ refinement
**************
Weight determination summary:
   number of chunks: 174
   random chunks:
   chunk 146 optimal weight:    0.5980
   chunk 131 optimal weight:    0.7980
   chunk 73 optimal weight:    3.9990
   chunk 45 optimal weight:    9.9990
   chunk 88 optimal weight:    4.9990
   chunk 70 optimal weight:    1.9990
   chunk 136 optimal weight:    0.5980
   chunk 52 optimal weight:    1.9990
   chunk 82 optimal weight:    5.9990
   chunk 101 optimal weight:    0.6980
   chunk 158 optimal weight:    1.9990
   overall best weight:    0.9382

-------------------------------------------------------------------------------
NQH flips
*********
Analyzing N/Q/H residues for possible flip corrections...

Flipped N/Q/H residues before XYZ refinement:
** A 349 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
** A 374 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
** A 776 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
A 784 HIS
A1686 GLN
A1837 GLN
A1854 GLN

Total number of N/Q/H flips: 4

-------------------------------------------------------------------------------
ADP refinement
**************
|-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------|
| r_work = 0.4241 r_free = 0.4241 target = 0.200881 restraints weight = None  |
|-----------------------------------------------------------------------------|
|-group b-factor refinement (macro cycle = 1; iterations = 46)----------------|
| r_work = 0.3635 r_free = 0.3635 target = 0.139770 restraints weight = 33087.214|
|-----------------------------------------------------------------------------|
r_work (start): 0.3603 rms_B_bonded: 3.25
r_work: 0.3434 rms_B_bonded: 3.16 restraints_weight: 0.5000
r_work (final): 0.3434
-------------------------------------------------------------------------------
Occupancy refinement
********************
r_start: 0.3478
|-occupancy refinement: start-------------------------------------------------|
| r_work = 0.3478 r_free = 0.3478    target_work(ls_wunit_k1) = 0.129         |
| occupancies: max = 1.00  min = 0.50   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
|-occupancy refinement: end---------------------------------------------------|
| r_work = 0.3478 r_free = 0.3478    target_work(ls_wunit_k1) = 0.128         |
| occupancies: max = 1.00  min = 0.29   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
r_final: 0.3478
-------------------------------------------------------------------------------
Overall statistics
******************
model-to-map fit, CC_mask: 0.7058
moved from start:          0.0875

Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2
Deviations from Ideal Values - rmsd, rmsZ for bonds and angles.
  Bond      :  0.004   0.048  12826  Z= 0.183
  Angle     :  0.607   7.924  17400  Z= 0.321
  Chirality :  0.041   0.211   2102
  Planarity :  0.005   0.054   2313
  Dihedral  :  4.358  29.814   1824
  Min Nonbonded Distance : 2.522

Molprobity Statistics.
  All-atom Clashscore : 7.95
  Ramachandran Plot:
    Outliers :  0.00 %
    Allowed  :  0.92 %
    Favored  : 99.08 %
  Rotamer:
    Outliers :  1.96 %
    Allowed  :  7.61 %
    Favored  : 90.44 %
  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.86 (0.19), residues: 1745
  helix:  2.45 (0.13), residues: 1440
  sheet:  0.09 (0.72), residues: 53
  loop : -1.16 (0.37), residues: 252

Max deviation from planes:
   Type  MaxDev  MeanDev LineInFile
 TRP   0.008   0.001   TRP A 351 
 HIS   0.015   0.002   HIS A 784 
 PHE   0.019   0.002   PHE A1390 
 TYR   0.022   0.002   TYR A1893 
 ARG   0.006   0.001   ARG A1184 

Details of bonding type rmsd
  hydrogen bonds       : bond        0.05439 ( 1226)
  hydrogen bonds       : angle       3.78341 ( 3672)
  covalent geometry    : bond        0.00415 (12826)
  covalent geometry    : angle       0.60698 (17400)

*********************** REFINEMENT MACRO_CYCLE 2 OF 10 ************************

-------------------------------------------------------------------------------
Update Rama plot phi/psi targets (oldfield only)
************************************************

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Re-set Ramachandran plot restraints
***********************************
  favored: oldfield
  allowed: oldfield
  outlier: oldfield

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Optimize residue side-chains
****************************
Evaluate side-chains
  138 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 22
    poor density    : 116
  time to evaluate  : 1.426 
Fit side-chains
REVERT: A  772 THR cc_start: 0.7460 (OUTLIER) cc_final: 0.7248 (t)
REVERT: A 1184 ARG cc_start: 0.6615 (mtm110) cc_final: 0.6325 (mtp-110)
REVERT: A 1304 ASN cc_start: 0.7071 (m-40) cc_final: 0.6860 (m-40)
REVERT: A 1397 MET cc_start: 0.6980 (OUTLIER) cc_final: 0.6716 (mtm)
REVERT: A 1455 GLN cc_start: 0.7614 (OUTLIER) cc_final: 0.7294 (mt0)
REVERT: A 1531 GLU cc_start: 0.7543 (tp30) cc_final: 0.7006 (tp30)
REVERT: A 1910 GLU cc_start: 0.7286 (mm-30) cc_final: 0.6794 (mt-10)
  outliers start: 22
  outliers final: 10
  residues processed: 126
  average time/residue: 1.3694
  time to fit residues: 188.3891
Evaluate side-chains
  124 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 13
    poor density    : 111
  time to evaluate  : 1.402 
-------------------------------------------------------------------------------
Set rotamer restraints
**********************
Chi-restraints excluded: chain A residue  314 LEU
Chi-restraints excluded: chain A residue  380 VAL
Chi-restraints excluded: chain A residue  384 GLU
Chi-restraints excluded: chain A residue  703 ILE
Chi-restraints excluded: chain A residue  772 THR
Chi-restraints excluded: chain A residue 1333 SER
Chi-restraints excluded: chain A residue 1345 SER
Chi-restraints excluded: chain A residue 1397 MET
Chi-restraints excluded: chain A residue 1449 VAL
Chi-restraints excluded: chain A residue 1455 GLN
Chi-restraints excluded: chain A residue 1520 THR
Chi-restraints excluded: chain A residue 1905 VAL
Chi-restraints excluded: chain A residue 2111 VAL
Rotamers are restrained with sigma=4.50
-------------------------------------------------------------------------------
XYZ refinement
**************
Weight determination summary:
   number of chunks: 174
   random chunks:
   chunk 159 optimal weight:    6.9990
   chunk 76 optimal weight:    2.9990
   chunk 59 optimal weight:    0.9980
   chunk 60 optimal weight:    0.9980
   chunk 55 optimal weight:    0.9990
   chunk 155 optimal weight:    0.6980
   chunk 13 optimal weight:    2.9990
   chunk 97 optimal weight:    0.4980
   chunk 44 optimal weight:    5.9990
   chunk 17 optimal weight:    0.8980
   chunk 92 optimal weight:    1.9990
   overall best weight:    0.8180

-------------------------------------------------------------------------------
NQH flips
*********
Analyzing N/Q/H residues for possible flip corrections...

Flipped N/Q/H residues before XYZ refinement:
** A 349 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
** A 374 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
A 784 HIS
A1686 GLN

Total number of N/Q/H flips: 2

-------------------------------------------------------------------------------
ADP refinement
**************
|-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------|
| r_work = 0.4247 r_free = 0.4247 target = 0.201548 restraints weight = None  |
|-----------------------------------------------------------------------------|
|-group b-factor refinement (macro cycle = 1; iterations = 46)----------------|
| r_work = 0.3632 r_free = 0.3632 target = 0.139397 restraints weight = 35514.426|
|-----------------------------------------------------------------------------|
r_work (start): 0.3600 rms_B_bonded: 3.21
r_work: 0.3435 rms_B_bonded: 3.18 restraints_weight: 0.5000
r_work (final): 0.3435
-------------------------------------------------------------------------------
Occupancy refinement
********************
r_start: 0.3472
|-occupancy refinement: start-------------------------------------------------|
| r_work = 0.3472 r_free = 0.3472    target_work(ls_wunit_k1) = 0.128         |
| occupancies: max = 1.00  min = 0.29   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
|-occupancy refinement: end---------------------------------------------------|
| r_work = 0.3472 r_free = 0.3472    target_work(ls_wunit_k1) = 0.128         |
| occupancies: max = 1.00  min = 0.26   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
r_final: 0.3472
-------------------------------------------------------------------------------
Overall statistics
******************
model-to-map fit, CC_mask: 0.7044
moved from start:          0.1107

Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2
Deviations from Ideal Values - rmsd, rmsZ for bonds and angles.
  Bond      :  0.004   0.043  12826  Z= 0.162
  Angle     :  0.565   8.759  17400  Z= 0.295
  Chirality :  0.039   0.205   2102
  Planarity :  0.004   0.054   2313
  Dihedral  :  4.162  26.888   1824
  Min Nonbonded Distance : 2.518

Molprobity Statistics.
  All-atom Clashscore : 7.76
  Ramachandran Plot:
    Outliers :  0.00 %
    Allowed  :  0.98 %
    Favored  : 99.02 %
  Rotamer:
    Outliers :  2.18 %
    Allowed  :  8.73 %
    Favored  : 89.08 %
  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:  3.39 (0.19), residues: 1745
  helix:  2.84 (0.12), residues: 1445
  sheet:  0.01 (0.71), residues: 53
  loop : -1.18 (0.38), residues: 247

Max deviation from planes:
   Type  MaxDev  MeanDev LineInFile
 TRP   0.007   0.001   TRP A 351 
 HIS   0.010   0.001   HIS A 784 
 PHE   0.017   0.002   PHE A1390 
 TYR   0.022   0.002   TYR A1893 
 ARG   0.005   0.000   ARG A 827 

Details of bonding type rmsd
  hydrogen bonds       : bond        0.04926 ( 1226)
  hydrogen bonds       : angle       3.56449 ( 3672)
  covalent geometry    : bond        0.00362 (12826)
  covalent geometry    : angle       0.56531 (17400)

*********************** REFINEMENT MACRO_CYCLE 3 OF 10 ************************

-------------------------------------------------------------------------------
Update Rama plot phi/psi targets (oldfield only)
************************************************

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Re-set Ramachandran plot restraints
***********************************
  favored: oldfield
  allowed: oldfield
  outlier: oldfield

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Optimize residue side-chains
****************************
Evaluate side-chains
  144 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 25
    poor density    : 119
  time to evaluate  : 1.314 
Fit side-chains
REVERT: A  784 HIS cc_start: 0.7172 (OUTLIER) cc_final: 0.5763 (p90)
REVERT: A 1153 ASP cc_start: 0.7082 (m-30) cc_final: 0.6754 (m-30)
REVERT: A 1184 ARG cc_start: 0.6540 (mtm110) cc_final: 0.6271 (mtp-110)
REVERT: A 1304 ASN cc_start: 0.7049 (m-40) cc_final: 0.6832 (m-40)
REVERT: A 1397 MET cc_start: 0.6982 (OUTLIER) cc_final: 0.6688 (mtm)
REVERT: A 1431 LYS cc_start: 0.8096 (tttt) cc_final: 0.6863 (tmtt)
REVERT: A 1455 GLN cc_start: 0.7661 (OUTLIER) cc_final: 0.7355 (mt0)
REVERT: A 1531 GLU cc_start: 0.7478 (tp30) cc_final: 0.6944 (tp30)
REVERT: A 1644 ILE cc_start: 0.7121 (mt) cc_final: 0.6632 (mp)
  outliers start: 25
  outliers final: 12
  residues processed: 131
  average time/residue: 1.2699
  time to fit residues: 182.2898
Evaluate side-chains
  125 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 15
    poor density    : 110
  time to evaluate  : 1.433 
-------------------------------------------------------------------------------
Set rotamer restraints
**********************
Chi-restraints excluded: chain A residue  314 LEU
Chi-restraints excluded: chain A residue  380 VAL
Chi-restraints excluded: chain A residue  384 GLU
Chi-restraints excluded: chain A residue  703 ILE
Chi-restraints excluded: chain A residue  784 HIS
Chi-restraints excluded: chain A residue 1064 LEU
Chi-restraints excluded: chain A residue 1200 VAL
Chi-restraints excluded: chain A residue 1333 SER
Chi-restraints excluded: chain A residue 1345 SER
Chi-restraints excluded: chain A residue 1397 MET
Chi-restraints excluded: chain A residue 1449 VAL
Chi-restraints excluded: chain A residue 1455 GLN
Chi-restraints excluded: chain A residue 1520 THR
Chi-restraints excluded: chain A residue 1573 SER
Chi-restraints excluded: chain A residue 1905 VAL
Rotamers are restrained with sigma=4.00
-------------------------------------------------------------------------------
XYZ refinement
**************
Weight determination summary:
   number of chunks: 174
   random chunks:
   chunk 86 optimal weight:    2.9990
   chunk 153 optimal weight:    2.9990
   chunk 21 optimal weight:    1.9990
   chunk 0 optimal weight:    2.9990
   chunk 162 optimal weight:   10.0000
   chunk 135 optimal weight:    3.9990
   chunk 8 optimal weight:    9.9990
   chunk 20 optimal weight:    1.9990
   chunk 94 optimal weight:    0.7980
   chunk 95 optimal weight:    4.9990
   chunk 52 optimal weight:    0.9990
   overall best weight:    1.7588

-------------------------------------------------------------------------------
NQH flips
*********
Analyzing N/Q/H residues for possible flip corrections...

Flipped N/Q/H residues before XYZ refinement:
** A 349 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
** A 374 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
A 665 GLN
** A 784 HIS ** both conformations clash, **PLEASE CHECK MANUALLY**
** A1058 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
A1686 GLN
A2035 GLN

Total number of N/Q/H flips: 3

-------------------------------------------------------------------------------
ADP refinement
**************
|-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------|
| r_work = 0.4220 r_free = 0.4220 target = 0.198470 restraints weight = None  |
|-----------------------------------------------------------------------------|
|-group b-factor refinement (macro cycle = 1; iterations = 48)----------------|
| r_work = 0.3599 r_free = 0.3599 target = 0.137104 restraints weight = 26374.672|
|-----------------------------------------------------------------------------|
r_work (start): 0.3595 rms_B_bonded: 2.70
r_work: 0.3440 rms_B_bonded: 2.99 restraints_weight: 0.5000
r_work: 0.3313 rms_B_bonded: 4.80 restraints_weight: 0.2500
r_work (final): 0.3313
-------------------------------------------------------------------------------
Occupancy refinement
********************
r_start: 0.3331
|-occupancy refinement: start-------------------------------------------------|
| r_work = 0.3331 r_free = 0.3331    target_work(ls_wunit_k1) = 0.116         |
| occupancies: max = 1.00  min = 0.26   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
|-occupancy refinement: end---------------------------------------------------|
| r_work = 0.3331 r_free = 0.3331    target_work(ls_wunit_k1) = 0.116         |
| occupancies: max = 1.00  min = 0.30   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
r_final: 0.3331
-------------------------------------------------------------------------------
Overall statistics
******************
model-to-map fit, CC_mask: 0.7192
moved from start:          0.1168

Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2
Deviations from Ideal Values - rmsd, rmsZ for bonds and angles.
  Bond      :  0.006   0.068  12826  Z= 0.252
  Angle     :  0.658   8.792  17400  Z= 0.345
  Chirality :  0.045   0.248   2102
  Planarity :  0.005   0.060   2313
  Dihedral  :  4.363  22.634   1823
  Min Nonbonded Distance : 2.505

Molprobity Statistics.
  All-atom Clashscore : 7.76
  Ramachandran Plot:
    Outliers :  0.00 %
    Allowed  :  1.38 %
    Favored  : 98.62 %
  Rotamer:
    Outliers :  2.64 %
    Allowed  :  9.19 %
    Favored  : 88.18 %
  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.80 (0.19), residues: 1745
  helix:  2.40 (0.12), residues: 1444
  sheet: -0.10 (0.71), residues: 53
  loop : -1.16 (0.38), residues: 248

Max deviation from planes:
   Type  MaxDev  MeanDev LineInFile
 TRP   0.005   0.002   TRP A 359 
 HIS   0.024   0.002   HIS A 784 
 PHE   0.023   0.002   PHE A1390 
 TYR   0.025   0.003   TYR A1893 
 ARG   0.008   0.001   ARG A 898 

Details of bonding type rmsd
  hydrogen bonds       : bond        0.05765 ( 1226)
  hydrogen bonds       : angle       3.74408 ( 3672)
  covalent geometry    : bond        0.00603 (12826)
  covalent geometry    : angle       0.65846 (17400)

*********************** REFINEMENT MACRO_CYCLE 4 OF 10 ************************

-------------------------------------------------------------------------------
Update Rama plot phi/psi targets (oldfield only)
************************************************

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Re-set Ramachandran plot restraints
***********************************
  favored: oldfield
  allowed: oldfield
  outlier: oldfield

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Optimize residue side-chains
****************************
Evaluate side-chains
  145 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 31
    poor density    : 114
  time to evaluate  : 1.409 
Fit side-chains
   revert: symmetry clash
REVERT: A  670 GLN cc_start: 0.7902 (pp30) cc_final: 0.7598 (pp30)
REVERT: A  776 GLN cc_start: 0.7194 (mp10) cc_final: 0.6190 (mm-40)
REVERT: A 1184 ARG cc_start: 0.6654 (mtm110) cc_final: 0.6341 (mtp-110)
REVERT: A 1304 ASN cc_start: 0.7071 (m-40) cc_final: 0.6854 (m-40)
REVERT: A 1397 MET cc_start: 0.7120 (OUTLIER) cc_final: 0.6847 (mtm)
REVERT: A 1431 LYS cc_start: 0.8175 (tttt) cc_final: 0.6953 (tmtt)
REVERT: A 1455 GLN cc_start: 0.7934 (OUTLIER) cc_final: 0.7657 (mt0)
REVERT: A 1531 GLU cc_start: 0.7703 (tp30) cc_final: 0.7064 (tp30)
REVERT: A 1644 ILE cc_start: 0.7085 (mt) cc_final: 0.6537 (mp)
  outliers start: 31
  outliers final: 18
  residues processed: 133
  average time/residue: 1.3458
  time to fit residues: 195.8477
Evaluate side-chains
  131 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 20
    poor density    : 111
  time to evaluate  : 1.469 
-------------------------------------------------------------------------------
Set rotamer restraints
**********************
Chi-restraints excluded: chain A residue  314 LEU
Chi-restraints excluded: chain A residue  380 VAL
Chi-restraints excluded: chain A residue  384 GLU
Chi-restraints excluded: chain A residue  577 VAL
Chi-restraints excluded: chain A residue  703 ILE
Chi-restraints excluded: chain A residue  784 HIS
Chi-restraints excluded: chain A residue  790 THR
Chi-restraints excluded: chain A residue 1064 LEU
Chi-restraints excluded: chain A residue 1147 VAL
Chi-restraints excluded: chain A residue 1200 VAL
Chi-restraints excluded: chain A residue 1316 LEU
Chi-restraints excluded: chain A residue 1345 SER
Chi-restraints excluded: chain A residue 1384 ILE
Chi-restraints excluded: chain A residue 1397 MET
Chi-restraints excluded: chain A residue 1449 VAL
Chi-restraints excluded: chain A residue 1455 GLN
Chi-restraints excluded: chain A residue 1520 THR
Chi-restraints excluded: chain A residue 1573 SER
Chi-restraints excluded: chain A residue 1650 SER
Chi-restraints excluded: chain A residue 1905 VAL
Rotamers are restrained with sigma=3.50
-------------------------------------------------------------------------------
XYZ refinement
**************
Weight determination summary:
   number of chunks: 174
   random chunks:
   chunk 61 optimal weight:    0.4980
   chunk 37 optimal weight:    4.9990
   chunk 45 optimal weight:    9.9990
   chunk 15 optimal weight:    0.5980
   chunk 41 optimal weight:    0.5980
   chunk 36 optimal weight:    4.9990
   chunk 94 optimal weight:    1.9990
   chunk 72 optimal weight:    0.7980
   chunk 81 optimal weight:    0.9980
   chunk 46 optimal weight:    0.8980
   chunk 160 optimal weight:    8.9990
   overall best weight:    0.6780

-------------------------------------------------------------------------------
NQH flips
*********
Analyzing N/Q/H residues for possible flip corrections...

Flipped N/Q/H residues before XYZ refinement:
** A 349 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
** A 374 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
** A 784 HIS ** both conformations clash, **PLEASE CHECK MANUALLY**
A1686 GLN

Total number of N/Q/H flips: 1

-------------------------------------------------------------------------------
ADP refinement
**************
|-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------|
| r_work = 0.4248 r_free = 0.4248 target = 0.201500 restraints weight = None  |
|-----------------------------------------------------------------------------|
|-group b-factor refinement (macro cycle = 1; iterations = 44)----------------|
| r_work = 0.3624 r_free = 0.3624 target = 0.138701 restraints weight = 35457.876|
|-----------------------------------------------------------------------------|
r_work (start): 0.3643 rms_B_bonded: 3.15
r_work: 0.3486 rms_B_bonded: 3.19 restraints_weight: 0.5000
r_work (final): 0.3486
-------------------------------------------------------------------------------
Occupancy refinement
********************
r_start: 0.3484
|-occupancy refinement: start-------------------------------------------------|
| r_work = 0.3484 r_free = 0.3484    target_work(ls_wunit_k1) = 0.129         |
| occupancies: max = 1.00  min = 0.30   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
|-occupancy refinement: end---------------------------------------------------|
| r_work = 0.3484 r_free = 0.3484    target_work(ls_wunit_k1) = 0.129         |
| occupancies: max = 1.00  min = 0.26   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
r_final: 0.3484
-------------------------------------------------------------------------------
Overall statistics
******************
model-to-map fit, CC_mask: 0.6991
moved from start:          0.1277

Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2
Deviations from Ideal Values - rmsd, rmsZ for bonds and angles.
  Bond      :  0.003   0.042  12826  Z= 0.152
  Angle     :  0.553   8.518  17400  Z= 0.287
  Chirality :  0.038   0.181   2102
  Planarity :  0.004   0.056   2313
  Dihedral  :  4.053  21.086   1823
  Min Nonbonded Distance : 2.525

Molprobity Statistics.
  All-atom Clashscore : 7.64
  Ramachandran Plot:
    Outliers :  0.00 %
    Allowed  :  1.09 %
    Favored  : 98.91 %
  Rotamer:
    Outliers :  2.03 %
    Allowed  : 10.32 %
    Favored  : 87.65 %
  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:  3.58 (0.19), residues: 1745
  helix:  2.98 (0.12), residues: 1442
  sheet: -0.00 (0.72), residues: 53
  loop : -1.13 (0.38), residues: 250

Max deviation from planes:
   Type  MaxDev  MeanDev LineInFile
 TRP   0.006   0.001   TRP A 351 
 HIS   0.010   0.001   HIS A 784 
 PHE   0.016   0.001   PHE A1390 
 TYR   0.021   0.002   TYR A1893 
 ARG   0.004   0.000   ARG A 827 

Details of bonding type rmsd
  hydrogen bonds       : bond        0.04725 ( 1226)
  hydrogen bonds       : angle       3.47844 ( 3672)
  covalent geometry    : bond        0.00335 (12826)
  covalent geometry    : angle       0.55281 (17400)

*********************** REFINEMENT MACRO_CYCLE 5 OF 10 ************************

-------------------------------------------------------------------------------
Update Rama plot phi/psi targets (oldfield only)
************************************************

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Re-set Ramachandran plot restraints
***********************************
  favored: oldfield
  allowed: oldfield
  outlier: oldfield

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Optimize residue side-chains
****************************
Evaluate side-chains
  138 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 23
    poor density    : 115
  time to evaluate  : 1.368 
Fit side-chains
   revert: symmetry clash
REVERT: A  670 GLN cc_start: 0.7598 (pp30) cc_final: 0.7337 (pp30)
REVERT: A 1304 ASN cc_start: 0.6868 (m-40) cc_final: 0.6651 (m-40)
REVERT: A 1397 MET cc_start: 0.7016 (OUTLIER) cc_final: 0.6720 (mtm)
REVERT: A 1431 LYS cc_start: 0.8069 (tttt) cc_final: 0.6819 (tmtt)
REVERT: A 1455 GLN cc_start: 0.7656 (OUTLIER) cc_final: 0.7359 (mt0)
REVERT: A 1531 GLU cc_start: 0.7461 (tp30) cc_final: 0.6881 (tp30)
REVERT: A 1644 ILE cc_start: 0.7051 (mt) cc_final: 0.6525 (mp)
  outliers start: 23
  outliers final: 14
  residues processed: 128
  average time/residue: 1.4110
  time to fit residues: 197.0611
Evaluate side-chains
  125 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 16
    poor density    : 109
  time to evaluate  : 1.312 
-------------------------------------------------------------------------------
Set rotamer restraints
**********************
Chi-restraints excluded: chain A residue  337 VAL
Chi-restraints excluded: chain A residue  380 VAL
Chi-restraints excluded: chain A residue  384 GLU
Chi-restraints excluded: chain A residue  703 ILE
Chi-restraints excluded: chain A residue  784 HIS
Chi-restraints excluded: chain A residue  790 THR
Chi-restraints excluded: chain A residue 1345 SER
Chi-restraints excluded: chain A residue 1384 ILE
Chi-restraints excluded: chain A residue 1397 MET
Chi-restraints excluded: chain A residue 1449 VAL
Chi-restraints excluded: chain A residue 1455 GLN
Chi-restraints excluded: chain A residue 1520 THR
Chi-restraints excluded: chain A residue 1571 ASN
Chi-restraints excluded: chain A residue 1573 SER
Chi-restraints excluded: chain A residue 1650 SER
Chi-restraints excluded: chain A residue 1905 VAL
Rotamers are restrained with sigma=3.00
-------------------------------------------------------------------------------
XYZ refinement
**************
Weight determination summary:
   number of chunks: 174
   random chunks:
   chunk 166 optimal weight:    0.9980
   chunk 115 optimal weight:   10.0000
   chunk 15 optimal weight:    0.5980
   chunk 4 optimal weight:    0.4980
   chunk 157 optimal weight:    2.9990
   chunk 44 optimal weight:    5.9990
   chunk 132 optimal weight:    1.9990
   chunk 130 optimal weight:    0.8980
   chunk 139 optimal weight:    0.6980
   chunk 119 optimal weight:    1.9990
   chunk 124 optimal weight:    5.9990
   overall best weight:    0.7380

-------------------------------------------------------------------------------
NQH flips
*********
Analyzing N/Q/H residues for possible flip corrections...

Flipped N/Q/H residues before XYZ refinement:
** A 349 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
** A 374 GLN ** both conformations clash, **PLEASE CHECK MANUALLY**
** A 784 HIS ** both conformations clash, **PLEASE CHECK MANUALLY**

Total number of N/Q/H flips: 0

-------------------------------------------------------------------------------
ADP refinement
**************
|-group b-factor refinement (macro cycle = 0; iterations = 0)-----------------|
| r_work = 0.4249 r_free = 0.4249 target = 0.201646 restraints weight = None  |
|-----------------------------------------------------------------------------|
|-group b-factor refinement (macro cycle = 1; iterations = 48)----------------|
| r_work = 0.3624 r_free = 0.3624 target = 0.138587 restraints weight = 26418.251|
|-----------------------------------------------------------------------------|
r_work (start): 0.3640 rms_B_bonded: 2.57
r_work: 0.3515 rms_B_bonded: 2.88 restraints_weight: 0.5000
r_work: 0.3389 rms_B_bonded: 4.61 restraints_weight: 0.2500
r_work (final): 0.3389
-------------------------------------------------------------------------------
Occupancy refinement
********************
r_start: 0.3377
|-occupancy refinement: start-------------------------------------------------|
| r_work = 0.3377 r_free = 0.3377    target_work(ls_wunit_k1) = 0.120         |
| occupancies: max = 1.00  min = 0.26   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
|-occupancy refinement: end---------------------------------------------------|
| r_work = 0.3378 r_free = 0.3378    target_work(ls_wunit_k1) = 0.120         |
| occupancies: max = 1.00  min = 0.24   number of occupancies < 0.1:        0 |
|-----------------------------------------------------------------------------|
r_final: 0.3378
-------------------------------------------------------------------------------
Overall statistics
******************
model-to-map fit, CC_mask: 0.7076
moved from start:          0.1377

Geometry Restraints Library: GeoStd + Monomer Library + CDL v1.2
Deviations from Ideal Values - rmsd, rmsZ for bonds and angles.
  Bond      :  0.003   0.041  12826  Z= 0.153
  Angle     :  0.544   7.629  17400  Z= 0.283
  Chirality :  0.038   0.174   2102
  Planarity :  0.004   0.056   2313
  Dihedral  :  3.984  20.644   1823
  Min Nonbonded Distance : 2.527

Molprobity Statistics.
  All-atom Clashscore : 7.33
  Ramachandran Plot:
    Outliers :  0.00 %
    Allowed  :  0.98 %
    Favored  : 99.02 %
  Rotamer:
    Outliers :  2.03 %
    Allowed  : 10.39 %
    Favored  : 87.58 %
  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:  3.75 (0.19), residues: 1745
  helix:  3.10 (0.12), residues: 1444
  sheet: -0.01 (0.74), residues: 53
  loop : -1.05 (0.39), residues: 248

Max deviation from planes:
   Type  MaxDev  MeanDev LineInFile
 TRP   0.006   0.001   TRP A 351 
 HIS   0.011   0.001   HIS A 784 
 PHE   0.017   0.002   PHE A1390 
 TYR   0.022   0.002   TYR A1893 
 ARG   0.004   0.000   ARG A 827 

Details of bonding type rmsd
  hydrogen bonds       : bond        0.04644 ( 1226)
  hydrogen bonds       : angle       3.42440 ( 3672)
  covalent geometry    : bond        0.00341 (12826)
  covalent geometry    : angle       0.54395 (17400)

*********************** REFINEMENT MACRO_CYCLE 6 OF 10 ************************

-------------------------------------------------------------------------------
Update Rama plot phi/psi targets (oldfield only)
************************************************

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Re-set Ramachandran plot restraints
***********************************
  favored: oldfield
  allowed: oldfield
  outlier: oldfield

  3490 Ramachandran restraints generated.
    1745 Oldfield, 0 Emsley, 1745 emsley8k and 0 Phi/Psi/2.
-------------------------------------------------------------------------------
Optimize residue side-chains
****************************
Evaluate side-chains
  137 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 23
    poor density    : 114
  time to evaluate  : 1.389 
Fit side-chains
   revert: symmetry clash
REVERT: A  670 GLN cc_start: 0.7861 (pp30) cc_final: 0.7546 (pp30)
REVERT: A 1397 MET cc_start: 0.7056 (OUTLIER) cc_final: 0.6763 (mtm)
REVERT: A 1431 LYS cc_start: 0.8090 (tttt) cc_final: 0.6827 (tmtt)
REVERT: A 1455 GLN cc_start: 0.7791 (OUTLIER) cc_final: 0.7504 (mt0)
REVERT: A 1531 GLU cc_start: 0.7432 (tp30) cc_final: 0.6905 (tp30)
REVERT: A 1644 ILE cc_start: 0.7027 (mt) cc_final: 0.6512 (mp)
  outliers start: 23
  outliers final: 16
  residues processed: 128
  average time/residue: 1.3745
  time to fit residues: 191.9985
Evaluate side-chains
  132 residues out of total 1324 (non-[ALA,GLY,PRO]) need to be fit.
    rotamer outliers: 18
    poor density    : 114
  time to evaluate  : 1.393 
-------------------------------------------------------------------------------
Set rotamer restraints
**********************
Chi-restraints excluded: chain A residue  314 LEU
Chi-restraints excluded: chain A residue  330 LEU
Chi-restraints excluded: chain A residue  337 VAL
Chi-restraints excluded: chain A residue  380 VAL
Chi-restraints excluded: chain A residue  384 GLU
Chi-restraints excluded: chain A residue  784 HIS
Chi-restraints excluded: chain A residue  790 THR
Chi-restraints excluded: chain A residue 1200 VAL
Chi-restraints excluded: chain A residue 1316 LEU
Chi-restraints excluded: chain A residue 1345 SER
Chi-restraints excluded: chain A residue 1384 ILE
Chi-restraints excluded: chain A residue 1397 MET
Chi-restraints excluded: chain A residue 1449 VAL
Chi-restraints excluded: chain A residue 1455 GLN
Chi-restraints excluded: chain A residue 1520 THR
Chi-restraints excluded: chain A residue 1573 SER
Chi-restraints excluded: chain A residue 1650 SER
Chi-restraints excluded: chain A residue 1905 VAL
Rotamers are restrained with sigma=2.50
-------------------------------------------------------------------------------
XYZ refinement
**************
Traceback (most recent call last):
  File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/build/../modules/phenix/phenix/command_line/real_space_refine.py", line 8, in <module>
    run_program(real_space_refine.Program)
  File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/cctbx_project/iotbx/cli_parser.py", line 994, in run_program
    task.run()
  File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/phenix/phenix/programs/real_space_refine.py", line 210, in run
    log         = self.logger)
  File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/modules/phenix/phenix/refinement/rsr/wrappers.py", line 58, in __init__
    log                   = log)
  File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/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-5710/modules/phenix/phenix/refinement/macro_cycle_real_space.py", line 321, in caller
    func()
  File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/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-5710/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-5710/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-5710/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-5710/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-5710/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-5710/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-5710/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-5710/modules/cctbx_project/scitbx/lbfgs/__init__.py", line 278, in run
    line_search)
  File "/net/cci-filer2/raid1/xp/phenix/phenix-dev-5710/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-5710/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-5710/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-5710/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-5710/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: 58.7438 > 50:
  distance: 56 - 78:  13.252
  distance: 60 - 89:  18.894
  distance: 68 - 94:  10.071
  distance: 73 - 78:  15.718
  distance: 74 - 103:  12.679
  distance: 78 - 79:   7.765
  distance: 79 - 80:  12.531
  distance: 79 - 82:   7.972
  distance: 80 - 81:   3.705
  distance: 80 - 89:   5.868
  distance: 81 - 109:  17.710
  distance: 82 - 83:  11.474
  distance: 83 - 84:  12.216
  distance: 84 - 85:   8.300
  distance: 85 - 86:   8.328
  distance: 86 - 87:   4.774
  distance: 86 - 88:  11.791
  distance: 89 - 90:   5.012
  distance: 90 - 91:   4.069
  distance: 90 - 93:  24.107
  distance: 91 - 92:  11.741
  distance: 91 - 94:  17.182
  distance: 92 - 120:   6.848
  distance: 94 - 95:   4.351
  distance: 95 - 96:  10.738
  distance: 95 - 98:   7.645
  distance: 96 - 97:   9.701
  distance: 96 - 103:  16.505
  distance: 97 - 125:  26.900
  distance: 98 - 99:  18.979
  distance: 99 - 100:   9.249
  distance: 100 - 101:   9.866
  distance: 100 - 102:   9.565
  distance: 103 - 104:  10.576
  distance: 104 - 105:  17.631
  distance: 104 - 107:   9.765
  distance: 105 - 106:   5.035
  distance: 105 - 109:  17.174
  distance: 106 - 130:  21.349
  distance: 107 - 108:  19.478
  distance: 109 - 110:  12.720
  distance: 110 - 111:   3.544
  distance: 110 - 113:  13.224
  distance: 111 - 112:  24.209
  distance: 111 - 120:  19.095
  distance: 113 - 114:  21.396
  distance: 114 - 115:  25.465
  distance: 115 - 116:  10.301
  distance: 116 - 117:  12.719
  distance: 117 - 118:  12.350
  distance: 117 - 119:  11.000
  distance: 120 - 121:   8.893
  distance: 121 - 122:  15.663
  distance: 121 - 124:   4.495
  distance: 122 - 123:  12.867
  distance: 122 - 125:  22.042
  distance: 125 - 126:  21.465
  distance: 126 - 127:   9.435
  distance: 126 - 129:  15.543
  distance: 127 - 128:  17.866
  distance: 127 - 130:  21.838
  distance: 130 - 131:   8.645
  distance: 131 - 132:  19.037
  distance: 131 - 134:  15.494
  distance: 132 - 133:  42.166
  distance: 132 - 137:  16.057
  distance: 134 - 135:  23.313
  distance: 134 - 136:  23.997
  distance: 138 - 139:  24.307
  distance: 138 - 141:  41.357
  distance: 139 - 140:  39.502
  distance: 139 - 142:  23.405
  distance: 140 - 163:  31.294
  distance: 142 - 143:  10.772
  distance: 142 - 148:   7.282
  distance: 143 - 144:  14.035
  distance: 143 - 146:   7.784
  distance: 144 - 145:  15.259
  distance: 144 - 149:   3.241
  distance: 145 - 172:  14.900
  distance: 146 - 147:  21.596
  distance: 147 - 148:  22.017