Hi all,
I am trying to analyze my RNA-seq data with RSEM by using STAR aligner.
For this I downloaded the gtf file, hg19 genome and is using the following code
rsem-prepare-reference --gtf hg19-iGenomes.gtf --star /data/RSEM/hg19 hg19_rsem
based on RSEM usage "rsem-prepare-reference [options] reference_fasta_file(s) reference_name"
where "/data/RSEM/hg19" contains the .fa file (fasta files) for all the chromosomes and hg_rsem is the reference name or --outFileNamePrefix.
It is generating the output files hg19_rsem.n2g.idx.fa, .idx.fa, .seq, .transcript.fa, .chrlist, grp, ti and chrLength.txt, chrName.txt, chrNameLength.txt, chrStart.txt, genomeParameters.txt, hg19_rsemLog.out, hg19_rsem_STARtmp.
However, it is ending with:
Aug 22 13:34:46 ..... started STAR run
Aug 22 13:34:46 ... starting to generate Genome files
Aug 22 13:35:57 ... starting to sort Suffix Array. This may take a long time...
Aug 22 13:36:22 ... sorting Suffix Array chunks and saving them to disk...
"STAR --runThreadN 1 --runMode genomeGenerate --genomeDir . --genomeFastaFiles /data//RSEM/hg19/chr1.fa, ........ /data/RSEM/hg19/chrY.fa --sjdbGTFfile hg19-iGenomes.gtf --sjdbOverhang 100 --outFileNamePrefix hg19_rsem" failed! Plase check if you provide correct parameters/options for the pipeline!
I am unable to understand where I am going wrong. Please help me sort this out.
Thanks Debbie
