1 #include "hsp.h"
2
3
4
5
6
main(int argc,char ** argv)7 int main(int argc,char ** argv)
8 {
9 char * seq1s = "MVNSNQNQNGNSNGHDDDFPQDSITEPEHMRKLFIGGLDYRTTDENLKAHRRKRRR";
10 char * seq2s = "MHKSEAPNEPEQLRKLFIGGLSFETTDESLREHFEQWGTLTDCVVMRDPNSKRSRGFGFV";
11 int query_coord = 33;
12 int target_coord = 16;
13
14 Sequence * seq1;
15 Sequence * seq2;
16 HSP * hsp;
17 CompMat * mat;
18
19
20 seq1 = new_Sequence_from_strings("seq1",seq1s);
21 seq2 = new_Sequence_from_strings("seq2",seq2s);
22
23
24 mat = read_Blast_file_CompMat("blosum62.bla");
25
26
27 hsp = new_HSP(seq1,seq2,query_coord,target_coord,mat,20);
28
29 show_HSP(hsp,80,stdout);
30
31 }
32