Mercurial > repos > siyuan > prada
comparison pyPRADA_1.2/tools/bwa-0.5.7-mh/bwtgap.h @ 0:acc2ca1a3ba4
Uploaded
author | siyuan |
---|---|
date | Thu, 20 Feb 2014 00:44:58 -0500 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:acc2ca1a3ba4 |
---|---|
1 #ifndef BWTGAP_H_ | |
2 #define BWTGAP_H_ | |
3 | |
4 #include "bwt.h" | |
5 #include "bwtaln.h" | |
6 | |
7 typedef struct { // recursion stack | |
8 u_int32_t info; // score<<21 | a<<20 | i | |
9 u_int32_t n_mm:8, n_gapo:8, n_gape:8, state:2, n_seed_mm:6; | |
10 bwtint_t k, l; // (k,l) is the SA region of [i,n-1] | |
11 int last_diff_pos; | |
12 } gap_entry_t; | |
13 | |
14 typedef struct { | |
15 int n_entries, m_entries; | |
16 gap_entry_t *stack; | |
17 } gap_stack1_t; | |
18 | |
19 typedef struct { | |
20 int n_stacks, best, n_entries; | |
21 gap_stack1_t *stacks; | |
22 } gap_stack_t; | |
23 | |
24 #ifdef __cplusplus | |
25 extern "C" { | |
26 #endif | |
27 | |
28 gap_stack_t *gap_init_stack(int max_mm, int max_gapo, int max_gape, const gap_opt_t *opt); | |
29 void gap_destroy_stack(gap_stack_t *stack); | |
30 bwt_aln1_t *bwt_match_gap(bwt_t *const bwt[2], int len, const ubyte_t *seq[2], bwt_width_t *w[2], | |
31 bwt_width_t *seed_w[2], const gap_opt_t *opt, int *_n_aln, gap_stack_t *stack); | |
32 void bwa_aln2seq(int n_aln, const bwt_aln1_t *aln, bwa_seq_t *s); | |
33 | |
34 #ifdef __cplusplus | |
35 } | |
36 #endif | |
37 | |
38 #endif |