Skip to content
GitLab
Menu
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
e60b100f
Commit
e60b100f
authored
Apr 03, 2017
by
akaljuvee
Browse files
-
parent
fbc8e4ab
Changes
1
Hide whitespace changes
Inline
Side-by-side
shiva/src/main/scala/nl/lumc/sasc/biopet/pipelines/shiva/ShivaSvCalling.scala
View file @
e60b100f
...
...
@@ -16,7 +16,7 @@ package nl.lumc.sasc.biopet.pipelines.shiva
import
htsjdk.variant.vcf.VCFFileReader
import
nl.lumc.sasc.biopet.core.summary.
{
Summarizable
,
SummaryQScript
}
import
nl.lumc.sasc.biopet.core.
{
PipelineCommand
,
Reference
,
SampleLibraryTag
}
import
nl.lumc.sasc.biopet.core.
{
PipelineCommand
,
Reference
}
import
nl.lumc.sasc.biopet.extensions.Pysvtools
import
nl.lumc.sasc.biopet.pipelines.shiva.svcallers._
import
nl.lumc.sasc.biopet.utils.config.Configurable
...
...
@@ -28,7 +28,7 @@ import org.broadinstitute.gatk.queue.QScript
*
* Created by pjvan_thof on 2/26/15.
*/
class
ShivaSvCalling
(
val
parent
:
Configurable
)
extends
QScript
with
SummaryQScript
with
SampleLibraryTag
with
Reference
{
class
ShivaSvCalling
(
val
parent
:
Configurable
)
extends
QScript
with
SummaryQScript
with
Reference
{
qscript
=>
def
this
()
=
this
(
null
)
...
...
@@ -100,7 +100,7 @@ class ShivaSvCalling(val parent: Configurable) extends QScript with SummaryQScri
for
((
sample
,
mergedResultFile
)
<-
outputMergedVCFbySample
)
{
lazy
val
counts
=
getVariantCounts
(
mergedResultFile
,
ShivaSvCallingReport
.
histogramBinBoundaries
)
addSummarizable
(
new
Summarizable
{
def
summaryFiles
=
Map
.
empty
def
summaryFiles
=
Map
(
"output_vcf"
->
mergedResultFile
)
def
summaryStats
=
counts
},
"variantsBySizeAndType"
,
Some
(
sample
))
}
...
...
@@ -119,7 +119,7 @@ class ShivaSvCalling(val parent: Configurable) extends QScript with SummaryQScri
def
summarySettings
=
Map
(
"sv_callers"
->
configCallers
.
toList
)
/** Files for the summary */
def
summaryFiles
:
Map
[
String
,
File
]
=
outputMergedVCFbySample
++
(
if
(
inputBams
.
size
>
1
)
Map
(
"final_mergedvcf"
->
outputMergedVCF
)
else
Nil
)
def
summaryFiles
:
Map
[
String
,
File
]
=
if
(
inputBams
.
size
>
1
)
Map
(
"final_mergedvcf"
->
outputMergedVCF
)
else
Map
.
empty
def
getVariantCounts
(
vcfFile
:
File
,
breaks
:
Array
[
Int
])
:
Map
[
String
,
Any
]
=
{
val
delCounts
,
insCounts
,
dupCounts
,
invCounts
=
Array
.
fill
(
breaks
.
size
+
1
)
{
0
}
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a 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