Question: STAR alignment error: quantmod not found
0
gravatar for prabin.dm
6 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 • 306 views
ADD COMMENTlink modified 6 months ago by Devon Ryan90k • written 6 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 6 months ago by RamRS22k
2
gravatar for Devon Ryan
6 months ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

-- genomeDir has a space in it.

ADD COMMENTlink written 6 months ago by Devon Ryan90k

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

ADD REPLYlink modified 6 months ago • written 6 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 6 months ago by Devon Ryan90k
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: 1607 users visited in the last hour