## ----setup, include = FALSE-------------------------------------------------------------------------------------------
options(width = 120)
knitr::opts_chunk$set(
   collapse = TRUE,
   eval = interactive(),
   echo = TRUE,
   comment = "#>"
)

## ----eval=FALSE-------------------------------------------------------------------------------------------------------
# library(igvR)
# igv <- igvR()
# setBrowserWindowTitle(igv, "hg38 explicit")
# setCustomGenome(igv,
#                 id = "hg38",
#                 genomeName = "Human (GRCh38/hg38)",
#                 fastaURL = "https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg38/hg38.fa",
#                 fastaIndexURL = "https://s3.amazonaws.com/igv.broadinstitute.org/genomes/seq/hg38/hg38.fa.fai",
#                 cytobandURL = "https://s3.amazonaws.com/igv.broadinstitute.org/annotations/hg38/cytoBandIdeo.txt",
#                 chromosomeAliasURL = NA,
#                 geneAnnotationName = "Refseq Genes",
#                 geneAnnotationURL = "https://s3.amazonaws.com/igv.org.genomes/hg38/refGene.txt.gz",
#                 geneAnnotationTrackHeight = 500,
#                 geneAnnotationTrackColor = "darkBlue",
#                 initialLocus = "chr5:88,621,308-89,001,037",
#                 visibilityWindow = 5000000)

## ----eval=FALSE-------------------------------------------------------------------------------------------------------
# base.url <- "https://gladki.pl/igvR/testFiles/sarsGenome"
# fasta.file <- sprintf("%s/%s", base.url, "Sars_cov_2.ASM985889v3.dna.toplevel.fa")
# fastaIndex.file <- sprintf("%s/%s", base.url, "Sars_cov_2.ASM985889v3.dna.toplevel.fa.fai")
# annotation.file <- sprintf("%s/%s", base.url, "Sars_cov_2.ASM985889v3.101.gff3")
# 
# Sys.sleep(2)
# 
# setCustomGenome(igv,
#                 id = "Sars_cov_2",
#                 genomeName = "Sars_cov_2.ASM985889v3",
#                 fastaURL = fasta.file,
#                 fastaIndexURL = fastaIndex.file,
#                 geneAnnotationURL = annotation.file,
#                 geneAnnotationName = "ASM985889v3",
#                 geneAnnotationTrackHeight = 500,
#                 geneAnnotationTrackColor = "darkBlue",
#                 visibilityWindow = 30000)

## ----eval=TRUE--------------------------------------------------------------------------------------------------------
sessionInfo()