# ifndef __SEQ_H__
# define __SEQ_H__
#ifdef __cplusplus
extern "C" {
#endif

# include "util.h"

/* Maximum number of positions that   */
/* can be selected in a "globin style" selection.  */

# define MAX_SEQ_SELECTION 300

# define SEQ_BADVOL_FLAG 1 

typedef struct {
  char aa[600] ;
  char id[32] ;
  int  access ;
  double vol;
  double wt ;
  double VolDev, NatomsDev;
  int natoms ;
  int flag ;
  char * struc ;
  char chain;
  float homology ;
} Sequence  ;

/* # define SEQ_START_POS 8 */
typedef struct {
  Sequence * seq[2000] ;
  int nseq, SeqSz,LengthSeg1, SeqStartPos;
  double nselect ;
  StringIndex * SI ;
  Statistic vol, natoms;
  Statistic AVarWt, VVarWt;
  Statistic IndVol ;
  char * HomolVs ;
  char * dVVs ;
} Sequences ;

# define NUM_AA 30
typedef struct {
  StringIndex * SI ;
  char aa[NUM_AA];
  char aa3[NUM_AA][7];
  double vol[NUM_AA];
  int natoms[NUM_AA];
  double frac[NUM_AA];
} VolumeLookupTable ;

# define IS_AA(c) (\
		   (c == 'B') || \
		   (c == 'Z') || \
		   (c == 'G') || \
		   (c == 'A') || \
		   (c == 'V') || \
		   (c == 'L') || \
		   (c == 'I') || \
		   (c == 'P') || \
		   (c == 'M') || \
		   (c == 'F') || \
		   (c == 'Y') || \
		   (c == 'W') || \
		   (c == 'S') || \
		   (c == 'T') || \
		   (c == 'N') || \
		   (c == 'Q') || \
		   (c == 'C') || \
		   (c == 'H') || \
		   (c == 'E') || \
		   (c == 'D') || \
		   (c == 'R') || \
		   (c == 'K'))

# include "SetOfInt.h"
# include "Seq.prototypes.h"

#ifdef __cplusplus
}
#endif
#endif /* !__SEQ_H__ */