File tree 1 file changed +8
-2
lines changed
1 file changed +8
-2
lines changed Original file line number Diff line number Diff line change 12
12
hificoveragefile = sys .argv [8 ]
13
13
min_hifi_solid_coverage = int (sys .argv [9 ])
14
14
15
+
16
+ length_solid_node_threshold = 200000
17
+
15
18
def find (s , parent ):
16
19
while parent [s ] != parent [parent [s ]]:
17
20
parent [s ] = parent [parent [s ]]
@@ -56,19 +59,22 @@ def canon(left, right):
56
59
if parts [0 ] in unique_nodes : continue
57
60
hifi_solid_nodes .add (parts [0 ])
58
61
62
+ length_solid_nodes = set ()
59
63
ont_solid_nodes = set ()
60
64
with open (ontcoveragefile ) as f :
61
65
for l in f :
62
66
parts = l .strip ().split ('\t ' )
63
67
if parts [2 ] == "coverage" : continue
64
- if float (parts [2 ]) < min_ont_solid_coverage : continue
65
68
if parts [0 ] in unique_nodes : continue
66
- ont_solid_nodes .add (parts [0 ])
69
+ if float (parts [2 ]) >= min_ont_solid_coverage : ont_solid_nodes .add (parts [0 ])
70
+ if int (parts [1 ]) >= length_solid_node_threshold : length_solid_nodes .add (parts [0 ])
67
71
68
72
solid_nodes = set ()
69
73
for node in hifi_solid_nodes :
70
74
if node in ont_solid_nodes :
71
75
solid_nodes .add (node )
76
+ for node in length_solid_nodes :
77
+ solid_nodes .add (node )
72
78
73
79
solid_edges = set ()
74
80
for key in edge_coverage :
You can’t perform that action at this time.
0 commit comments