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
e79bb410
Commit
e79bb410
authored
Oct 26, 2016
by
Sander Bollen
Browse files
Merge branch 'feature-single_bed_sort' into 'develop'
For for BIOPET-389 See merge request !464
parents
e6199e01
abb1f25e
Changes
1
Hide whitespace changes
Inline
Side-by-side
bammetrics/src/main/scala/nl/lumc/sasc/biopet/pipelines/bammetrics/BamMetrics.scala
View file @
e79bb410
...
...
@@ -166,11 +166,15 @@ class BamMetrics(val root: Configurable) extends QScript
addSummarizable
(
biopetFlagstatLoose
,
targetName
+
"_biopet_flagstat_loose"
)
add
(
new
BiopetFifoPipe
(
this
,
List
(
biLoose
,
biopetFlagstatLoose
)))
val
sorter
=
new
BedtoolsSort
(
this
)
sorter
.
input
=
intervals
.
bed
sorter
.
output
=
swapExt
(
targetDir
,
intervals
.
bed
,
".bed"
,
".sorted.bed"
)
add
(
sorter
)
val
bedCov
=
BedtoolsCoverage
(
this
,
sorter
.
output
,
inputBam
,
depth
=
true
)
val
sortedBed
=
BamMetrics
.
sortedbedCache
.
getOrElse
(
intervals
.
bed
,
{
val
sorter
=
new
BedtoolsSort
(
this
)
sorter
.
input
=
intervals
.
bed
sorter
.
output
=
swapExt
(
targetDir
,
intervals
.
bed
,
".bed"
,
".sorted.bed"
)
add
(
sorter
)
BamMetrics
.
sortedbedCache
+=
intervals
.
bed
->
sorter
.
output
sorter
.
output
})
val
bedCov
=
BedtoolsCoverage
(
this
,
sortedBed
,
inputBam
,
depth
=
true
)
val
covStats
=
CoverageStats
(
this
,
targetDir
,
inputBam
.
getName
.
stripSuffix
(
".bam"
)
+
".coverage"
)
covStats
.
title
=
Some
(
"Coverage Plot"
)
covStats
.
subTitle
=
Some
(
s
"for file '$targetName.bed'"
)
...
...
@@ -198,4 +202,6 @@ object BamMetrics extends PipelineCommand {
bamMetrics
.
biopetScript
()
bamMetrics
}
private
var
sortedbedCache
:
Map
[
File
,
File
]
=
Map
()
}
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