Skip to content
GitLab
Explore
Sign in
Primary navigation
Search or go to…
Project
B
biopet.biopet
Manage
Activity
Members
Code
Repository
Branches
Commits
Tags
Repository graph
Compare revisions
Build
Pipelines
Jobs
Pipeline schedules
Artifacts
Deploy
Releases
Container Registry
Model registry
Operate
Environments
Analyze
Contributor analytics
CI/CD analytics
Repository analytics
Model experiments
Help
Help
Support
GitLab documentation
Compare GitLab plans
Community forum
Contribute to GitLab
Provide feedback
Keyboard shortcuts
?
Snippets
Groups
Projects
Show more breadcrumbs
Mirrors
biopet.biopet
Commits
c4c5dc60
Commit
c4c5dc60
authored
9 years ago
by
bow
Browse files
Options
Downloads
Patches
Plain Diff
Set default htseq-count sorting order to pos
parent
0906f14b
No related branches found
No related tags found
No related merge requests found
Changes
1
Hide whitespace changes
Inline
Side-by-side
Showing
1 changed file
public/gentrap/src/main/scala/nl/lumc/sasc/biopet/pipelines/gentrap/measures/FragmentsPerGene.scala
+2
-4
2 additions, 4 deletions
.../biopet/pipelines/gentrap/measures/FragmentsPerGene.scala
with
2 additions
and
4 deletions
public/gentrap/src/main/scala/nl/lumc/sasc/biopet/pipelines/gentrap/measures/FragmentsPerGene.scala
+
2
−
4
View file @
c4c5dc60
...
@@ -26,6 +26,8 @@ import org.broadinstitute.gatk.queue.QScript
...
@@ -26,6 +26,8 @@ import org.broadinstitute.gatk.queue.QScript
class
FragmentsPerGene
(
val
root
:
Configurable
)
extends
QScript
with
Measurement
with
AnnotationGtf
{
class
FragmentsPerGene
(
val
root
:
Configurable
)
extends
QScript
with
Measurement
with
AnnotationGtf
{
def
mergeArgs
=
MergeArgs
(
idCols
=
List
(
1
),
valCol
=
2
,
numHeaderLines
=
0
,
fallback
=
"0"
)
def
mergeArgs
=
MergeArgs
(
idCols
=
List
(
1
),
valCol
=
2
,
numHeaderLines
=
0
,
fallback
=
"0"
)
override
def
fixedValues
:
Map
[
String
,
Any
]
=
Map
(
"htseqcount"
->
Map
(
"order"
->
"pos"
))
/** Pipeline itself */
/** Pipeline itself */
def
biopetScript
()
:
Unit
=
{
def
biopetScript
()
:
Unit
=
{
val
jobs
=
bamFiles
.
map
{
val
jobs
=
bamFiles
.
map
{
...
@@ -38,10 +40,6 @@ class FragmentsPerGene(val root: Configurable) extends QScript with Measurement
...
@@ -38,10 +40,6 @@ class FragmentsPerGene(val root: Configurable) extends QScript with Measurement
job
.
output
=
new
File
(
outputDir
,
s
"$id.$name.counts"
)
job
.
output
=
new
File
(
outputDir
,
s
"$id.$name.counts"
)
job
.
format
=
Option
(
"bam"
)
job
.
format
=
Option
(
"bam"
)
add
(
job
)
add
(
job
)
// We are forcing the sort order to be ID-sorted, since HTSeq-count often chokes when using position-sorting due
// to its buffer not being large enough.
//TODO: ID sorting job
//job.order = Option("name")
id
->
job
id
->
job
}
}
...
...
This diff is collapsed.
Click to expand it.
Preview
0%
Loading
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!
Save comment
Cancel
Please
register
or
sign in
to comment