Skip to content

output annotation meta columns #201

@shawpa

Description

@shawpa

I'm wondering if there is a way to link the meta data to the annotation result for non-gene part annotation. For example, if you do the following, you can get the feature name linked to the intersection data.

head(data.frame([email protected]))
target.row dist.to.feature feature.name feature.strand
12504 1 226 NR_149138 +
1250.1102 2 234 NR_149138 +
1250.2102 3 236 NR_149138 +
1250.3102 4 248 NR_149138 +
12518 5 -251 NR_110795 +
1251.1102 6 -249 NR_110795 +

If I have annotation from the ensembl regulatory build, I'd also like to know not only if my site overlaps with a feature but also the name of the feature it is overlapping. I know how to read in the information for the regulatory build and do the intersection, but can I link it to the ensembl regulatory feature name instead of it just showing 0 or 1. I've looked at the slots available in the class and there doesn't seem to be a way to do it.

       seqnames              ranges strand |            name    category
          <Rle>           <IRanges>  <Rle> |     <character> <character>
   [1]     chrX   10001201-10001601      * | ENSR00001290757    enhancer
   [2]     chrX 100081001-100081200      * | ENSR00000911315    enhancer
   [3]     chrX   10029001-10029200      * | ENSR00001290758    enhancer
   [4]     chrX 100382401-100383000      * | ENSR00000911348    enhancer
   [5]     chrX 100438200-100438601      * | ENSR00001161921    enhancer

Thanks,

Annie

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions