addGeneAligns {spliceSites} | R Documentation

## Reads a bamRange object for a given `bamReader`, `refGenome` and gene name.

### Description

Locates gene in genome via `refGenome` and reads a `bamRange` from the determined region.

### Usage

        addGeneAligns(x)

### Arguments

`x` | gapSites. The result contains a copy of the passed object.

### Details

The function adds a `gene_aligns` column to the contained `data.frame`.

### Value

gapSites

### Author(s)

Wolfgang Kaisers

### Examples

    # A) Read gapSites
    bam <- system.file("extdata", "rna_fem.bam", package="spliceSites")
    reader <- bamReader(bam, idx=TRUE)
    ga <- alignGapList(reader)
    bamClose(reader)

    # B) Annotate
    ucf <- system.file("extdata", "uc_small_junc.RData", package="spliceSites")
    ucj <- loadGenome(ucf)
    annotation(ga) <- annotate(ga, ucj)

    # C) align part
    gal <- addGeneAligns(ga)
    gal

[Package spliceSites version 1.28.0 Index]