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
43cd6298
Commit
43cd6298
authored
Mar 17, 2016
by
Wai Yi Leung
Browse files
Merge branch 'fix-pValue' into 'develop'
Added wigle and tdf files to summary See merge request !361
parents
18184194
48b7822f
Changes
2
Hide whitespace changes
Inline
Side-by-side
public/bam2wig/src/main/scala/nl/lumc/sasc/biopet/pipelines/bamtobigwig/Bam2Wig.scala
View file @
43cd6298
...
...
@@ -38,6 +38,10 @@ class Bam2Wig(val root: Configurable) extends QScript with BiopetQScript {
inputFiles
:+=
new
InputFile
(
bamFile
)
}
def
outputWigleFile
=
new
File
(
outputDir
,
bamFile
.
getName
+
".wig"
)
def
outputTdfFile
=
new
File
(
outputDir
,
bamFile
.
getName
+
".tdf"
)
def
outputBwFile
=
new
File
(
outputDir
,
bamFile
.
getName
+
".bw"
)
def
biopetScript
()
:
Unit
=
{
val
bs
=
new
BamToChromSizes
(
this
)
bs
.
bamFile
=
bamFile
...
...
@@ -48,14 +52,14 @@ class Bam2Wig(val root: Configurable) extends QScript with BiopetQScript {
val
igvCount
=
new
IGVToolsCount
(
this
)
igvCount
.
input
=
bamFile
igvCount
.
genomeChromSizes
=
bs
.
chromSizesFile
igvCount
.
wig
=
Some
(
new
File
(
outputDir
,
bamFile
.
getName
+
".wig"
)
)
igvCount
.
tdf
=
Some
(
new
File
(
outputDir
,
bamFile
.
getName
+
".tdf"
)
)
igvCount
.
wig
=
Some
(
outputWigleFile
)
igvCount
.
tdf
=
Some
(
outputTdfFile
)
add
(
igvCount
)
val
wigToBigWig
=
new
WigToBigWig
(
this
)
wigToBigWig
.
inputWigFile
=
igvCount
.
wig
.
get
wigToBigWig
.
inputChromSizesFile
=
bs
.
chromSizesFile
wigToBigWig
.
outputBigWig
=
new
File
(
outputDir
,
bamFile
.
getName
+
".bw"
)
wigToBigWig
.
outputBigWig
=
outputBwFile
add
(
wigToBigWig
)
}
}
...
...
public/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala
View file @
43cd6298
...
...
@@ -109,7 +109,14 @@ class Mapping(val root: Configurable) extends QScript with SummaryQScript with S
/** File to add to the summary */
def
summaryFiles
:
Map
[
String
,
File
]
=
Map
(
"output_bam"
->
finalBamFile
,
"input_R1"
->
inputR1
,
"reference"
->
referenceFasta
())
++
(
if
(
inputR2
.
isDefined
)
Map
(
"input_R2"
->
inputR2
.
get
)
else
Map
())
(
if
(
inputR2
.
isDefined
)
Map
(
"input_R2"
->
inputR2
.
get
)
else
Map
())
++
(
bam2wig
match
{
case
Some
(
b
)
=>
Map
(
"output_wigle"
->
b
.
outputWigleFile
,
"output_tdf"
->
b
.
outputTdfFile
,
"output_bigwig"
->
b
.
outputBwFile
)
case
_
=>
Map
()
})
/** Settings to add to summary */
def
summarySettings
=
Map
(
...
...
@@ -270,12 +277,15 @@ class Mapping(val root: Configurable) extends QScript with SummaryQScript with S
add
(
gears
)
}
if
(
config
(
"generate_wig"
,
default
=
false
).
asBoolean
)
addAll
(
Bam2Wig
(
this
,
finalBamFile
).
functions
)
bam2wig
.
foreach
(
add
(
_
))
addSummaryJobs
()
}
protected
lazy
val
bam2wig
=
if
(
config
(
"generate_wig"
,
default
=
false
))
{
Some
(
Bam2Wig
(
this
,
finalBamFile
))
}
else
None
/** Add bwa aln jobs */
def
addBwaAln
(
R1
:
File
,
R2
:
Option
[
File
],
output
:
File
)
:
File
=
{
val
bwaAlnR1
=
new
BwaAln
(
this
)
...
...
Write
Preview
Markdown
is supported
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