Hello All! I am trying to map the RNA seq data with the reference genome using STAR aligner tool in cluster. I used the reference fasta file and gtf file to create the genomeDirectory that has the genome index. Now, I want to map the RNA seq data using that star index, however I get the following error:
**EXITING: FATAL INPUT ERROR: empty value for parameter "genomeDir" in input "Command-Line"
SOLUTION: use non-empty value for this parameter.****
Could you please help me with this? Here is the batch script that I ran to do the mapping. genom_dir which is star_index is the directory that has all the files created from the star index step.
#!/bin/bash
#
# job standard output will go to the file slurm-%j.out (where %j is the job ID)
#SBATCH --time=3:00:00 # walltime limit (HH:MM:SS)
#SBATCH --nodes=1 # number of nodes
#SBATCH --ntasks-per-node=16 # 16 processor core(s) per node
#SBATCH --mem=100G # maximum memory per node
#SBATCH --mail-user=username@college.edu # email address
#SBATCH --mail-type=END
#SBATCH --output=/output_directoy_path/map_out/"s-%j-0716_temp_map.out" # job standard output file (%j replaced by job id)
# LOAD MODULES, INSERT CODE, AND RUN YOUR PROGRAMS HERE
# Load modules
#=====================================================#
#--| STAR
module load star/2.7.6a-wq7xoea
#=====================================================#
# Assign varialbles
#=====================================================#
#--| PATH
pKS="/working directory path/directory_name"
genom_dir= ${pKS}/"star_index"
inputFile = ${pKS}/"job"
#--| FILE
finput=${inputFile}/"temp.fasta"
#=====================================================#
# Run
#=====================================================#
#--| Make Mapping
STAR \
--runThreadN 8 --runMode genomeGenerate \
--genomeDir ${genom_dir} \
--readFilesIn ${finput} \
--quantMode GeneCounts \
--outFilterMismatchNmax 0
~