Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Mirrors
biopet.biopet
Commits
ff0c3a5e
Commit
ff0c3a5e
authored
Mar 20, 2015
by
Peter van 't Hof
Browse files
Added some errors when there is a field conflict
parent
2c3cf84b
Changes
1
Hide whitespace changes
Inline
Side-by-side
public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/VCFWithVCF.scala
View file @
ff0c3a5e
...
...
@@ -53,6 +53,11 @@ object VCFWithVCF extends ToolCommand {
val
secondHeader
=
secondaryReader
.
getFileHeader
for
(
x
<-
commandArgs
.
fields
)
{
if
(
header
.
hasInfoLine
(
x
.
outputField
))
throw
new
IllegalArgumentException
(
"Field '"
+
x
.
outputField
+
"' already exist in input vcf"
)
if
(!
secondHeader
.
hasInfoLine
(
x
.
inputField
))
throw
new
IllegalArgumentException
(
"Field '"
+
x
.
inputField
+
"' does not exist in secondary vcf"
)
val
oldHeaderLine
=
secondHeader
.
getInfoHeaderLine
(
x
.
inputField
)
val
newHeaderLine
=
new
VCFInfoHeaderLine
(
x
.
outputField
,
VCFHeaderLineCount
.
UNBOUNDED
,
...
...
@@ -65,7 +70,7 @@ object VCFWithVCF extends ToolCommand {
writer
.
writeHeader
(
header
)
var
idx
=
0
for
(
record
:
VariantContext
<-
reader.iterator
()
)
{
for
(
record
<-
reader
)
{
val
secondaryRecords
=
if
(
commandArgs
.
matchAllele
)
{
secondaryReader
.
query
(
record
.
getChr
,
record
.
getStart
,
record
.
getEnd
).
toList
.
filter
(
x
=>
record
.
getAlternateAlleles
.
exists
(
x
.
hasAlternateAllele
(
_
)))
...
...
@@ -97,6 +102,7 @@ object VCFWithVCF extends ToolCommand {
logger
.
info
(
s
"""Processed $idx records"""
)
}
}
logger
.
info
(
s
"""Processed $idx records"""
)
logger
.
debug
(
"Closing readers"
)
writer
.
close
()
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment