Skip to content
Snippets Groups Projects
Commit 97d99695 authored by Sander Bollen's avatar Sander Bollen
Browse files

set sequence dicttionary for vepnormalizer

parent e7ed7e5f
No related branches found
No related tags found
No related merge requests found
...@@ -90,6 +90,11 @@ object VepNormalizer extends ToolCommand { ...@@ -90,6 +90,11 @@ object VepNormalizer extends ToolCommand {
logger.debug("Parsing header") logger.debug("Parsing header")
val new_infos = parseCsq(header) val new_infos = parseCsq(header)
header.setWriteCommandLine(true) header.setWriteCommandLine(true)
val writer = new AsyncVariantContextWriter(new VariantContextWriterBuilder().
setOutputFile(output).setReferenceDictionary(header.getSequenceDictionary)
build())
for (info <- new_infos) { for (info <- new_infos) {
val tmpheaderline = new VCFInfoHeaderLine(info, VCFHeaderLineCount.UNBOUNDED, VCFHeaderLineType.String, "A VEP annotation") val tmpheaderline = new VCFInfoHeaderLine(info, VCFHeaderLineCount.UNBOUNDED, VCFHeaderLineType.String, "A VEP annotation")
header.addMetaDataLine(tmpheaderline) header.addMetaDataLine(tmpheaderline)
...@@ -97,9 +102,7 @@ object VepNormalizer extends ToolCommand { ...@@ -97,9 +102,7 @@ object VepNormalizer extends ToolCommand {
logger.debug("Header parsing done") logger.debug("Header parsing done")
logger.debug("Writing header to file") logger.debug("Writing header to file")
val writer = new AsyncVariantContextWriter(new VariantContextWriterBuilder().
setOutputFile(output).
build())
writer.writeHeader(header) writer.writeHeader(header)
logger.debug("Wrote header to file") logger.debug("Wrote header to file")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment