Hi, I have a
GRanges object e.g.
GRanges object with 6 ranges and 1 metadata column: seqnames ranges strand | gene_id <Rle> <IRanges> <Rle> | <character> 15S_rRNA MT [ 6546, 8194] + | 15S_rRNA 21S_rRNA MT [58009, 62447] + | 21S_rRNA Q0010 MT [ 3952, 4338] + | Q0010 Q0017 MT [ 4254, 4415] + | Q0017 Q0032 MT [11667, 11957] + | Q0032 Q0045 MT [13818, 26701] + | Q0045 ------- seqinfo: 17 sequences (1 circular) from an unspecified genome; no seqlengths
Now I would like to extract from it only specific rows. Lets say I want to have the genes
I have tried to do it iwht the subset command, but it only takes one pattern
subset(genes.mod, gene_id =="Q0017
Is there a way to add multiple patterns w.o. doing it multiple times for each gene?