genomation icon indicating copy to clipboard operation
genomation copied to clipboard

output annotation meta columns

Open shawpa opened this issue 3 years ago • 0 comments

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

shawpa avatar May 18 '21 19:05 shawpa