We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 181d2c5 commit 7e77044Copy full SHA for 7e77044
sc2ts/utils.py
@@ -1320,7 +1320,7 @@ def get_recombinant_samples(ts):
1320
# Search the subtree for a causal sample.
1321
for v in tree.nodes(u, order="levelorder"):
1322
child = ts.node(v)
1323
- if child.is_sample() and child.metadata["date"] == recomb_date:
+ if child.is_sample() and child.metadata["date"] <= recomb_date:
1324
edge = ts.edge(tree.edge(v))
1325
assert edge.left == 0 and edge.right == ts.sequence_length
1326
causal_sample = child
0 commit comments