Skip to content
Permalink
Browse files

load seqinfo from BSgenome.Hsapiens

  • Loading branch information
gavinha committed Apr 18, 2019
1 parent df3d74d commit 2f660f374ebdf04765285b51b92e03298b0e159f
Showing with 6 additions and 1 deletion.
  1. +6 −1 code/titanCNA_v1.21.0_TenX.R
@@ -144,7 +144,12 @@ if (is.null(outplot)){
outImage <- gsub(".titan.txt", ".RData", outfile)

## set up chromosome naming convention ##
seqinfo <- Seqinfo(genome=genomeBuild)
bsg <- paste0("BSgenome.Hsapiens.UCSC.", genomeBuild)
if (!require(bsg, character.only=TRUE, quietly=TRUE, warn.conflicts=FALSE)) {
seqinfo <- Seqinfo(genome=genomeBuild)
} else {
seqinfo <- seqinfo(get(bsg))
}
seqlevelsStyle(seqinfo) <- genomeStyle
seqlevelsStyle(chrs) <- genomeStyle

0 comments on commit 2f660f3

Please sign in to comment.
You can’t perform that action at this time.