Skip to content

[Bug] Promoter intervals exceed chromosome boundaries on chrM (hg38) #73

@XinzeWu

Description

@XinzeWu

version

GenomicFeatures_1.58.0

Description

Generated promoter intervals include coordinates outside chromosome boundaries. For example, on hg38 chrM (length = 16,569 bp):

chrM   -924   1076
chrM   -853   1147
chrM  15523  17523

The first two have negative starts, and the last exceeds the chromosome length.

Expected

All BED intervals should satisfy:
0 <= start < end <= chromSize.

Suggestion

When extending TSS to promoters, clip intervals to chromosome sizes (e.g. with bedtools slop -g hg38.sizes or GenomicRanges::trim() in R). Add validation or warnings for out-of-bound coordinates.


Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions