Reduce 3' UTRs to the gene level across all transcripts — .reduce_features • SeedMatchR

Reduce 3' UTRs to the gene level across all transcripts

Usage

.reduce_features(gr.list, txdb, filter)

Arguments

gr.list: Input GenomicRanges::GRangesList file with features grouped by transcript
txdb: An ensembldb::EnsDb object
filter: An AnnotationFilter::AnnotationFilterList for selecting transcripts

Value

A gtf of unique feature intervals per gene

Examples

if (FALSE) { # interactive()
anno.db = SeedMatchR::load_species_anno_db("hg38")

features = SeedMatchR::get_features_gtf(anno.db$txdb)

gtf = SeedMatchR::.reduce_features(features, anno.db$txdb)
}