diff --git a/generate-indexes/src/main/scala/nl/lumc/sasc/biopet/pipelines/generateindexes/GenerateIndexes.scala b/generate-indexes/src/main/scala/nl/lumc/sasc/biopet/pipelines/generateindexes/GenerateIndexes.scala index 5041dbb559ac4885cf2e5fc9c69623daf2c0c0fe..1f39db6b2079bb8b0d6f2c51ddce08ab6d5fafd6 100644 --- a/generate-indexes/src/main/scala/nl/lumc/sasc/biopet/pipelines/generateindexes/GenerateIndexes.scala +++ b/generate-indexes/src/main/scala/nl/lumc/sasc/biopet/pipelines/generateindexes/GenerateIndexes.scala @@ -196,6 +196,21 @@ class GenerateIndexes(val root: Configurable) extends QScript with BiopetQScript outputConfig += "dbsnp" -> cv.out } + val gtfFile: Option[File] = genomeConfig.get("gtf_uri").map { gtfUri => + val curl = new Curl(this) + curl.url = gtfUri.toString + curl.output = new File(annotationDir, new File(curl.url).getName) + add(curl) + outputConfig += "annotation_gtf" -> curl.output + curl.output + } + + val refFlatFile: Option[File] = gtfFile.map { gtf => + val refFlat = new File(gtf + ".refFlat") + //TODO: gtf to refFlat conversion + refFlat + } + // Bwa index val bwaIndex = new BwaIndex(this) bwaIndex.reference = createLinks(new File(genomeDir, "bwa"))