STAR mapping with genomeLoad issue bash script
0
0
Entering edit mode
18 months ago

Hello,

I am trying map RNA seq data in STAR bash script with the below command when I am running script its going in sleep mode

script which I am running

 date
STAR --genomeLoad LoadAndExit --genomeDir /media/index/
for f in `ls -1 *_1.fastq.gz | sed 's/_1.fastq.gz//' `
do
echo  ""${f}
STAR --readFilesIn ${f}_1.fastq.gz ${f}_2.fastq.gz\
     --outFileNamePrefix ${f}. \
     --readFilesCommand zcat \
     --runThreadN 8 --outSAMtype BAM \
     --genomeLoad LoadAndKeep \
     --genomeDir /media/index/ \
     --sjdbGTFfile /media/mm10.gtf

STAR --genomeLoad Remove --genomeDir /media/index/
done 
date

please help me to figure it out what is the mistake in the script

RNA-Seq STAR • 489 views
ADD COMMENT
1
Entering edit mode

No answer validated in 2 years?

ADD REPLY

Login before adding your answer.

Traffic: 1534 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6