PyMOL>load ./selSi.pdb,obj ObjectMolReadPDBStr: read MODEL 1 ObjectMolReadPDBStr: read MODEL 2 ObjectMolReadPDBStr: read MODEL 3 ObjectMolReadPDBStr: read MODEL 4 ObjectMolReadPDBStr: read MODEL 5 ObjectMolReadPDBStr: read MODEL 6 ObjectMolReadPDBStr: read MODEL 7 ObjectMolReadPDBStr: read MODEL 8 ObjectMolReadPDBStr: read MODEL 9 ObjectMolReadPDBStr: read MODEL 10 ObjectMolReadPDBStr: read MODEL 11 ObjectMolReadPDBStr: read MODEL 12 ObjectMolReadPDBStr: read MODEL 13 ObjectMolReadPDBStr: read MODEL 14 ObjectMolReadPDBStr: read MODEL 15 ObjectMolReadPDBStr: read MODEL 16 ObjectMolReadPDBStr: read MODEL 17 ObjectMolReadPDBStr: read MODEL 18 ObjectMolReadPDBStr: read MODEL 19 ObjectMolReadPDBStr: read MODEL 20 CmdLoad: "./selSi.pdb" loaded as "obj". PyMOL>hide everything, all PyMOL>dss PyMOL>show cartoon, all PyMOL>set bg_rgb, [1,1,1] Setting: bg_rgb set to [ 1.00000, 1.00000, 1.00000 ]. PyMOL>color red, ss h Executive: Colored 92 atoms. PyMOL>color cyan,ss s PyMOL>color gray,ss ""+l Executive: Colored 2586 atoms. PyMOL>set cartoon_highlight_color,yellow Setting: cartoon_highlight_color set to yellow. PyMOL>orient PyMOL>set ray_trace_mode, 3 Setting: ray_trace_mode set to 3. PyMOL>set antialias, 2 Setting: antialias set to 2. PyMOL>ray 600,600 Ray: render time: 1.59 sec. = 2267.2 frames/hour (1.59 sec. accum.). PyMOL>png molecule.png ScenePNG: wrote 600x600 pixel image to file "molecule.png". PyMOL>set stereo,1 Setting: stereo set to on. PyMOL>set all_states,1 Setting: all_states set to on. PyMOL>stereo walleye Setting: stereo_mode set to 3. PyMOL>select bb,n. c+ca+n Selector: selection "bb" defined with 480 atoms. PyMOL>hide everything,all PyMOL>show lines,bb PyMOL>color blue,bb Executive: Colored 480 atoms. PyMOL>set line_width,0.5 Setting: line_width set to 0.50000. PyMOL>set bg_rgb, [1,1,1] Setting: bg_rgb set to [ 1.00000, 1.00000, 1.00000 ]. PyMOL>zoom complete=1 PyMOL>set ray_trace_mode, 0 Setting: ray_trace_mode set to 0. PyMOL>ray 1200,600 Ray: render time: 6.62 sec. = 544.0 frames/hour (8.21 sec. accum.). PyMOL>png stereo.png ScenePNG: wrote 1200x600 pixel image to file "stereo.png". PyMOL>run /data/PSVS/PSVS-v1.4-dev/strcalc//average3d.py PyMOL>from average3d import avgStates Traceback (most recent call last): File "/farm/software/pymol-1.1r1-x64/modules/pymol/parser.py", line 455, in parse exec(layer.com2+"\n",self.pymol_names,self.pymol_names) File "", line 1, in ? ImportError: No module named average3d PyMOL>set all_states,1 Setting: all_states set to on. PyMOL>set stereo,0 Setting: stereo set to off. PyMOL>hide everything, all PyMOL>avgStates('obj','name CA',1,0,'yesfit_ALL','yes',1,0,0) -------------------------------------------------------------------------------- Averaging 20 states [1, 20] of object "obj and name CA" (20 total states) to get new single-state object "yesfit_ALL"... -------------------------------------------------------------------------------- -> proceeding WITH FITTING to reference state... intrafit_rmsds = [-1.0, 4.6467714309692383, 12.263684272766113, 9.1835832595825195, 6.3323593139648438, 10.615207672119141, 4.620795726776123, 6.0782604217529297, 7.0329079627990723, 5.9390196800231934, 4.3088283538818359, 13.932190895080566, 3.8001627922058105, 5.9407286643981934, 7.1109766960144043, 6.6996183395385742, 6.805389404296875, 6.3814840316772461, 10.017776489257812, 4.6687173843383789] -> atom index_map = [None, 2, 21, 32, 47, 66, 86, 97, 113, 127, 148, 168, 187, 202, 216, 235, 257, 274, 291, 310, 325, 342, 356, 373, 384, 399, 411, 433, 452, 469, 479, 496, 510, 521, 542, 563, 587, 598, 614, 630, 641, 655, 674, 690, 707, 719, 736, 755, 769, 791, 805, 815, 831, 847, 866, 888, 912, 931, 948, 965, 984, 996, 1011, 1021, 1042, 1056, 1078, 1094, 1116, 1140, 1147, 1162, 1173, 1195, 1214, 1229, 1246, 1263, 1280, 1297, 1314, 1332, 1351, 1362, 1377, 1396, 1416, 1427, 1443, 1457, 1478, 1498, 1517, 1532, 1546, 1565, 1587, 1604, 1621, 1640, 1655, 1672, 1686, 1703, 1714, 1729, 1741, 1763, 1782, 1799, 1809, 1826, 1840, 1851, 1872, 1893, 1917, 1928, 1944, 1960, 1971, 1985, 2004, 2020, 2037, 2049, 2066, 2085, 2099, 2121, 2135, 2145, 2161, 2177, 2196, 2218, 2242, 2261, 2278, 2295, 2314, 2326, 2341, 2351, 2372, 2386, 2408, 2424, 2446, 2470, 2477, 2492, 2503, 2525, 2544, 2559, 2576, 2593, 2610, 2627, 2644] "obj and name CA" index = 1 rmsd = 7.868 "obj and name CA" index = 2 rmsd = 6.934 "obj and name CA" index = 3 rmsd = 6.326 "obj and name CA" index = 4 rmsd = 5.359 "obj and name CA" index = 5 rmsd = 5.728 "obj and name CA" index = 6 rmsd = 5.127 "obj and name CA" index = 7 rmsd = 3.969 "obj and name CA" index = 8 rmsd = 3.720 "obj and name CA" index = 9 rmsd = 3.896 "obj and name CA" index = 10 rmsd = 3.036 "obj and name CA" index = 11 rmsd = 2.569 "obj and name CA" index = 12 rmsd = 3.194 "obj and name CA" index = 13 rmsd = 3.065 "obj and name CA" index = 14 rmsd = 2.541 "obj and name CA" index = 15 rmsd = 3.052 "obj and name CA" index = 16 rmsd = 3.538 "obj and name CA" index = 17 rmsd = 3.449 "obj and name CA" index = 18 rmsd = 3.606 "obj and name CA" index = 19 rmsd = 4.511 "obj and name CA" index = 20 rmsd = 5.132 "obj and name CA" index = 21 rmsd = 5.375 "obj and name CA" index = 22 rmsd = 6.017 "obj and name CA" index = 23 rmsd = 5.674 "obj and name CA" index = 24 rmsd = 4.649 "obj and name CA" index = 25 rmsd = 4.518 "obj and name CA" index = 26 rmsd = 4.166 "obj and name CA" index = 27 rmsd = 4.112 "obj and name CA" index = 28 rmsd = 4.169 "obj and name CA" index = 29 rmsd = 3.414 "obj and name CA" index = 30 rmsd = 3.030 "obj and name CA" index = 31 rmsd = 3.728 "obj and name CA" index = 32 rmsd = 3.614 "obj and name CA" index = 33 rmsd = 3.033 "obj and name CA" index = 34 rmsd = 3.291 "obj and name CA" index = 35 rmsd = 3.778 "obj and name CA" index = 36 rmsd = 3.593 "obj and name CA" index = 37 rmsd = 3.806 "obj and name CA" index = 38 rmsd = 4.397 "obj and name CA" index = 39 rmsd = 4.737 "obj and name CA" index = 40 rmsd = 4.780 "obj and name CA" index = 41 rmsd = 5.457 "obj and name CA" index = 42 rmsd = 6.159 "obj and name CA" index = 43 rmsd = 6.221 "obj and name CA" index = 44 rmsd = 6.461 "obj and name CA" index = 45 rmsd = 7.113 "obj and name CA" index = 46 rmsd = 7.297 "obj and name CA" index = 47 rmsd = 7.975 "obj and name CA" index = 48 rmsd = 8.136 "obj and name CA" index = 49 rmsd = 7.825 "obj and name CA" index = 50 rmsd = 7.952 "obj and name CA" index = 51 rmsd = 7.256 "obj and name CA" index = 52 rmsd = 6.414 "obj and name CA" index = 53 rmsd = 6.505 "obj and name CA" index = 54 rmsd = 6.375 "obj and name CA" index = 55 rmsd = 5.444 "obj and name CA" index = 56 rmsd = 5.040 "obj and name CA" index = 57 rmsd = 5.305 "obj and name CA" index = 58 rmsd = 4.577 "obj and name CA" index = 59 rmsd = 3.972 "obj and name CA" index = 60 rmsd = 4.286 "obj and name CA" index = 61 rmsd = 3.873 "obj and name CA" index = 62 rmsd = 2.950 "obj and name CA" index = 63 rmsd = 3.252 "obj and name CA" index = 64 rmsd = 3.656 "obj and name CA" index = 65 rmsd = 2.844 "obj and name CA" index = 66 rmsd = 2.660 "obj and name CA" index = 67 rmsd = 3.819 "obj and name CA" index = 68 rmsd = 4.003 "obj and name CA" index = 69 rmsd = 3.735 "obj and name CA" index = 70 rmsd = 4.427 "obj and name CA" index = 71 rmsd = 5.950 "obj and name CA" index = 72 rmsd = 7.041 "obj and name CA" index = 73 rmsd = 8.524 "obj and name CA" index = 74 rmsd = 8.957 "obj and name CA" index = 75 rmsd = 9.850 "obj and name CA" index = 76 rmsd = 10.545 "obj and name CA" index = 77 rmsd = 11.879 "obj and name CA" index = 78 rmsd = 12.463 "obj and name CA" index = 79 rmsd = 13.951 "obj and name CA" index = 80 rmsd = 15.177 "obj and name CA" index = 81 rmsd = 7.869 "obj and name CA" index = 82 rmsd = 6.934 "obj and name CA" index = 83 rmsd = 6.326 "obj and name CA" index = 84 rmsd = 5.359 "obj and name CA" index = 85 rmsd = 5.728 "obj and name CA" index = 86 rmsd = 5.127 "obj and name CA" index = 87 rmsd = 3.969 "obj and name CA" index = 88 rmsd = 3.721 "obj and name CA" index = 89 rmsd = 3.896 "obj and name CA" index = 90 rmsd = 3.036 "obj and name CA" index = 91 rmsd = 2.569 "obj and name CA" index = 92 rmsd = 3.193 "obj and name CA" index = 93 rmsd = 3.065 "obj and name CA" index = 94 rmsd = 2.541 "obj and name CA" index = 95 rmsd = 3.052 "obj and name CA" index = 96 rmsd = 3.538 "obj and name CA" index = 97 rmsd = 3.449 "obj and name CA" index = 98 rmsd = 3.606 "obj and name CA" index = 99 rmsd = 4.511 "obj and name CA" index = 100 rmsd = 5.132 "obj and name CA" index = 101 rmsd = 5.375 "obj and name CA" index = 102 rmsd = 6.017 "obj and name CA" index = 103 rmsd = 5.674 "obj and name CA" index = 104 rmsd = 4.649 "obj and name CA" index = 105 rmsd = 4.518 "obj and name CA" index = 106 rmsd = 4.166 "obj and name CA" index = 107 rmsd = 4.112 "obj and name CA" index = 108 rmsd = 4.169 "obj and name CA" index = 109 rmsd = 3.414 "obj and name CA" index = 110 rmsd = 3.030 "obj and name CA" index = 111 rmsd = 3.728 "obj and name CA" index = 112 rmsd = 3.614 "obj and name CA" index = 113 rmsd = 3.033 "obj and name CA" index = 114 rmsd = 3.291 "obj and name CA" index = 115 rmsd = 3.778 "obj and name CA" index = 116 rmsd = 3.593 "obj and name CA" index = 117 rmsd = 3.806 "obj and name CA" index = 118 rmsd = 4.397 "obj and name CA" index = 119 rmsd = 4.737 "obj and name CA" index = 120 rmsd = 4.780 "obj and name CA" index = 121 rmsd = 5.456 "obj and name CA" index = 122 rmsd = 6.159 "obj and name CA" index = 123 rmsd = 6.221 "obj and name CA" index = 124 rmsd = 6.461 "obj and name CA" index = 125 rmsd = 7.113 "obj and name CA" index = 126 rmsd = 7.297 "obj and name CA" index = 127 rmsd = 7.975 "obj and name CA" index = 128 rmsd = 8.136 "obj and name CA" index = 129 rmsd = 7.825 "obj and name CA" index = 130 rmsd = 7.952 "obj and name CA" index = 131 rmsd = 7.255 "obj and name CA" index = 132 rmsd = 6.415 "obj and name CA" index = 133 rmsd = 6.505 "obj and name CA" index = 134 rmsd = 6.375 "obj and name CA" index = 135 rmsd = 5.443 "obj and name CA" index = 136 rmsd = 5.040 "obj and name CA" index = 137 rmsd = 5.305 "obj and name CA" index = 138 rmsd = 4.577 "obj and name CA" index = 139 rmsd = 3.972 "obj and name CA" index = 140 rmsd = 4.286 "obj and name CA" index = 141 rmsd = 3.873 "obj and name CA" index = 142 rmsd = 2.950 "obj and name CA" index = 143 rmsd = 3.252 "obj and name CA" index = 144 rmsd = 3.655 "obj and name CA" index = 145 rmsd = 2.844 "obj and name CA" index = 146 rmsd = 2.660 "obj and name CA" index = 147 rmsd = 3.819 "obj and name CA" index = 148 rmsd = 4.004 "obj and name CA" index = 149 rmsd = 3.735 "obj and name CA" index = 150 rmsd = 4.427 "obj and name CA" index = 151 rmsd = 5.950 "obj and name CA" index = 152 rmsd = 7.041 "obj and name CA" index = 153 rmsd = 8.524 "obj and name CA" index = 154 rmsd = 8.957 "obj and name CA" index = 155 rmsd = 9.850 "obj and name CA" index = 156 rmsd = 10.545 "obj and name CA" index = 157 rmsd = 11.879 "obj and name CA" index = 158 rmsd = 12.463 "obj and name CA" index = 159 rmsd = 13.951 "obj and name CA" index = 160 rmsd = 15.177 PyMOL>cartoon putty, yesfit_ALL and name CA PyMOL>show cartoon, yesfit_ALL and name CA PyMOL>color red, ss h Executive: Colored 276 atoms. PyMOL>color cyan,ss s PyMOL>color gray,ss ""+l Executive: Colored 5222 atoms. PyMOL>set cartoon_highlight_color,yellow Setting: cartoon_highlight_color set to yellow. PyMOL>set bg_rgb, [1,1,1] Setting: bg_rgb set to [ 1.00000, 1.00000, 1.00000 ]. PyMOL>set ray_trace_mode, 3 Setting: ray_trace_mode set to 3. PyMOL>ray 600,600 Ray: render time: 1.90 sec. = 1898.8 frames/hour (10.10 sec. accum.). PyMOL>png sausage.png ScenePNG: wrote 600x600 pixel image to file "sausage.png". PyMOL>quit