update trio-binning

This commit is contained in:
chhylp123
2020-04-10 23:32:15 -04:00
parent ad1b79a6bb
commit 20c6c1d2e6
3 changed files with 17 additions and 31 deletions
+14 -13
View File
@@ -13467,9 +13467,10 @@ kvec_asg_arc_t_warp* new_rtg_edges)
delete_useless_nodes(ug);
deduplicate(*ug, read_g, reverse_sources, 20, 100, 0.05, 0.2, ruIndex, 1);
delete_useless_nodes(ug);
// purge_dups(*ug, read_g, coverage_cut, reverse_sources, ruIndex, new_rtg_edges, 0.75, 50, 50, 0.5, max_hang,
// min_ovlp, bubble_dist, drop_ratio, 1);
// deduplicate(*ug, read_g, reverse_sources, 20, 100, 0.05, 0.2, ruIndex, 1);
// delete_useless_nodes(ug);
update_unitig_graph((*ug), read_g, reverse_sources, ruIndex, flag, drop_rate);
@@ -13485,17 +13486,17 @@ kvec_asg_arc_t_warp* new_rtg_edges)
renew_utg(ug, read_g, new_rtg_edges);
update_unitig_graph((*ug), read_g, reverse_sources, ruIndex, flag, drop_rate);
renew_utg(ug, read_g, new_rtg_edges);
delete_useless_nodes(ug);
purge_dups(*ug, read_g, coverage_cut, reverse_sources, ruIndex, new_rtg_edges, 0.75, 50, 50, 0.5, max_hang,
min_ovlp, bubble_dist, drop_ratio, 1);
delete_useless_nodes(ug);
set_drop_trio_flag(*ug);
n_vtx = read_g->n_seq;
for (v = 0; v < n_vtx; v++)
{
if(read_g->seq[v].c == HAP_LABLE)
{
read_g->seq[v].c = PRIMARY_LABLE;
}
}
/**
kv_destroy(new_rtg_nodes.a);