sdparekh / zUMIs

zUMIs: A fast and flexible pipeline to process RNA sequencing data with UMIs
GNU General Public License v3.0
275 stars 67 forks source link

BUG Samtools related bug inDrops v3 samples #240

Closed droplet-lab closed 3 years ago

droplet-lab commented 3 years ago

Describe the bug I get a samtools error when running one of my samples, but it is fine for the other samples, any ideas where this could come from ?

I have tried a few things but nothing has been successful so far. Any help would greatly be appreciated :). Thanks a lot

Warning message: In data.table::fread(cmd = paste(samtools, "view", filtered_bams[1], : File '/tmp/RtmpQBqt4O/file59e963125f9c' has size 0. Returning a NULL data.table. Usage: samtools cat [options] [... ] samtools cat [options] [... ]

Concatenate BAM or CRAM files, first those in , then those on the command line.

Options: -b FILE list of input BAM/CRAM file names, one per line -h FILE copy the header from FILE [default is 1st input file] -o FILE output BAM/CRAM Usage: samtools cat [options] [... ] samtools cat [options] [... ]

Concatenate BAM or CRAM files, first those in , then those on the command line.

Options: -b FILE list of input BAM/CRAM file names, one per line -h FILE copy the header from FILE [default is 1st input file] -o FILE output BAM/CRAM rm: cannot remove '/PATH/tmp.SAMPLE.*': No such file or directory Warning message: In system(paste(STAR_command, "&", sammerge_command, "& wait")) : error in running command Error in gsub("SN:", "", chr) : object 'chr' not found Calls: .makeSAF ... .chromLengthFilter -> [ -> [.data.table -> eval -> eval -> gsub In addition: Warning message: In data.table::fread(bread, col.names = c("chr", "len"), header = F) : File '/tmp/RtmpcLgHAv/file5a4e6ca5199c' has size 0. Returning a NULL data.table. Execution halted Loading required package: yaml Loading required package: Matrix Error in gzfile(file, "rb") : cannot open the connection Calls: rds_to_loom -> readRDS -> gzfile In addition: Warning message: In gzfile(file, "rb") : cannot open compressed file '/PATH/Sample.dgecounts.rds', probable reason 'No such file or directory' Execution halted [E::hts_open_format] Failed to open file NA samtools view: failed to open "NA" for reading: No such file or directory Error in gzfile(file, "rb") : cannot open the connection Calls: readRDS -> gzfile In addition: Warning message: In gzfile(file, "rb") : cannot open compressed file '/PATH/Sample.dgecounts.rds', probable reason 'No such file or directory' Execution halted

To Reproduce project:SAMPLE sequence_files: file1: name: /PATH/SAMPLE_R1.fastq.gz base_definition: cDNA(1-61) file2: name: /PATH/SAMPLE_R2.fastq.gz base_definition: BC(1-8) file3: name: /PATH/SAMPLE_R3.fastq.gz base_definition: BC(1-8) file4: name: /PATH/SAMPLE_R4.fastq.gz base_definition:

Screenshots If applicable, add screenshots to help explain your problem.

Desktop (please complete the following information):

Additional context Add any other context about the problem here.

droplet-lab commented 3 years ago

Looks like it was another issue with my environment, re-installing everything seems to run.

Sorry about this