Protein Search

James Knight knight at quad.cs.ucdavis.edu
Sun Jun 30 20:59:21 EST 1996

If you pick up the SEQIO package, whose new version I just
released, then the following program will do what you want
(except that you need to fill in the actual sequence
analysis part):

#include <stdio.h>
#include <stdlib.h>
#include "seqio.h"

int main(int argc, char *argv[])
  int i, len;
  char *seq;
  SEQFILE *sfp;

  for (i=1; i < argc; i++) {
    if ((sfp = seqfopen2(argv[i])) == NULL)

    while ((seq = seqfgetseq(sfp, &len, 0)) != NULL) {
       * Fill in the sequence analysis on "seq", whose
       * length is "len", here.


  return 0;

Once compiled with the package, this program can read
one or more files, or read one, a couple or all of the
entries from the Swiss-Prot database.


More information about the Bio-soft mailing list

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