slamdunk all: No tcount files and no error

[andtsouris]

@t-neumann @jkobject While running slamdunk all of the jkobject/slamdunk fork on a set of 10 samples, slamdunk doesn't produce tcount files for 2 of them without any error or other indication of what went wrong.
The same behaviour was reproduced with both the forked and original versions, when running slamdunk count on the filtered bam file.
Do you know what went wrong and how to fix it?

[jkobject]

Hi, if this error happens even on the un-forked version, I am not sure it concerns the modifications I made..?

Do you have code and/or input data we could look at to understand or reproduce your error?

[andtsouris]

Thanks for the quick response! I too think that it's not due to the modifications you made.
My call of slamdunk is pretty basic slamdunk all -r genome.fa -b utr3.bed -o out/${sample} -t 12 ${sample}_R1_001.fastq.gz ${sample}_R2_001.fastq.gz.

Other than that, I'm only importing VarScan prior to calling slamdunk with alias varscan="java -jar /home/myuser/Tools/VarScan.v2.4.4.jar".

What is weird is that all samples are processed in exactly the same way prior to slamdunk but only some fail.

[t-neumann]

Hm that is indeed weird if it works for some but not for all. I second @jkobject I would probably need one of the filtered bam files to see what's going on there