Skip to contents

Check genome build

Source: R/check_genome_build.R
check_genome_build.Rd

Check if the genome build is valid and return the appropriate BSGenome object.

Usage

check_genome_build(genome_build)

Arguments

genome_build

A character string with the abbreviated genome build name, or a BSGenome object. At the moment, only hg38 and hg19 are supported as abbreviated input.

Value

A BSGenome object.

See also

BSgenome-class for more information on BSGenome objects.

Examples

if (requireNamespace("BSgenome.Hsapiens.UCSC.hg38", quietly = TRUE)) {
    check_genome_build("hg38")
}
#> | BSgenome object for Human
#> | - organism: Homo sapiens
#> | - provider: UCSC
#> | - genome: hg38
#> | - release date: 2023/01/31
#> | - 711 sequence(s):
#> |     chr1                    chr2                    chr3                   
#> |     chr4                    chr5                    chr6                   
#> |     chr7                    chr8                    chr9                   
#> |     chr10                   chr11                   chr12                  
#> |     chr13                   chr14                   chr15                  
#> |     ...                     ...                     ...                    
#> |     chr19_KV575256v1_alt    chr19_KV575257v1_alt    chr19_KV575258v1_alt   
#> |     chr19_KV575259v1_alt    chr19_KV575260v1_alt    chr19_MU273387v1_alt   
#> |     chr22_KN196485v1_alt    chr22_KN196486v1_alt    chr22_KQ458387v1_alt   
#> |     chr22_KQ458388v1_alt    chr22_KQ759761v1_alt    chrX_KV766199v1_alt    
#> |     chrX_MU273395v1_alt     chrX_MU273396v1_alt     chrX_MU273397v1_alt    
#> | 
#> | Tips: call 'seqnames()' on the object to get all the sequence names, call
#> | 'seqinfo()' to get the full sequence info, use the '$' or '[[' operator to
#> | access a given sequence, see '?BSgenome' for more information.