Question: STAR alignment error: quantmod not found
0
gravatar for prabin.dm
10 weeks ago by
prabin.dm30
prabin.dm30 wrote:

Hi all,

I am using STAR for aligment and gettting an error. I can't figure out why?

star 2.5.3a is located under /share/pkg/star/2.5.3a
/home/.lsbatch/1544131842.9569914.shell: line 22: --quantMode: command not found

Here is my code

module load star/2.5.3a

STAR -- genomeDir star_genome_mm10 \
        --readFilesCommand zcat \
        --readFilesIn L001_R1_001.fastq.gz \
        --outSAMtype BAM SortedByCoordinate \
        --outSAMunmapped Within \
        --twopassMode Basic \
        --outFilterMultimapNmax 1 \ 
        --quantMode TranscriptomeSAM GeneCounts \
        --runThreadN 6 \
        --outFileNamePrefix /STAR_output/test

Another issue I am having is, the output files are not stored in the STAR_ouput directory. It is stored in the directory where the code is executed ?

star alignment rnaseq • 181 views
ADD COMMENTlink modified 10 weeks ago by Devon Ryan88k • written 10 weeks ago by prabin.dm30

Maybe change this

--outFileNamePrefix /STAR_output/test

to

--outFileNamePrefix $PWD/STAR_output/test

unless there's actually a dir named STAR_output with proper permissions in the root directory?

ADD REPLYlink written 10 weeks ago by RamRS20k
2
gravatar for Devon Ryan
10 weeks ago by
Devon Ryan88k
Freiburg, Germany
Devon Ryan88k wrote:

-- genomeDir has a space in it.

ADD COMMENTlink written 10 weeks ago by Devon Ryan88k

correcting that resulted in the same output. Should I have the rsem-prepare-reference run before I use --quantMode TranscriptomeSAM ?

ADD REPLYlink modified 10 weeks ago • written 10 weeks ago by prabin.dm30

Your command has a random space or a missing \ somewhere that you're not showing us. Write it all on one line without any \.

ADD REPLYlink written 10 weeks ago by Devon Ryan88k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1265 users visited in the last hour