This is a test version of Biostars. For the public version, visit https://www.biostars.org.
STAR alignment error: quantmod not found

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 rnaseq alignment
• 3,099 views • link
updated by Devon Ryan Devon Ryan • written by prabin.dm prabin.dm

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?

• 1 views • link
written by Ram Ram

1 answer

-- genomeDir has a space in it.

• 1 views • link
written by Devon Ryan Devon Ryan

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

• 1 views • link
written by prabin.dm prabin.dm

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

• 1 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this