Skip to content

Commit

Permalink
v1.0.4
Browse files Browse the repository at this point in the history
  • Loading branch information
Yan Gao committed Oct 6, 2020
1 parent 681cee4 commit c92f8c5
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions sub_example.c
Original file line number Diff line number Diff line change
Expand Up @@ -104,12 +104,12 @@ int main(void) {
// perform abpoa-msa
abpoa_res_t res;
for (i = 0; i < n_seqs; ++i) {
res.graph_cigar = 0, res.n_cigar = 0;
res.graph_cigar = 0, res.n_cigar = 0; res.is_rc = 0;
abpoa_align_sequence_to_subgraph(ab, abpt, beg_end_id[i][0], beg_end_id[i][1], bseqs[i], seq_lens[i], &res);
int exc_beg, exc_end;
if (i != 0) abpoa_subgraph_nodes(ab, beg_end_id[i][0], beg_end_id[i][1], &exc_beg, &exc_end);
else exc_beg = 0, exc_end = 0;
abpoa_add_subgraph_alignment(ab, abpt, exc_beg, exc_end, bseqs[i], seq_lens[i], res.n_cigar, res.graph_cigar, i, n_seqs);
abpoa_add_subgraph_alignment(ab, abpt, exc_beg, exc_end, bseqs[i], seq_lens[i], res, i, n_seqs);
if (res.n_cigar) free(res.graph_cigar);
}
if (abpt->out_cons) {
Expand All @@ -118,7 +118,7 @@ int main(void) {
fprintf(stderr, "Warning: no consensus sequence generated.\n");
}
if (abpt->out_msa) {
abpoa_generate_rc_msa(ab, abpt, NULL, n_seqs, stdout, NULL, NULL);
abpoa_generate_rc_msa(ab, abpt, NULL, NULL, n_seqs, stdout, NULL, NULL);
}
// abpoa_reset_graph(ab, abpt, seq_lens[0]); // reset graph before re-use

Expand Down

0 comments on commit c92f8c5

Please sign in to comment.