Commit f6598ac4 authored by Peter van 't Hof's avatar Peter van 't Hof
Browse files

Fix dict check

parent 459afe7a
......@@ -132,7 +132,18 @@ trait MultisampleMappingTrait extends MultiSampleQScript
val header = inputSam.getFileHeader
val readGroups = header.getReadGroups
val referenceFile = new FastaSequenceFile(referenceFasta(), true)
val dictOke = header.getSequenceDictionary == referenceFile.getSequenceDictionary
val dictOke: Boolean = {
var oke = true
try {
header.getSequenceDictionary.assertSameDictionary(referenceFile.getSequenceDictionary)
} catch {
case e: AssertionError => {
println(e.getMessage)
oke = false
}
}
oke
}
inputSam.close()
referenceFile.close()
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment