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

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 22 months ago by RamRS30k
2
gravatar for Devon Ryan
22 months ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

-- genomeDir has a space in it.

ADD COMMENTlink written 22 months ago by Devon Ryan96k

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

ADD REPLYlink modified 22 months ago • written 22 months ago by prabin.dm180

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 22 months ago by Devon Ryan96k
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: 1551 users visited in the last hour