diff --git a/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/VepNormalizer.scala b/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/VepNormalizer.scala index a10882b817fa36b461446afd63a9e50e49261e99..7945cfbf9b79260debff150c063fb033d1b80aba 100644 --- a/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/VepNormalizer.scala +++ b/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/VepNormalizer.scala @@ -90,6 +90,11 @@ object VepNormalizer extends ToolCommand { logger.debug("Parsing header") val new_infos = parseCsq(header) header.setWriteCommandLine(true) + + val writer = new AsyncVariantContextWriter(new VariantContextWriterBuilder(). + setOutputFile(output).setReferenceDictionary(header.getSequenceDictionary) + build()) + for (info <- new_infos) { val tmpheaderline = new VCFInfoHeaderLine(info, VCFHeaderLineCount.UNBOUNDED, VCFHeaderLineType.String, "A VEP annotation") header.addMetaDataLine(tmpheaderline) @@ -97,9 +102,7 @@ object VepNormalizer extends ToolCommand { logger.debug("Header parsing done") logger.debug("Writing header to file") - val writer = new AsyncVariantContextWriter(new VariantContextWriterBuilder(). - setOutputFile(output). - build()) + writer.writeHeader(header) logger.debug("Wrote header to file")