Question: STAR alignment error: quantmod not found
0
gravatar for prabin.dm
12 months ago by
prabin.dm160
USA/Amherst/Umass
prabin.dm160 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 • 505 views
ADD COMMENTlink modified 12 months ago by Devon Ryan93k • written 12 months ago by prabin.dm160

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 12 months ago by RamRS25k
2
gravatar for Devon Ryan
12 months ago by
Devon Ryan93k
Freiburg, Germany
Devon Ryan93k wrote:

-- genomeDir has a space in it.

ADD COMMENTlink written 12 months ago by Devon Ryan93k

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

ADD REPLYlink modified 12 months ago • written 12 months ago by prabin.dm160

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 12 months ago by Devon Ryan93k
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: 1752 users visited in the last hour