static void sw_set_default_scoring() { scoring_system_default(&scoring); // Change slightly scoring.match = 2; scoring.mismatch = -2; scoring.gap_open = -2; scoring.gap_extend = -1; }
CallDecomp* call_decomp_init(htsFile *vcffh, bcf_hdr_t *vcfhdr) { CallDecomp *dc = ctx_calloc(1, sizeof(CallDecomp)); dc->nw_aligner = needleman_wunsch_new(); dc->aln = alignment_create(1024); dc->scoring = ctx_calloc(1, sizeof(dc->scoring[0])); scoring_system_default(dc->scoring); dc->vcffh = vcffh; dc->vcfhdr = vcfhdr; dc->v = bcf_init(); strbuf_alloc(&dc->sbuf, 256); return dc; }