Nothing
### =========================================================================
### Tabix support
### -------------------------------------------------------------------------
###
setMethod("import", "TabixFile",
function(con, format, text,
which = if (is.na(genome)) NULL
else as(seqinfoForGenome(genome), "GenomicRanges"),
genome = NA, header = TRUE, ...)
{
if (missing(format)) {
format <- file_ext(file_path_sans_ext(path(con)))
}
m <- manager()
buffer <- queryForResource(m, con, which, header = header)
on.exit(release(m, buffer))
file <- try(FileForFormat(buffer, format), silent = TRUE)
if (is(file, "try-error")) {
tabixHeader <- headerTabix(con)
args <- list(...)
if (is.null(args$comment.char))
args$comment.char <- tabixHeader$comment
skip <- tabixHeader$skip
if (header) {
skip <- skip - 1L
}
do.call(import.tabSeparated,
c(list(buffer), genome = genome, tabixHeader$indexColumns,
skip = skip, header = header, args))
} else import(file, genome = genome, ...)
})
setGeneric("exportToTabix",
function(object, con, ...) standardGeneric("exportToTabix"))
setMethod("exportToTabix", c("ANY", "character"),
function(object, con, ...) {
con <- TabSeparatedFile(con)
export(sort(object, ignore.strand=TRUE), con,
row.names=FALSE, col.names=TRUE, ...)
indexTrack(con, seq=1L, start=2L, end=3L, skip=1L)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.