buildCodonDB-methods: build a codon sqlLite DB

buildCodonDBR Documentation

build a codon sqlLite DB

Description

Index genome according to codon positions in coding sequence.
This can be useful when downstream analysis or output of alignements should only be applied to certain codons.
For example calculating genetic distance on only 4-fold sites or outputting fasta files based on codon position.

Usage

buildCodonDB(genome, exons, sqlDir = NULL, ...)

Arguments

genome

character or DNAStringSet
Import fasta reference genome from file by passsing a character or provide a preloaded (using Biostrings::readDNAStringSet()) DNAStringSet

exons

GrangesList
Generated using getFeatures() by passing "gene:cds" to the feature argument.

sqlDir

character set to sql output to run in lower memory mode.

position

character Options are "all", default, or any combination of c("first", "second", "third").
Output will contain only specified codon positions.

fourFoldCodon

character Options are NULL, "only" or "remove".
Determines what to do with 4-fold degenerate sites.

Value

An object of class GRanges containing codon ranges for positions specified.


CMWbio/geaR documentation built on April 22, 2023, 6:23 a.m.