Entering edit mode
3.1 years ago
Palgrave
▴
130
I am running multiple samples defined in files.txt. How should I load the genome correctly to avoid it being loaded for each iteration in the loop? I tried the following, but the --genomeLoad LoadAndKeep
required fastq-files to be loaded, in addition to the genome.
mapfile -s 1 -t files < files.txt
echo "${files[@]}"
STAR --genomeLoad LoadAndKeep --genomeDir /Genome/
for FILE in ${files[@]}; do
STAR --genomeLoad LoadAndKeep --quantMode GeneCounts ......;
done
STAR --genomeLoad Remove --genomeDir /Genome/