From 0b6b3fa8fc20461b41a368b82c62c2f9d7472b31 Mon Sep 17 00:00:00 2001 From: Justin Bedo Date: Wed, 19 Jun 2019 11:24:38 +1000 Subject: octopus: initial functions --- test-tnpair.nix | 4 ++++ 1 file changed, 4 insertions(+) (limited to 'test-tnpair.nix') diff --git a/test-tnpair.nix b/test-tnpair.nix index 76fcfab..f90e3a9 100644 --- a/test-tnpair.nix +++ b/test-tnpair.nix @@ -37,8 +37,10 @@ let processPair = { tumour, normal }: rec { alignments = mapAttrs (_: x: markdup (sort (fixmate (alignWithRG x.name x.files)))) { inherit normal tumour; }; variants = callVariants alignments; + octopusSomatic = octopus.callSomatic {} {inherit (alignments) normal; tumours = [ alignments.tumour ];}; glvariants = strelka.call {} (builtins.attrValues alignments); platypusVars = platypus.call {} (builtins.attrValues alignments); + octopusVars = octopus.call {} (builtins.attrValues alignments); shards = pipe [ (shard.fastQPair 2) (map (bwa.align {inherit ref;})) @@ -56,6 +58,8 @@ let (ln (facets.callCNV {} {vcf = tnpairResult.platypusVars; bams = with tnpairResult.alignments; [ normal tumour ];}) "facets") (ln cnvkitResults.cnvs "cnvkit") (ln cnvkitResults.plot "cnvkit.pdf") + (ln tnpairResult.octopusVars "octopus.vcf") + (ln tnpairResult.octopusSomatic "octopus-somatic.vcf") (ln tnpairResult.variants "strelka") (ln tnpairResult.glvariants "strelka-gl") (ln tnpairResult.variants.indels "strelka.indels.vcf") -- cgit v1.2.3