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
4785e578
Commit
4785e578
authored
May 20, 2014
by
Peter van 't Hof
Browse files
- Fixed input field on markduplicates in mapping module
parent
205355eb
Changes
1
Hide whitespace changes
Inline
Side-by-side
mapping/src/main/java/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala
View file @
4785e578
...
...
@@ -84,11 +84,11 @@ class Mapping(private var globalConfig: Config) extends QScript {
val
bwaCommand
=
new
Bwa
(
config
)
{
R1
=
fastq_R1
;
if
(
paired
)
R2
=
fastq_R2
;
RG
=
getReadGroup
;
output
=
new
File
(
outputDir
+
outputName
+
".sam"
)
}
add
(
bwaCommand
)
bamFile
=
add
AddOrReplaceReadGroups
(
List
(
bwaCommand
.
output
),
swapExt
(
outputDir
,
bwaCommand
.
output
,
".sam"
,
".bam"
),
outputDir
)
bamFile
=
add
SortSam
(
List
(
bwaCommand
.
output
),
swapExt
(
outputDir
,
bwaCommand
.
output
,
".sam"
,
".bam"
),
outputDir
)
}
else
if
(
aligner
==
"star"
)
{
val
starCommand
=
new
Star
(
config
)
{
R1
=
fastq_R1
;
if
(
paired
)
R2
=
fastq_R2
;
this
.
outputDir
=
outputDir
}
add
(
starCommand
)
bamFile
=
add
SortSam
(
List
(
starCommand
.
outputSam
),
swapExt
(
outputDir
,
starCommand
.
outputSam
,
".sam"
,
".bam"
),
outputDir
)
bamFile
=
add
AddOrReplaceReadGroups
(
List
(
starCommand
.
outputSam
),
swapExt
(
outputDir
,
starCommand
.
outputSam
,
".sam"
,
".bam"
),
outputDir
)
}
if
(!
skipMarkduplicates
)
bamFile
=
addMarkDuplicates
(
List
(
bamFile
),
swapExt
(
outputDir
,
bamFile
,
".bam"
,
".dedup.bam"
),
outputDir
)
...
...
@@ -110,12 +110,12 @@ class Mapping(private var globalConfig: Config) extends QScript {
def
addAddOrReplaceReadGroups
(
inputSam
:
List
[
File
],
outputFile
:
File
,
dir
:
String
)
:
File
=
{
val
addOrReplaceReadGroups
=
new
AddOrReplaceReadGroups
{
input
=
inputSam
output
=
outputFile
createIndex
=
true
memoryLimit
=
2
nCoresRequest
=
2
jobResourceRequests
:+=
"h_vmem=4G"
this
.
input
=
inputSam
this
.
output
=
outputFile
this
.
createIndex
=
true
this
.
memoryLimit
=
2
this
.
nCoresRequest
=
2
this
.
jobResourceRequests
:+=
"h_vmem=4G"
this
.
RGLB
=
RGLB
this
.
RGPL
=
RGPL
...
...
@@ -127,15 +127,15 @@ class Mapping(private var globalConfig: Config) extends QScript {
return
addOrReplaceReadGroups
.
output
}
def
addMarkDuplicates
(
inputBams
:
List
[
File
],
outputFile
:
File
,
dir
:
String
)
:
File
=
{
def
addMarkDuplicates
(
input
_
Bams
:
List
[
File
],
outputFile
:
File
,
dir
:
String
)
:
File
=
{
val
markDuplicates
=
new
MarkDuplicates
{
input
=
inputBams
output
=
outputFile
REMOVE_DUPLICATES
=
false
metrics
=
swapExt
(
dir
,
outputFile
,
".bam"
,
".metrics"
)
outputIndex
=
swapExt
(
dir
,
this
.
output
,
".bam"
,
".bai"
)
memoryLimit
=
2
jobResourceRequests
:+=
"h_vmem=4G"
this
.
input
=
input
_
Bams
this
.
output
=
outputFile
this
.
REMOVE_DUPLICATES
=
false
this
.
metrics
=
swapExt
(
dir
,
outputFile
,
".bam"
,
".metrics"
)
this
.
outputIndex
=
swapExt
(
dir
,
this
.
output
,
".bam"
,
".bai"
)
this
.
memoryLimit
=
2
this
.
jobResourceRequests
:+=
"h_vmem=4G"
}
add
(
markDuplicates
)
...
...
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