inst/extdata/GentlemanLab/BSgenome.Ptroglodytes.UCSC.panTro5-tools/sort_panTro5.2bit.R

library(Biostrings)
library(rtracklayer)
panTro5 <- import("panTro5.2bit")

seqlevels1 <- paste0("chr", c(1, "2A", "2B", 3:22, "X", "Y", "M"))

tmp <- CharacterList(strsplit(names(panTro5), "_"))
npart <- elementNROWS(tmp)
stopifnot(all(npart <= 4L))

idx1 <- which(npart == 1L)
stopifnot(setequal(names(panTro5)[idx1], seqlevels1))

idx2 <- which(npart == 2L)
stopifnot(length(idx2) == 0L)

idx4 <- which(npart == 4L)
m4 <- matrix(unlist(tmp[idx4]), ncol=4L, byrow=TRUE)
stopifnot(all(m4[ , 2L] == "NW"))
stopifnot(all(m4[ , 4L] == "random"))
chrom <- match(m4[ , 1L], seqlevels1)
stopifnot(!anyNA(chrom))
oo4 <- order(chrom, m4[ , 3L])
seqlevels4 <- names(panTro5)[idx4[oo4]]

idx3 <- which(npart == 3L)
m3 <- matrix(unlist(tmp[idx3]), ncol=3L, byrow=TRUE)
stopifnot(all(m3[ , 1L] == "chrUn"))
stopifnot(all(m3[ , 2L] == "NW"))
oo3 <- order(m3[ , 3L])
seqlevels3 <- names(panTro5)[idx3[oo3]]

panTro5_seqlevels <- c(seqlevels1, seqlevels4, seqlevels3)

export(panTro5[panTro5_seqlevels], "panTro5.sorted.2bit")

Try the BSgenome package in your browser

Any scripts or data that you put into this service are public.

BSgenome documentation built on Nov. 8, 2020, 7:48 p.m.