Commit c2be46d8 authored by cziegenhain's avatar cziegenhain
Browse files

barcode sharing fix

parent aaf6b5f3
Loading
Loading
Loading
Loading
+2 −0
Original line number Diff line number Diff line
@@ -31,6 +31,8 @@ We provide a script to convert zUMIs output into loom file automatically based o
zUMIs will try to automatically do this, otherwise convert zUMIs output to loom by simply running `Rscript rds2loom.R myRun.yaml`.

## Changelog
07 Apr 2021: zUMIs2.9.6: Fixed a potential crash when the input file to the barcode sharing feature does not contain an equal number of columns.

18 Feb 2021: zUMIs2.9.5: Mismatches for detecting Smart-seq3 UMI-read pattern are now user-settable in the YAML file as follows: `find_pattern: ATTGCGCAATG;2` would allow 2 mismatches. Usage as prior to this version (ie `find_pattern: ATTGCGCAATG`) will default to the previous value of 1 mismatch allowed.

18 Sept 2020 - 29 Nov Sept 2020: zUMIs.2.9.4b/c/d/e/f: Fix & speed up Smart-seq3 UMI read counting. Prevent crash when a chunk of cell BCs does not match any downsampling. Speed up barcode detection steps for some cases. Prevent too much CPU usage in UMI error correction. Take correct samtools executable in gene annotation parsing. Prevent crash in BC error correction for huge datasets.
+2 −1
Original line number Diff line number Diff line
@@ -298,11 +298,12 @@ BCbin <- function(bccount_file, bc_detected) {
  }

  if(!is.null(opt$barcodes$barcode_sharing)){
    share_table <- data.table::fread( opt$barcodes$barcode_sharing, header = F, skip = 1)
    share_table <- data.table::fread( opt$barcodes$barcode_sharing, header = F, skip = 1, fill = TRUE)
    if(ncol(share_table) > 2){ #flatten table more if necessary
      share_table <- data.table::melt(share_table, id.vars = "V1")[,variable := NULL]
    }
    setnames(share_table, c("main_bc","shared_bc"))
    share_table <- share_table[shared_bc != '']

    share_mode <- data.table::fread( opt$barcodes$barcode_sharing, header = F, nrows = 1)$V1
    share_mode <- as.numeric(unlist(strsplit(gsub(pattern = "#",replacement = "", x = share_mode),"-")))
+1 −1
Original line number Diff line number Diff line
@@ -3,7 +3,7 @@
# Pipeline to run UMI-seq analysis from fastq to read count tables.
# Authors: Swati Parekh, Christoph Ziegenhain, Beate Vieth & Ines Hellmann
# Contact: sparekh@age.mpg.de or christoph.ziegenhain@ki.se
vers=2.9.5
vers=2.9.6
currentv=$(curl -s https://raw.githubusercontent.com/sdparekh/zUMIs/main/zUMIs.sh | grep '^vers=' | cut -f2 -d "=")
if [ "$currentv" != "$vers" ] ; then
    echo -e "------------- \n\n Good news! A newer version of zUMIs is available at https://github.com/sdparekh/zUMIs \n\n-------------";