samtools / htsjdk

A Java API for high-throughput sequencing data (HTS) formats.
http://samtools.github.io/htsjdk/
283 stars 242 forks source link

CRAM should validate the SAM header against the reference being used #1244

Open cmnbroad opened 5 years ago

cmnbroad commented 5 years ago

When creating a CRAM file, the SAM header from the input is propagated to the output, but I don't think there is any checking to validate that the reference being used matches the sequence dictionary in the header. At a minimum we should issue a warning when this happens.

Somewhat related to https://github.com/samtools/htsjdk/issues/506.