Skip to content

Commit c248876

Browse files
committed
gene_caller: Avoid using primary transcript in doc string, since the anchor transcript does not define which one is primary
1 parent 2ededb2 commit c248876

1 file changed

Lines changed: 3 additions & 3 deletions

File tree

src/geneml/gene_caller.py

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -596,10 +596,10 @@ def select_gene_calls_per_group(group: list[tuple[float, list[GeneEvent]]], max_
596596

597597
@njit
598598
def split_into_genes(group: list[tuple[float, list[GeneEvent]]]) -> list:
599-
"""Split overlapping gene calls into gene groups anchored by primary transcripts.
599+
"""Split overlapping gene calls into gene groups.
600600
601-
Uses primary transcripts (selected by start position and score) as anchors to define
602-
gene loci boundaries. Candidates starting before an anchor's end are grouped together,
601+
Uses anchor transcripts (selected by start position and score) to define gene loci boundaries.
602+
Candidates starting before an anchor's end are grouped together,
603603
while candidates starting after initiate a new group. Transcripts overlapping multiple
604604
anchors are assigned to the first group.
605605

0 commit comments

Comments
 (0)