#!/usr/bin/env python
# coding: utf-8

# In[ ]:


from pymol import cmd
cmd.fetch("1a02")
cmd.remove("not (chain F+J or polymer.nucleic)")
cmd.select("set1"," resn DA+DT+DC+DG and name O3'+O4'+O5'")
cmd.select("set2"," resn DA+DT+DC+DG and name OP1+OP2")
cmd.select("set3"," resn DA+DT+DC+DG and name N1+N2+N3+N4+N6+N7+N9")

