Extract intergenic regions from processed GenBank annotations

Share:

Description

Extract the intergenic regions from a set of GenBank annotations.

Usage

1
2
## S4 method for signature 'GenBankRecord'
intergenic(x)

Arguments

x

A GenBankRecord object

Value

A GRanges for the intergenic regions, defined as regions not overlapping any genes defined in the annotations on either strand.

Examples

1
2
gb = readGenBank(system.file("sample.gbk", package="genbankr"))
intergenic(gb)