diff --git a/public/gentrap/src/main/scala/nl/lumc/sasc/biopet/pipelines/gentrap/measures/FragmentsPerGene.scala b/public/gentrap/src/main/scala/nl/lumc/sasc/biopet/pipelines/gentrap/measures/FragmentsPerGene.scala index 78224cc181f042083e87f759faf4acbc548f7906..bcfc417864cc16279e0a886306f0ebce1ca1816c 100644 --- a/public/gentrap/src/main/scala/nl/lumc/sasc/biopet/pipelines/gentrap/measures/FragmentsPerGene.scala +++ b/public/gentrap/src/main/scala/nl/lumc/sasc/biopet/pipelines/gentrap/measures/FragmentsPerGene.scala @@ -26,6 +26,8 @@ import org.broadinstitute.gatk.queue.QScript class FragmentsPerGene(val root: Configurable) extends QScript with Measurement with AnnotationGtf { 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 */ def biopetScript(): Unit = { val jobs = bamFiles.map { @@ -38,10 +40,6 @@ class FragmentsPerGene(val root: Configurable) extends QScript with Measurement job.output = new File(outputDir, s"$id.$name.counts") job.format = Option("bam") 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 }