#!/bin/sh

# $Id: getsequence,v 1.2 1991/07/08 07:45:25 schwartz Exp $
# Get sequences given a search key

SEQDIR=${SEQDIR-/import/msd/dat/sequence}

case $# in
1) ;;
*) echo 'Usage: getsequence key' 2>&1 ; exit 1;;
esac

egrep "$1" $SEQDIR/seqasn.index |
 while read line
 do
	eval $line
	getseq $file $pos | ncbitoascii | wrap 70
 done