Question: The two mates do not align to a same transcript -- RSEM error

Skip to first unread message

Abhilash KANNAN

Jan 22, 2021, 11:06:16 AM1/22/21
to RSEM Users

Hello RSEM Community,

I am trying to to run rsem-calculate-expression using STAR to aligned reads.

the command that i use is:

for preparing the reference

rsem-prepare-reference --gtf GENOME_data/Homo_sapiens.GRCh38.90.gtf --star GENOME_data/Homo_sapiens.GRCh38.dna.genome.fa --star-path /users/data/applications/STAR-2.5.3a/bin/ -p 28 GENOME_data/rsem_ref

for Quantification with RSEM

rsem-calculate-expression --p 28 --star --star-gzipped-file --star-path /users/data/applications/STAR-2.5.3a/bin --paired-end --seed 5 --output-genome-bam --sort-bam-by-read-name --forward-prob=1.0 ../R1.read1.fastq.gz ../R1.read2.fastq.gz GENOME_data/rsem_ref/ RNASEQ_data/R1

Now i get this error:

Read A00182:360:HJM25DRXX:1:1102:8703:35728: the two mates do not align to the same transcript. RSEM does not support discordant alignments.

"rsem-parse-alignments GENOME_data/rsem-ref RNASEQ_data/R1.temp/R1 RNASEQ_data/R1.stat RNASEQ_data/R1.temp/R1.sorted.bam 3 -t ag XM" failed. Please check if you provide correct parameters/options for the pipeline

Any help on this is highly appreciated. Thank you very much.

Abhilash Kannan

Reply all
Reply to author
0 new messages