read.genome: Read the genome of a given organism

Description Usage Arguments Details Value Author(s) Examples

Description

Read an organism specific genome stored in fasta/fastq file format.

Usage

1

Arguments

file

a character string specifying the path to the file storing the genome.

format

a character string specifying the file format used to store the genome, e.g. format = "fasta" or format = "fastq".

...

additional arguments that are used by the readDNAStringSet function.

Details

The read.genome function takes a string specifying the path to the genome file of interest as first argument.

Genomes stored in fasta files can be downloaded from http://ensemblgenomes.org/info/genomes.

Value

A data.frame storing the gene id in the first column, the corresponding sequence as string in the second column, and the sequence length in the third column.

Author(s)

Hajk-Georg Drost

Examples

1
2
3
4
5
# reading a genome stored in a fasta file
Ath.genome <- read.genome(system.file('seqs/ortho_thal_cds.fasta', package = 'seqreadr'),
                           format = "fasta")

dplyr::glimpse(Ath.genome)

HajkD/seqreadr documentation built on May 6, 2019, 10:55 p.m.