Skip to content

Picard vs Samtools queryname sort incompatibility is a common user pain point #8

@clintval

Description

@clintval

Picard/fgbio/htsjdk aligned files (such as to a wildtype genome) are of a different sort order to the often-used samtools approach of using samtools fastq and aligning with bwa. Thus this exception is hit when trying to neodisambiguate:

require(
templates.map(_.name).distinct.length <= 1,
"Templates with different names found! This can only occur if your SAM sources are queryname sorted using"
+ " different implementations, such as with Picard tools versus Samtools. If you have encountered this"
+ " exception, then please alert the maintainer!"
)

It might be nice to have a fallback mode that re-sorts the SAM streams if they:

  • Say they are queryname sorted
  • And appear out of order (parsimony is different sort implementations)

Metadata

Metadata

Assignees

Labels

No labels
No labels

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions