## Not run: 
# ## Read two pdbs
# a.pdb <- read.pdb("1bg2")
# b.pdb <- read.pdb("1goj")
# 
# seqs <- seqbind(aa321(a.pdb$atom[a.pdb$calpha,"resid"]),
#                 aa321(b.pdb$atom[b.pdb$calpha,"resid"]))
# 
# # seqaln(seqs)
# ## End(Not run)
Run the code above in your browser using DataLab