Getting samtools cram-to-bam error: [E::cram_get_ref] Failed to populate reference for id?

Question

[claudiadast]

I am trying to convert a cram to a bam using the following samtools command set-up:

samtools view -h -b -@ 16 file.cram -o file.bam -T reference.fa 2>&1 | tee -a log.txt

However, when I try this, I get the following error:

[E::cram_get_ref] Failed to populate reference for id 76
[E::cram_get_ref] Failed to populate reference for id 102
[E::cram_get_ref] Failed to populate reference for id 129
.
.
.
.
[continues for over 20 lines like this]

I'm not sure why I am getting this error. I've made sure to use the same reference file that was used to generate the cram (I was using a different reference file before) and have also included the reference.dict file.

It's also noteworthy that, although I'm getting this error, the cram file still gets converted to a bam and downstream processes like GATK are still able to run on the bam. It's just concerning to be getting such an error message. Any thoughts?

[GenoMax]

Hello claudiadast!

It appears that your post has been cross-posted to another site: https://stackoverflow.com/questions/48938569/getting-samtools-cram-to-bam-error-ecram-get-ref-failed-to-populate-referen

This is typically not recommended as it runs the risk of annoying people in both communities.