Entering edit mode
3.1 years ago
gonzalezb549
•
0
Hi very new to bioinformatics and am doing an analysis of single-cell sequencing fastq data from a previously published paper. I downloaded the FASTQ files and put them in their own directory. Then wrote this script based on an example given to me. The error I get is exit code 127 cellranger: command not found.
#!/bin/bash
#SBATCH --job-name=cellranger
#SBATCH --mail-type=END,FAIL
#SBATCH --ntasks=1
#SBATCH --mem=40gb
#SBATCH --cpus-per-task=4
#SBATCH --time=24:00:00
#SBATCH --output=cellrangerProject3.log
source ~/.bashrc
export PATH=/network/rit/misc/scratch/single_cell/cellranger/cellranger=3.0.2:$PATH
cellranger count --id=noinjury1 --transcriptome=/network/rit/misc/scratch/single_cell/refdata-cellranger-mm10-3.0.0 --fastqs=/network/rit/lab/bioinformaticslab/BGonzalez/Project3/Noinjury --sample=noinjury1
cellranger count --id=noinjury2 --transcriptome=/network/rit/misc/scratch/single_cell/refdata-cellranger-mm10-3.0.0 --fastqs=/network/rit/lab/bioinformaticslab/BGonzalez/Project3/Noinjury --sample=noninjury2
Well, did you make sure cellranger was installed?
cellranger is for sure installed. A sample script using the same commands just different fastq paths worked fine but for some reason this script will not
Are you sure this
PATH
value is correct/network/rit/misc/scratch/single_cell/cellranger/cellranger=3.0.2
? Perhaps there is abin
directory that you missed adding, so thePATH
should really be/network/rit/misc/scratch/single_cell/cellranger/cellranger=3.0.2/bin
.cellranger
executable should be in that directory. I would check on that first.The path used has worked on previous scripts but I tried what you said. It still resulted in the same error.
That does not make logical sense based on error which said command not found. If a prior script worked with same
PATH
but if this is generating an error then you should verify that thecellranger
executable is there. Do you see a file listed if you tryls -lh /network/rit/misc/scratch/single_cell/cellranger/cellranger=3.0.2/bin/cellranger
?