Skip to content

Commit

Permalink
Extend transcript boundaries instead of skipping the transcript
Browse files Browse the repository at this point in the history
This modifies the commit d6cf97b, it is better to include
such transcript and modify their boundaries instead of omitting
the transcripts entirely

Resolves #2323
  • Loading branch information
pd3 committed Dec 11, 2024
1 parent d6cf97b commit e686d37
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 3 deletions.
3 changes: 2 additions & 1 deletion NEWS
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,8 @@ Changes affecting specific commands:

* bcftools csq

- Check the input GFF for features outside transcript boundaries (#2323)
- Check the input GFF for features outside transcript boundaries and extend the transcript
to contain the feature fully (#2323)

* bcftools query

Expand Down
4 changes: 2 additions & 2 deletions gff.c
Original file line number Diff line number Diff line change
Expand Up @@ -1017,9 +1017,9 @@ int gff_parse(gff_t *gff)
if ( !gff->warned.ftr_out_of_bounds || gff->verbosity > 1 )
fprintf(stderr,"Warning: The GFF contains features outside the transcript boundaries .. %s\n",gff_id2string(gff,transcript,tr->id));
gff->warned.ftr_out_of_bounds++;
continue;
if ( ftr->beg < tr->beg ) tr->beg = ftr->beg;
if ( ftr->end > tr->end ) tr->end = ftr->end;
}

tr->used = 1;
tr->gene->used = 1;

Expand Down

0 comments on commit e686d37

Please sign in to comment.