last_thresh.c

序列对齐 Compare a protein sequence to a protein sequence database or a DNA sequence to a DNA sequenc

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#include "defs.h"
#include "mm_file.h"

#include "structs.h"
#include "param.h"

#ifndef PCOMPLIB
#include "mw.h"
#else
#include "msg.h"
#include "p_mw.h"

void do_stage2(struct beststr **bptr, int nbest, struct mngmsg m_msg0,
	       int s_func, struct qmng_str *qm_msp);
#endif

static char thresh_str[MAX_STR];

int E1_to_s(double e_val, int n0, int n1, int db_size, 	void *pu);

int
last_calc(
#ifndef PCOMPLIB
	  unsigned char **aa0, unsigned char *aa1, int maxn,
#endif	  
	  struct beststr **bptr, int nbest, 
	  struct mngmsg m_msg, struct pstruct *ppst
#ifdef PCOMPLIB
	  , struct qmng_str *qm_msp
#else
	  , void **f_str
#endif
	  , void *pstat_str)
{

  if (ppst->zdb_size < 0 ) ppst->zdb_size = m_msg.db.entries;
  ppst->repeat_thresh = E1_to_s(ppst->e_cut, m_msg.n0, bptr[0]->seq->n1, ppst->zdb_size, pstat_str);

  ppst->other_info = thresh_str;
  sprintf(thresh_str,"Threshold: E() < %.2g score: %d\n",ppst->e_cut, ppst->repeat_thresh);

  return nbest;
}