Skip to contents

The DNA sibling of tales_get_protein_seq(): reassembles each array's parts, in order, into one full-length nucleotide sequence, by pasting dna_seq together directly. See tales_get_protein_seq()'s Details for why position_in_array is the ordering key even for a tales_msa.

Usage

tales_get_dna_seq(x)

Arguments

x

A tales object carrying a dna_seq column.

Value

A Biostrings::DNAStringSet, named by array_id.

Examples

x <- tales_from_telltales(system.file("extdata", "tellTaleExampleOutput",
                                      package = "tantale"))
tales_get_dna_seq(x)[1]
#> DNAStringSet object of length 1:
#>     width seq                                               names               
#> [1]  4305 ATGGATCCCATTCGTCCGCGCAG...TGAGGGAGCTATTGCCTCAGTGA ROI_00001