genomation
genomation copied to clipboard
output annotation meta columns
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