yes, rsem already gives out the sample.genes.results the exact command is in the second comment. I know that if I add the transcript alignments it work. But what I actually want is to use the genome coordinates, since I use this alignment to filter reads. In the documentation of STAR it says it first uses the genomic alignment to to then match it to the transcriptome, its this step thats needed. Seems the solution is to just split them into fastq files and remap though.
using rsem sample.genome.STAR.genome.bam for calculating expression
Hello,
I have run rsem-calculate-expression on pair end sequencing. From the resulting BAM file mapped to the genome (sample.STAR.genome.bam), I would like to use that as in input to calculate expression of genes in the sample. When I try to use it as input in rsem-calculate-expression, it fails because the chromosome entries dont appear as a reference sequence( expected). Any help is appreciated!
1 answer
You need to post your command lines. Didn't running calculate-expression already give you the gene count files? And are you sure that genome alignments, and not transcriptome alignments are what's wanted as input?
See the examples here, for instance
https://ycl6.gitbook.io/rna-seq-data-analysis/rna-seq_analysis_workflow/quantification_using_rsem1
Log in to answer this question.
More posts like this
-
RSEM uses different gene lengths for each sample.
written by Apprentice
-
How to use RSEM to quantify genes and transcript expression from STAR outpute file (Aligned.toTrans…
written by otieno43 •
-
Change STAR Parameters in rsem-calculate-expression
written by patelk26
-
Error: RSEM Bam file for variant calling
written by kartikayprasad •
-
Error when using snakemake with RSEM
written by nattzy94 •
-
RSEM Error - rsem-calculate-expression
written by anoo •
-
RSEM with STAR aligment
written by prabin.dm
-
RSEM calculate expression troubleshoot
written by prabin.dm
-
RSEM error: RSEM currently does not support gapped alignments, sorry!
written by lfkhajavi •
-
The problem with RSEM
written by 498476732 •
If you could post what exact command you used and error you are getting will help more. One possibility may be reference fasta file used in generating bam files and using with rsem-calculate-expression are different in chromosome name.
Sorry for not being clear ! There is no "error", the
sample.STAR.genome.bamis part of the outputrsem-calculate-expression --star-output-genome-bam. I use this file to filter reads and want to reuse this file for gene quantification, so when I runrsem-calculate-expression -p 32 --paired-end --alignments my_bam.genome.bam GRCh38_gencode dedup_sampleI get the error (which makes sense, one is in genomic coordinates, and this requires it to be in transcriptome coordinates):
RSEM can not recognize reference sequence name chr1!From the documentation, it says STAR uses the genomic coordinate one, to establish concordance, this is the step im missing. I looked into running
STARwith--quantMode TranscriptomeSAM, the problem being that the input cant be a bam fileYou cannot use the STAR genome bam to generate counts. You need to FASTQ files (which you can get by extracting reads from the STAR aligned BAM file) as input to
rsem-calculate-expression.