Question: RNAseq error while running sem-calculate-expression
0
gravatar for shuksi1984
14 months ago by
shuksi198450
shuksi198450 wrote:

Hi,

I ran following command:

/path/to/rsem/rsem-calculate-expression -p 2 --output-genome-bam -star --star-path /path/to

/STAR_index_files --paired-end path/to/SRR925687_1.fastq path/to/SRR925687_2.fastq human.refseq.gff.gtf path/to/RNA.test

I got following error msg:

dFilesIn SRR925687_1.fastq SRR925687_2.fastq" failed! Plase check if you provide correct parameters/options for the pipeline!

STAR index folder consists of following files:

rwxrwxr-x 2 stuser stuser 4096 Apr 17 18:16 ./

drwxrwxr-x 7 stuser stuser 4096 Apr 18 11:12 ../

-rw-r--r-- 1 root root 1194 Apr 9 18:37 chrLength.txt

-rw-r--r-- 1 root root 3532 Apr 9 18:37 chrNameLength.txt

-rw-r--r-- 1 root root 2338 Apr 9 18:37 chrName.txt

-rw-r--r-- 1 root root 2108 Apr 9 18:37 chrStart.txt

-rw-r--r-- 1 root root 51816787 Apr 9 19:41 exonGeTrInfo.tab

-rw-r--r-- 1 root root 23213627 Apr 9 19:41 exonInfo.tab

-rw-r--r-- 1 root root 380243 Apr 9 19:41 geneInfo.tab

-rw-r--r-- 1 root root 3202741093 Apr 9 21:34 Genome

-rw-r--r-- 1 root root 858 Apr 9 21:25 genomeParameters.txt

-rw-r--r-- 1 root root 12693443564 Apr 9 22:47 SA

-rw-r--r-- 1 root root 1565873619 Apr 9 22:54 SAindex

-rw-r--r-- 1 root root 9414412 Apr 9 19:41 sjdbInfo.txt

-rw-r--r-- 1 root root 10724637 Apr 9 19:41 sjdbList.fromGTF.out.tab

-rw-r--r-- 1 root root 10723973 Apr 9 19:41 sjdbList.out.tab

-rw-r--r-- 1 root root 7215534 Apr 9 19:41 transcriptInfo.tab

I have provided correct path to the files generated by STAR indexing. Dont understand the problem.

ADD COMMENTlink modified 12 months ago by Biostar ♦♦ 20 • written 14 months ago by shuksi198450
1

It's likely that human.refseq.gff.gtf is an incorrect reference name. Did you use that with rsem-prepare-reference? BTW, your error message looks truncated.

ADD REPLYlink modified 14 months ago • written 14 months ago by Devon Ryan90k

Yes, human.refseq.gff.gtf was generated with rsem-prepare-reference. The same file was used during indexing with STAR.

ADD REPLYlink written 14 months ago by shuksi198450

Following shows the complete error.

sh: 1: /data/ngs/mona/RNAseq/genome/genome_file/STAR: not found
"/data/ngs/mona/RNAseq/genome/genome_file/STAR --genomeDir .  --outSAMunmapped Within  --outFilterType BySJout  --outSAMattributes NH HI AS NM MD  --outFilterMultimapNmax 20  --outFilterMismatchNmax 999  --outFilterMismatchNoverLmax 0.04  --alignIntronMin 20  --alignIntronMax 1000000  --alignMatesGapMax 1000000  --alignSJoverhangMin 8  --alignSJDBoverhangMin 1  --sjdbScore 1  --runThreadN 2  --genomeLoad NoSharedMemory  --outSAMtype BAM Unsorted  --quantMode TranscriptomeSAM  --outSAMheaderHD \@HD VN:1.4 SO:unsorted  --outFileNamePrefix RNA.test.temp/RNA.test  --readFilesIn SRR925687_1.fastq SRR925687_2.fastq" failed! Plase check if you provide correct parameters/options for the pipeline!
stuser@beta:/data/ngs/mona/RNAseq/genome$ cd ../../../
ADD REPLYlink modified 14 months ago by Devon Ryan90k • written 14 months ago by shuksi198450
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: 1885 users visited in the last hour