/usr/lib/R/site-library/rtracklayer/scripts/twobit.R is in r-bioc-rtracklayer 1.38.0-1build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | ## code to get twobit support working
library(rtracklayer)
library(Biostrings)
system.time(chr11 <-
read.DNAStringSet("~/genomes/HuRef/seqs/hs_alt_HuRef_chr11.fa"))
export(chr11, "chr11.2bit")
seqinfo(TwoBitFile("chr11.2bit"))
system.time(chr11_2bit <- import("chr11.2bit"))
identical(chr11_2bit, chr11)
export(subseq(chr11, 1, 100), "test.2bit")
|