diff --git a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/basty/Basty.scala b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/basty/Basty.scala index 9e73788f7027026f34960998ae8d48a9c1159db9..abd41b5ddbc97515cae8875716a5af937e8452cd 100644 --- a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/basty/Basty.scala +++ b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/basty/Basty.scala @@ -47,6 +47,11 @@ class Basty(val root: Configurable) extends QScript with MultiSampleQScript { val catVariantsSnps = Cat(this, refVariantSnps.variants :: samplesOutput.map(_._2.outputSnps.variants).toList, outputDir + "fastas/variant.snps_only.fasta") add(catVariantsSnps) + val catConsensus = Cat(this, refVariants.consensus :: samplesOutput.map(_._2.output.consensus).toList, outputDir + "fastas/consensus.fasta") + add(catConsensus) + val catConsensusSnps = Cat(this, refVariantSnps.consensus :: samplesOutput.map(_._2.outputSnps.consensus).toList, outputDir + "fastas/consensus.snps_only.fasta") + add(catConsensusSnps) + val catConsensusVariants = Cat(this, refVariants.consensusVariants :: samplesOutput.map(_._2.output.consensusVariants).toList, outputDir + "fastas/consensus.variant.fasta") add(catConsensusVariants) val catConsensusVariantsSnps = Cat(this, refVariantSnps.consensusVariants :: samplesOutput.map(_._2.outputSnps.consensusVariants).toList, outputDir + "fastas/consensus.variant.snps_only.fasta")