Hi All,
I'm mapping RNA-Seq reads onto the human transcriptome (gencode set) and I would to translate transcript coordinates in genomic coordinates.
I tried RSEM and, in particular, the tool rsem-tbam2gbam but I got the following error:
rsem-tbam2gbam: BamConverter.h:154: void BamConverter:rocess(): Assertion `cqname != ((char*)((b)->data))' failed.
and I don't know how to fix it.
Is there an alternative way?
Many Thanks,
Ernesto
I'm mapping RNA-Seq reads onto the human transcriptome (gencode set) and I would to translate transcript coordinates in genomic coordinates.
I tried RSEM and, in particular, the tool rsem-tbam2gbam but I got the following error:
rsem-tbam2gbam: BamConverter.h:154: void BamConverter:rocess(): Assertion `cqname != ((char*)((b)->data))' failed.
and I don't know how to fix it.
Is there an alternative way?
Many Thanks,
Ernesto
Comment