-
Notifications
You must be signed in to change notification settings - Fork 1
Open
Description
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:
neodisambiguate/neodisambiguate/src/io/cvbio/neodisambiguate/bam/Bams.scala
Lines 68 to 73 in db7a512
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