Buz Barstow <b...@mac.com> wrote: > Thanks to everyone who replied to my query about extracting an amino > acid sequence from a PDB file! > > Here is a summary of responses of my query; > > 1. Use SwissPDBViewer ... > 6. Use PDBSET (part of CCP4)
t_coffee comes with s perl script "unpack_extract_from_pdb". For instance, in the following it extracts all of the sequences from the entire PDB dataset: nice -15 find /u4/pdb -name '*.pdb' \ -exec ./extract_from_pdb -chain ALL \ -infile '{}' -mode fasta -seq_field ATOM \; > all_pdb_seq.aa Note, this takes a long time! It would probably be significantly faster if a persistent perl was used, instead of restarting the perl script from scratch for each file. Regards, David Mathog mat...@caltech.edu Manager, Sequence Analysis Facility, Biology Division, Caltech ------------------------------------------------------------------------------ Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day trial. Simplify your report design, integration and deployment - and focus on what you do best, core application coding. Discover what's new with Crystal Reports now. http://p.sf.net/sfu/bobj-july _______________________________________________ PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net) Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net