out: Hversion out: 2 out: Hnewname out: 0 out: Hmovname out: 11 out: bz_census.c in: Tversion in: 0 in: Tstartcmdfile in: 134601792 out: Hnewname out: 1 out: Hbindname out: 134601792 out: Hcurrent out: 1 out: Hmovname out: 7 out: ~~sam~~ out: Hcheck0 out: 1 out: Hack in: Tcheck out: Hcheck out: 1 in: Tack out: Hunlock in: Tstartfile out: Hbindname out: 134711368 out: Hcurrent out: 0 in: 0 out: Hgrow out: 0 out: 2436 out: Hcheck0 out: 0 out: Hcheck0 out: 0 out: Hack in: Trequest in: 0 in: 512 out: Hdata out: 0 out: 0 out: 512 out: #include "util.h" #include "seq.h" #include "bz_all.h" #include "bz_main.h" #include "bz_census.h" in: Tcheck out: Hcheck out: 0 in: Tcheck out: Hcheck out: 0 in: Tack out: Hgrow out: 0 out: 16 out: Hcheck0 out: 1 out: Hack in: Trequest in: 512 in: 512 out: Hdata out: 0 out: 512 out: 512 out: c(census,n,i); } void census_mask_align(align_t *a, int n, uchar *fwd, uchar *rev, census_t *cens in: Tcheck out: Hcheck out: 0 in: Trequest in: 512 in: 512 out: Hdata out: 0 out: 512 out: 0 out: in: Tcheck out: Hcheck out: 0 in: Trequest in: 0 in: 16 out: Hdata out: 1 out: 0 out: 16 out: +. bz_census.c in: Tcheck out: Hcheck out: 1 in: Tack out: Hsetdot out: 16 out: 16 out: Hunlock in: Trequest in: 1024 in: 512 out: Hdata out: 0 out: 1024 out: 512 out: += len; break; case EDIT_OP_REP: //fprintf(stderr, "%d %d -> ", x, y); for (i=x, j=x+len in: Tcheck out: Hcheck out: 0 in: Trequest in: 1536 in: 512 out: Hdata out: 0 out: 1536 out: 512 out: ensus_mask_interval(int n, uchar *fwd, uchar *rev, int a, int z, census_t census[], int k) { int i in: Tcheck out: Hcheck out: 0 in: Trequest in: 2048 in: 388 out: Hdata out: 0 out: 2048 out: 388 out: =n; ++i) { if (census[i] >= k) fprintf(fp, "%d %d\n", i, census[i]); } fprintf(fp, "}\n") in: Tcheck out: Hcheck out: 0 in: Torigin in: 19 out: Horigin out: 1223 in: Torigin in: 19 out: Horigin out: 696 in: Torigin in: 19 out: Horigin out: 302 in: Torigin in: 18 out: Horigin out: 0 in: Torigin in: 18 out: Horigin out: 0 in: Torigin in: 18 out: Horigin out: 0 in: Torigin in: 29 out: Horigin out: 1909 in: Torigin in: 29 out: Horigin out: 1251 in: Torigin in: 29 out: Horigin out: 527 in: Torigin in: 8 out: Horigin out: 364 in: Torigin in: 8 out: Horigin out: 201 in: Torigin in: 8 out: Horigin out: 76 in: Torigin in: 6 out: Horigin out: 1998 in: Torigin in: 6 out: Horigin out: 1932 in: Tworkfile in: 0 in: 0 in: 0 in: Ttype in: 16 in: /inter out: Hsetdot out: 1178 out: 1183 out: Hmoveto out: 1178 out: Hsetpat out: 5 out: inter out: Hunlock in: Torigin in: 2 out: Horigin out: 1140 in: Tworkfile in: 0 in: 1178 in: 1183 in: Ttype in: 23 in: / out: Hsetdot out: 1547 out: 1552 out: Hmoveto out: 1547 out: Hunlock in: Tworkfile in: 0 in: 1547 in: 1552 in: Ttype in: 25 in: / out: Hsetdot out: 1178 out: 1183 out: Hmoveto out: 1178 out: Hunlock in: Tworkfile in: 0 in: 1178 in: 1183 in: Ttype in: 27 in: / out: Hsetdot out: 1547 out: 1552 out: Hmoveto out: 1547 out: Hunlock in: Tworkfile in: 0 in: 2151 in: 2151 in: Ttype in: 2151 in: out: Hdirty out: 0 in: Ttype in: 2152 in: msp_ in: Tworkfile in: 0 in: 2156 in: 2156 in: Ttype in: 29 in: B #inc in: Tcheck out: Hcheck out: 3 in: Trequest in: 512 in: 512 out: Hdata out: 3 out: 512 out: 512 out: len, pos1, pos2; score_t score, cum_score; int filter; } msp_t; typedef struct msp_table { int in: Tcheck out: Hcheck out: 3 in: Trequest in: 1024 in: 450 out: Hdata out: 3 out: 1024 out: 450 out: bt, ss_t ss, int X, int K, int T); void blast_table_free(blast_table_t *bt); void msp_free_table(m in: Tcheck out: Hcheck out: 3 in: Tworkfile in: 3 in: 646 in: 646 in: Tdclick out: Hsetdot out: 643 out: 654 out: Hunlockfile out: 3 in: Tworkfile in: 0 in: 2156 in: 2156 in: Ttype in: 2156 in: table_t *census_intervals() in: Ttype in: 2184 in: { in: Ttype in: 2186 in: } in: Twrite in: 0 out: Hclean out: 0 out: Hgrowdata out: 64 out: 19 out: 19 out: bz_census.c: #2473 out: Hcheck0 out: 1 out: Hack in: Torigin in: 2 out: Horigin out: 64 in: Tcheck out: Hcheck out: 1 in: Tack out: Hsetdot out: 83 out: 83 out: Hunlock in: Tworkfile in: 0 in: 1703 in: 1703 in: Tworkfile in: 0 in: 1873 in: 1873 in: Tworkfile in: 0 in: 1931 in: 1931 in: Tworkfile in: 0 in: 1932 in: 1932 in: Tworkfile in: 0 in: 2006 in: 2006 in: Tworkfile in: 0 in: 2182 in: 2182 in: Ttype in: 2182 in: uchar census[], int n out: Hdirty out: 0 in: Tworkfile in: 0 in: 2206 in: 2206 in: Torigin in: 20 out: Horigin out: 1291 in: Torigin in: 20 out: Horigin out: 734 in: Torigin in: 20 out: Horigin out: 304 in: Torigin in: 32 out: Horigin out: 0 in: Ttype in: 2206 in: in: Ttype in: 2207 in: in: Tworkfile in: 0 in: 2273 in: 2493 in: Tsnarf in: Tworkfile in: 0 in: 2207 in: 2207 in: Tpaste in: 2207 out: Hgrow out: 2207 out: 220 out: Hcheck0 out: 0 out: Hack in: Trequest in: 2207 in: 220 out: Hdata out: 0 out: 2207 out: 220 out: int i=0, j=0, k; int in = 1; for (k=0; k 2) { print_align_header(sf1, sf2, ss, &gs, K, L); print_msp_l in: Tcheck out: Hcheck out: 4 in: Trequest in: 3276 in: 512 out: Hdata out: 4 out: 3276 out: 512 out: *argv[]) { SEQ *sf1, *sf2; uchar *rf1; blast_table_t *bt; gap_scores_t gs; int W, X, Y, K, L, in: Tcheck out: Hcheck out: 4 in: Tworkfile in: 4 in: 3083 in: 3089 in: Tsnarf in: Tsnarf in: Tsnarf in: Tsnarf in: Tsnarf in: Tworkfile in: 4 in: 3160 in: 3160 in: Ttype in: 3160 in: out: Hdirty out: 4 in: Tworkfile in: 4 in: 3161 in: 3161 in: Ttype in: 3161 in: in: Tworkfile in: 4 in: 2941 in: 2941 in: Ttype in: 2941 in: ; in: Tworkfile in: 4 in: 2945 in: 2946 in: Tworkfile in: 4 in: 2945 in: 2945 in: Ttype in: 2945 in: a in: Tworkfile in: 4 in: 2942 in: 2942 in: Ttype in: 2942 in: in: Ttype in: 2943 in: msp_table_t *m; in: Tworkfile in: 4 in: 3103 in: 3103 in: Ttype in: 3103 in: m = in: Tworkfile in: 4 in: 3187 in: 3187 in: Tworkfile in: 4 in: 3107 in: 3107 in: Tcut in: 3106 in: 3107 in: Tcut in: 3105 in: 3106 in: Tcut in: 3104 in: 3105 in: Tcut in: 3103 in: 3104 in: Tworkfile in: 4 in: 3180 in: 3180 in: Tworkfile in: 4 in: 3183 in: 3183 in: Ttype in: 3183 in: m = census_intervals(census, in: Tworkfile in: 4 in: 3124 in: 3137 in: Tsnarf in: Tworkfile in: 4 in: 3211 in: 3211 in: Tpaste in: 3211 out: Hgrowdata out: 3211 out: 13 out: 13 out: SEQ_LEN(sf1), out: Hcheck0 out: 4 out: Hack in: Tcheck out: Hcheck out: 4 in: Tack out: Hsetdot out: 3211 out: 3224 out: Hunlockfile out: 4 in: Tworkfile in: 4 in: 3211 in: 3211 in: Ttype in: 3211 in: in: Tworkfile in: 4 in: 3225 in: 3225 in: Ttype in: 3225 in: mask_thresh); in: Ttype in: 3240 in: in: Tworkfile in: 4 in: 3187 in: 3187 in: Tcut in: 3186 in: 3187 in: Tcut in: 3185 in: 3186 in: Tcut in: 3184 in: 3185 in: Tcut in: 3183 in: 3184 in: Ttype in: 3183 in: print_ in: Tworkfile in: 4 in: 3244 in: 3244 in: Torigin in: 5 out: Horigin out: 2683 in: Trequest in: 2683 in: 81 out: Hdata out: 4 out: 2683 out: 81 out: p->filter = 1 - p->cum_score; msp_compress(mt); msp_sort_pos1(mt); } in: Tcheck out: Hcheck out: 4 in: Torigin in: 5 out: Horigin out: 2492 in: Trequest in: 2492 in: 191 out: Hdata out: 4 out: 2492 out: 191 out: if (chain == 1 || chain == 2) { msp_t *p; (void)msp_make_chain(mt, bz_flags.G, bz_flags.G, M in: Tcheck out: Hcheck out: 4 in: Torigin in: 5 out: Horigin out: 2423 in: Trequest in: 2423 in: 69 out: Hdata out: 4 out: 2423 out: 69 out: { msp_table_t *mt; mt = blast_search(sf1, sf2, bt, sss, X, K, P); in: Tcheck out: Hcheck out: 4 in: Torigin in: 5 out: Horigin out: 2141 in: Trequest in: 2141 in: 282 out: Hdata out: 4 out: 2141 out: 282 out: /* strand1 -- process one sequence from the second file in one orientation */ static void strand1(S in: Tcheck out: Hcheck out: 4 in: Torigin in: 5 out: Horigin out: 2057 in: Trequest in: 2057 in: 84 out: Hdata out: 4 out: 2057 out: 84 out: substitutions*bz_flags.R : -substitutions*scale*ss[(uchar)'A'][(uchar)'A']); } in: Tcheck out: Hcheck out: 4 in: Tworkfile in: 4 in: 2943 in: 2943 in: Tdclick out: Hsetdot out: 2943 out: 2961 out: Hunlockfile out: 4 in: Tsnarf in: Tcut in: 2943 in: 2961 in: Tworkfile in: 4 in: 3015 in: 3015 in: Ttype in: 3015 in: { in: Ttype in: 3017 in: in: Tworkfile in: 4 in: 3067 in: 3067 in: Ttype in: 3067 in: in: Tworkfile in: 4 in: 3091 in: 3091 in: Ttype in: 3091 in: in: Tworkfile in: 4 in: 3172 in: 3172 in: Ttype in: 3172 in: in: Tworkfile in: 4 in: 3231 in: 3231 in: Tworkfile in: 4 in: 3234 in: 3234 in: Ttype in: 3234 in: } in: Twrite in: 4 out: Hclean out: 4 out: Hgrowdata out: 255 out: 17 out: 17 out: bz_main.c: #5860 out: Hcheck0 out: 1 out: Hack in: Tcheck out: Hcheck out: 1 in: Tack out: Hsetdot out: 272 out: 272 out: Hunlock in: Tworkfile in: 4 in: 3235 in: 3235 in: Tworkfile in: 4 in: 3231 in: 3231 in: Tworkfile in: 4 in: 3178 in: 3178 in: Tworkfile in: 4 in: 3231 in: 3231 in: Ttype in: 3231 in: out: Hdirty out: 4 in: Tworkfile in: 4 in: 3232 in: 3232 in: Tdclick out: Hsetdot out: 3170 out: 3233 out: Hunlockfile out: 4 in: Tworkfile in: 4 in: 3232 in: 3232 in: Tworkfile in: 4 in: 3177 in: 3177 in: Tworkfile in: 4 in: 3092 in: 3092 in: Tworkfile in: 4 in: 3016 in: 3016 in: Ttype in: 3016 in: in: Ttype in: 3017 in: int n; in: Tworkfile in: 4 in: 3102 in: 3102 in: Ttype in: 3102 in: n = in: Tworkfile in: 4 in: 3183 in: 3183 in: Ttype in: 3183 in: in: Ttype in: 3184 in: if (n) printf(" x %d\n", n); in: Tworkfile in: 4 in: 3219 in: 3219 in: Tdclick out: Hunlockfile out: 4 in: Tworkfile in: 4 in: 3184 in: 3184 in: Tworkfile in: 4 in: 3217 in: 3217 in: Tdclick out: Hsetdot out: 3217 out: 3280 out: Hunlockfile out: 4 in: Tsnarf in: Tcut in: 3217 in: 3280 in: Tworkfile in: 4 in: 3187 in: 3220 in: Tworkfile in: 4 in: 3187 in: 3194 in: Tsnarf in: Tcut in: 3187 in: 3194 in: Tworkfile in: 4 in: 3227 in: 3227 in: Tdclick out: Hsetdot out: 3216 out: 3231 out: Hunlockfile out: 4 in: Tworkfile in: 4 in: 3191 in: 3191 in: Twrite in: 4 out: Hclean out: 4 out: Hgrowdata out: 272 out: 17 out: 17 out: bz_main.c: #5838 out: Hcheck0 out: 1 out: Hack in: Torigin in: 2 out: Horigin out: 272 in: Tcheck out: Hcheck out: 1 in: Tack out: Hsetdot out: 289 out: 289 out: Hunlock in: Tworkfile in: 0 in: 2445 in: 2445 in: Ttype in: 289 in: /mask_ali out: Hsetdot out: 543 out: 551 out: Hmoveto out: 543 out: Hsetpat out: 8 out: mask_ali out: Hunlock in: Torigin in: 2 out: Horigin out: 530 in: Tworkfile in: 0 in: 534 in: 534 in: Tdclick out: Hsetdot out: 531 out: 535 out: Hunlockfile out: 0 in: Tsnarf in: Tcut in: 531 in: 535 out: Hdirty out: 0 in: Ttype in: 531 in: int in: Tworkfile in: 0 in: 638 in: 638 in: Ttype in: 638 in: in: Ttype in: 639 in: int cnt; in: Tcut in: 648 in: 649 in: Tworkfile in: 0 in: 648 in: 648 in: Tworkfile in: 0 in: 665 in: 665 in: Ttype in: 665 in: in: Tworkfile in: 0 in: 666 in: 666 in: Ttype in: 666 in: in: Ttype in: 667 in: cnt = 0; in: Tworkfile in: 0 in: 1186 in: 1186 in: Ttype in: 1186 in: cnt += in: Tworkfile in: 0 in: 1637 in: 1637 in: Tdclick out: Hsetdot out: 1637 out: 1643 out: Hunlockfile out: 0 in: Tworkfile in: 0 in: 1558 in: 1558 in: Tdclick out: Hsetdot out: 1557 out: 1561 out: Hunlockfile out: 0 in: Tsnarf in: Tcut in: 1557 in: 1561 in: Ttype in: 1557 in: int in: Tworkfile in: 0 in: 1663 in: 1663 in: Ttype in: 1663 in: int cnt; in: Ttype in: 1673 in: in: Ttype in: 1674 in: cnt = 0; in: Tworkfile in: 0 in: 1973 in: 1973 in: Ttype in: 1973 in: in: Ttype in: 1974 in: ++cnt; in: Tworkfile in: 0 in: 1990 in: 1990 in: Ttype in: 1990 in: in: Ttype in: 1991 in: return cnt; in: Tworkfile in: 0 in: 1569 in: 1569 in: Tdclick out: Hsetdot out: 1561 out: 1581 out: Hunlockfile out: 0 in: Tlook in: 1561 in: 1581 out: Hsetdot out: 1193 out: 1213 out: Hmoveto out: 1193 out: Hunlock in: Torigin in: 2 out: Horigin out: 1160 in: Tworkfile in: 0 in: 1435 in: 1435 in: Ttype in: 1435 in: in: Ttype in: 1436 in: return cnt; in: Torigin in: 17 out: Horigin out: 716 in: Torigin in: 17 out: Horigin out: 377 in: Torigin in: 17 out: Horigin out: 76 in: Tworkfile in: 0 in: 1448 in: 1448 in: Ttype in: 299 in: B bz_census.h out: Hnewname out: 5 out: Hmovname out: 11 out: bz_census.h out: Hgrowdata out: 313 out: 16 out: 16 out: -. bz_census.h out: Hcheck0 out: 1 out: Hack in: Tcheck out: Hcheck out: 1 in: Tack out: Hsetdot out: 329 out: 329 out: Hcurrent out: 5 out: Hunlock in: Tstartfile out: Hbindname out: 134734872 out: Hcurrent out: 5 in: 5 out: Hgrow out: 0 out: 425 out: Hcheck0 out: 5 out: Hmoveto out: 0 out: Hunlock in: Trequest in: 0 in: 425 out: Hdata out: 5 out: 0 out: 425 out: #ifndef BZ_CENSUS_H #define BZ_CENSUS_H typedef unsigned char census_t; void msp_census(census_t ce in: Tcheck out: Hcheck out: 5 in: Tworkfile in: 5 in: 212 in: 212 in: Tworkfile in: 5 in: 133 in: 133 in: Tdclick out: Hsetdot out: 132 out: 136 out: Hunlockfile out: 5 in: Tsnarf in: Tcut in: 132 in: 136 out: Hdirty out: 5 in: Ttype in: 132 in: int in: Tworkfile in: 5 in: 229 in: 229 in: Tdclick out: Hsetdot out: 228 out: 232 out: Hunlockfile out: 5 in: Tsnarf in: Tcut in: 228 in: 232 in: Ttype in: 228 in: int in: Tworkfile in: 5 in: 423 in: 423 in: Tworkfile in: 4 in: 3191 in: 3191 in: Ttype in: 329 in: X/'/ w out: Hclean out: 0 out: Hcurrent out: 5 out: Hgrowdata out: 336 out: 47 out: 47 out: bz_census.c: #2799 ?can't create "bz_census.h" out: Hcheck0 out: 1 out: Hack in: Tcheck out: Hcheck out: 1 in: Tack out: Hsetdot out: 383 out: 383 out: Hsetpat out: 1 out: ' out: Hunlock in: Torigin in: 2 out: Horigin out: 355 in: Ttype in: 383 in: !co -l in: Tdclick out: Hsetdot out: 370 out: 379 out: Hunlockfile out: 1 in: Tdclick out: Hsetdot out: 370 out: 381 out: Hunlockfile out: 1 in: Tworkfile in: 5 in: 423 in: 423 in: Tsend out: Hsnarflen out: Hgrowdata out: 390 out: 12 out: 12 out: bz_census.h out: Hcheck0 out: 1 out: Hack in: Tcheck out: Hcheck out: 1 in: Tack out: Hsetdot out: 402 out: 402 out: Hgrow out: 402 out: 65 out: Hcheck0 out: 1 out: Hack in: Trequest in: 402 in: 65 out: Hdata out: 1 out: 402 out: 65 out: RCS/bz_census.h,v --> bz_census.h revision 1.1 (locked) done ! in: Tcheck out: Hcheck out: 1 in: Tack out: Hsetdot out: 467 out: 467 out: Hunlock in: Tworkfile in: 5 in: 423 in: 423 in: Ttype in: 467 in: X w