File: lba.c

package info (click to toggle)
wise 2.4.1-21
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, buster, sid
  • size: 27,140 kB
  • sloc: ansic: 276,365; makefile: 1,003; perl: 886; lex: 93; yacc: 81; sh: 24
file content (120 lines) | stat: -rw-r--r-- 3,677 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
#include "version.h"
#include "largeblockdp.h"
#include "seqaligndisplay.h"


char * program_name = "lba";

void show_version(FILE * ofp)
{
  fprintf(ofp,"%s\nVersion: %s\nReleased: %s\nCompiled: %s\n",program_name,VERSION_NUMBER,RELEASE_DAY,COMPILE_DATE);
  fprintf(ofp,"\nThis program is freely distributed under a GNU General Public License\n");
  fprintf(ofp,"The source code is copyright (c) EMBL and others\n");
  fprintf(ofp,"There is no warranty, implied or otherwise on the performance of this program\n");
  fprintf(ofp,"For more information read the GNULICENSE file in the distribution\n\n");
  fprintf(ofp,"Credits: Ewan Birney <birney@ebi.ac.uk>\n");
  exit(63);   
}


void show_help(FILE * ofp)
{
  fprintf(ofp,"%s query-seq target-seq\n",program_name);
  fprintf(ofp,"Parameters (all as probabilities)\n");
  fprintf(ofp,"  -matchp [0.75] probability of match in block\n");
  fprintf(ofp,"  -gap    [0.1]  probability of opening a gap in a block\n");
  fprintf(ofp,"  -block_open [0.0005] probability of a block occuring\n");
  fprintf(ofp,"  -block_ext [0.8] probability ratio of block extension vs null\n");
  fprintf(ofp,"  -undual    [1.0] probability ratio of un-matched extension vs null, both seqs\n");
  fprintf(ofp,"  -unsingle   [1.0] probability ratio of un-matched extension vs null, one seq\n");
  fprintf(ofp,"Output format\n");
  fprintf(ofp,"  -[no]alb    show alb format (default no)\n");
  fprintf(ofp,"  -[no]pretty show pretty format (default no)\n");


  show_help_DPRunImpl(ofp);


  show_standard_options(ofp);

}



int main(int argc,char ** argv)
{
  Sequence * one;
  Sequence * two;
  PackAln  * pal;
  AlnBlock * alb;
  DnaProbMatrix * dmp;
  DnaMatrix * mat;

  Probability match    = 0.75;
  Probability gap_open = 0.1;
  Probability gap_ext  = 0.3;
  Probability block_open = 0.0005;
  Probability un_dual  = 1.0;
  Probability un_single = 1.0;
  Probability real_ext  = 0.6;

  boolean show_pretty = TRUE;
  boolean show_alb    = FALSE;

  ComplexSequence *cone, *ctwo;  
  ComplexSequenceEvalSet *cses;  

  DPRunImpl * dpri;


  dpri = new_DPRunImpl_from_argv(&argc,argv);

  strip_out_float_argument(&argc,argv,"matchp",&match);
  strip_out_float_argument(&argc,argv,"gap",&gap_open);
  strip_out_float_argument(&argc,argv,"ext",&gap_ext);
  strip_out_float_argument(&argc,argv,"block_open",&block_open);
  strip_out_float_argument(&argc,argv,"block_ext",&real_ext);
  strip_out_float_argument(&argc,argv,"undual",&un_dual);
  strip_out_float_argument(&argc,argv,"unsingle",&un_single);

  strip_out_boolean_def_argument(&argc,argv,"pretty",&show_pretty);
  strip_out_boolean_def_argument(&argc,argv,"alb",&show_alb);

  strip_out_standard_options(&argc,argv,show_help,show_version);


  if( argc != 3 ) {
    show_help(stdout);
    exit(12);
  }

  one = read_fasta_file_Sequence(argv[1]);
  two = read_fasta_file_Sequence(argv[2]);

  dmp = DnaProbMatrix_from_match(match,NMaskType_BANNED);
  flat_null_DnaProbMatrix(dmp);  
  mat = DnaMatrix_from_DnaProbMatrix(dmp);  

  cses = default_DNA_ComplexSequenceEvalSet();  
  cone = new_ComplexSequence(one,cses);  
  ctwo = new_ComplexSequence(two,cses);  


  pal = PackAln_bestmemory_LargeBlockAligner(cone,ctwo,mat,
					     Probability2Score(real_ext),
					     Probability2Score(block_open),
					     Probability2Score(un_dual),
					     Probability2Score(un_single),
					     Probability2Score(gap_open),
					     Probability2Score(gap_ext),NULL,dpri);

  alb = convert_PackAln_to_AlnBlock_LargeBlockAligner(pal);

  if( show_pretty == TRUE )
    write_pretty_seq_align(alb,one,two,12,60,stdout);

  if( show_alb == TRUE )
    dump_ascii_AlnBlock(alb,stdout);

 
}