Skip to content

samtools stats only checksums half the quality values #2187

Description

@daviesrob

Specifically, the first half on each alignment record, due to a cut and paste error from the line above.

As quality values are stored one per byte, the line should be:

    stats->checksum.quals += crc32(0L, qual, seq_len);

and not (seq_len+1)/2 as written.

If this is fixed, the change will need to be documented as quality checksums will change between samtools versions.

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Fields

No fields configured for issues without a type.

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions