IUBio Biosequences .. Software .. Molbio soft .. Network News .. FTP

Automatic Extraction/Joining of Exons

Paula Burch pburch at bcm.tmc.edu
Thu Jul 25 17:49:53 EST 1996

John Watson <watson_j at bms.com> writes:
>Does anyone know of a utility or script that will read a genbank-type 
>file and automatically extract the coding regions as listed in the 

I have a Unix script, called gb2pep, that creates a peptide sequence
file from the peptide sequence included in the comments of a genbank
file that's been turned into gcg format with 'fromgenbank'. It was
written for us by Scott Rose, then at GCG, using gawk; you have to
'reformat' the results before using them.  

Mike Hogan told me that we can redistribute the script. He says that
you may prefer the GCG program 'lookup', which can extract the coding
regions of the DNA sequence, ready for translation. I like gb2pep
because it's so simple and obvious to use. The script is available on
the web from

We're not providing any support for it, of course, but it's extremely
easy to use, in our experience.

Paula E. Burch, Ph.D.           Molecular Biology Computational Resource
Baylor College of Medicine      phone: (713)798-6023  fax: (713)798-4279
Houston, Texas 77030                        internet: pburch at bcm.tmc.edu
http://mbcr.bcm.tmc.edu/pburch.html                 www at mbcr.bcm.tmc.edu 

More information about the Info-gcg mailing list

Send comments to us at biosci-help [At] net.bio.net