Hello chimerax-users,

I am comparing different structures that all look mostly identical, and I have to generate figures of some interesting place, but they of course don't all have the same residue numbers and chain IDs (would be too easy...). I managed to work around this problem this way:

- open reference structure in which I know the chain IDs and residue numbers I want to look at

- open another structure, align it to the reference with mmaker

- use 'select zone' with option 'residues true' and an appropriate distance cutoff from the reference structure to select the equivalent residues in the structure under scrutiny: works nicely, I simply have to tweak the distance cutoff interactively depending on how different the current structure is compared to the reference, but it is easier than looking up residue numbers and chain IDs and keeping track of all this different info per structure

I am missing one thing: I would like to save these selections as named selections, so I can easily select them again further down in the script. So I need a way to say 'name someName sel' except I can't use that last 'sel' keyword because it is never expanded to the actual atom-spec that would make this trick work, and instead always points to whatever is currently selected.

Is there a way to do this? I looked at 'info sel level residue' that prints to the Log something close to what I need (although more than just the atom-spec). Ideally I would like this to work from a command script, so I can automate the generation of the picture of that interesting location from many structures.

Thank you in advance,