Question: RSEM error: RSEM can not recognize reference sequence name chrM!
2
gravatar for moxu
2.5 years ago by
moxu430
moxu430 wrote:

I have a set of BAM files from somewhere else which I don't know how they were generated. I have RSEM reference files generated with RSEM (bowtie, gtf, genomic fa). When running rsem-calcuate-expression, I got the following:

rsem-calculate-expression --bam mysample1.bam ~/softspace/reference/RSEM_bowtie_hg19_genes_reference test
rsem-parse-alignments /ark/home/mx010/softspace/reference/RSEM_bowtie_hg19_genes_reference test.temp/test test.stat/test KP9-1_Ctl2.bam 1 -tag XM
Warning: The SAM/BAM file declares less reference sequences (35) than RSEM knows (41970)!
RSEM can not recognize reference sequence name chrM!

The same reference files work well for BAM files generated by HISAT/Cufflinks.

Thanks in advance for your help!

ADD COMMENTlink written 2.5 years ago by moxu430

Are you using a BAM file with genomic alignments?

ADD REPLYlink written 2.5 years ago by Devon Ryan89k

Not sure. They are supposed to be RNAseq data. If the BAM files were with genomic alignments, what should I do?

Thanks!

ADD REPLYlink written 2.5 years ago by moxu430

Redo the alignments to the transcriptome that you're using with RSEM

ADD REPLYlink written 2.5 years ago by Devon Ryan89k

That would be nice but we don't have the fastq files. All we have are the BAM files.

ADD REPLYlink written 2.5 years ago by moxu430

You can convert them back to fastq with samtools bam2fq.

ADD REPLYlink written 2.5 years ago by Devon Ryan89k
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: 1165 users visited in the last hour