load 3QTP.pdb, 3QTP load 3qtp.omap, el-dens remove chain B remove resn HOH remove resn 2PG remove resn SO4 center resi 424 select chain A, chain A set line_width, 10 color orange isomesh new_map, el-dens, 1, chain A and (name ca+c+n+o), carve=2