Skip to content

Commit

Permalink
#31 added output message and correct use of unzipped annotation file …
Browse files Browse the repository at this point in the history
…| fixed some minor log messages
  • Loading branch information
OliverVoogd committed Aug 8, 2024
1 parent ea9d870 commit ad3a218
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 3 deletions.
2 changes: 2 additions & 0 deletions R/find_isoform.R
Original file line number Diff line number Diff line change
Expand Up @@ -54,8 +54,10 @@ find_isoform_bambu <- function(annotation, genome_fa, genome_bam, outdir, config
useTempAnnot <- FALSE
bambuTempAnnot <- ""
if (stringr::str_ends(annotation, ".gz")) {
cat("Unzipping annotation file for bambu\n")
useTempAnnot <- TRUE
bambuTempAnnot <- R.utils::gunzip(annotation, remove=FALSE)
annotation <- bambuTempAnnot # override using zipped annotation file
}
bambuAnnotations <- bambu::prepareAnnotations(annotation)
# Tmp fix: remove withr if bambu imports seqlengths properly
Expand Down
6 changes: 3 additions & 3 deletions R/sc_long_pipeline.R
Original file line number Diff line number Diff line change
Expand Up @@ -154,8 +154,8 @@ sc_long_pipeline <-
# run flexiplex
cat(format(Sys.time(), "%X %a %b %d %Y"), "Demultiplexing using flexiplex...\n")
if (!file.exists(barcodes_file)) {
stop("The brcodes_file ",barcodes_file, " doesn't exists. Please check the path.",
" If you do not have a barcode file, please set `brcodes_file = NULL` to run BLAZE.")
stop("The barcodes_file ",barcodes_file, " doesn't exists. Please check the path.",
" If you do not have a barcode file, please set `barcodes_file = NULL` to run BLAZE.")
}
cat("Matching cell barcodes...\n")
infq <- file.path(outdir, "matched_reads.fastq")
Expand All @@ -179,7 +179,7 @@ sc_long_pipeline <-
} else {
infq <- fastq
cat("Skipping Demultiplexing step...\n")
cat("Please make sure the `",infq,"`` is the the demultiplexing output from previous FLAEMS call.\n")
cat("Please make sure the `",infq,"`` is the demultiplexing output from previous FLAMES call.\n")
} # requesting to not match barcodes implies `fastq` has already been run through the
# function in a previous FLAMES call
cat("Running FLAMES pipeline...\n")
Expand Down

0 comments on commit ad3a218

Please sign in to comment.