diff options
author | l-d-s <distefano.l@wehi.edu.au> | 2018-12-11 14:17:22 +1100 |
---|---|---|
committer | l-d-s <distefano.l@wehi.edu.au> | 2018-12-11 14:17:22 +1100 |
commit | f0357e09b8c473b55d22c4632692fcc79021570d (patch) | |
tree | 9c92993d1f0c788a162159e6440ea958eae3e35d /test-tnpair.nix | |
parent | 6742864b8f92f5d7d56da4901486b7e83a271e23 (diff) | |
parent | 10bdecf6c1338d7f531ddf7b41da14dfe4a4ac33 (diff) |
Merge branch 'master' of https://github.com/PapenfussLab/bionix
Diffstat (limited to 'test-tnpair.nix')
-rw-r--r-- | test-tnpair.nix | 11 |
1 files changed, 7 insertions, 4 deletions
diff --git a/test-tnpair.nix b/test-tnpair.nix index a6ea163..b01d61e 100644 --- a/test-tnpair.nix +++ b/test-tnpair.nix @@ -21,7 +21,9 @@ let fetchfq = attrs: types.tagFiletype (types.filetype.fq {}) (fetchlocal attrs); fetchfa = attrs: types.tagFiletype (types.filetype.fa {}) (fetchlocal attrs); - alignWithRG = rg: bwa.align { ref = fetchfa ./example/ref.fa; flags = "-R'@RG\\tID:${rg}\\tSM:${rg}'";}; + ref = fetchfa ./example/ref.fa; + + alignWithRG = rg: bwa.align { inherit ref; flags = "-R'@RG\\tID:${rg}\\tSM:${rg}'";}; sort = samtools.sort {}; flagstat = samtools.flagstat {}; check = fastqc.check {}; @@ -35,7 +37,7 @@ let }; normal = {name = "mysample2"; files = { input1 = fetchfq ./example/sample2-1.fq; - input2 = fetchfq ./example/sample2-1.fq; + input2 = fetchfq ./example/sample2-2.fq; }; }; }; @@ -53,8 +55,9 @@ let mkdir $out ln -s ${tnpairResult.variants} $out/strelka mkdir $out/alignments - ln -s ${gridss.callVariants {} (with tnpairResult.alignments; [tumour])} $out/gridss - ln -s ${gridss.call (with tnpairResult.alignments; [tumour])} $out/gridss2 + ln -s ${bowtie.align {inherit ref;} tnpair.normal.files} $out/alignments/bowtie-normal.bam + ln -s ${gridss.callVariants {} (with tnpairResult.alignments; [normal tumour])} $out/gridss + ln -s ${gridss.call (with tnpairResult.alignments; [normal tumour])} $out/gridss2 ln -s ${samtools.merge {} [tnpairResult.alignments.tumour tnpairResult.alignments.normal]} $out/alignments/merged.bam ln -s ${samtools.view { outfmt = types.toCram; } (tnpairResult.alignments.tumour)} $out/alignments/${tnpair.tumour.name}.cram ln -s ${samtools.view { outfmt = types.toCram; } (tnpairResult.alignments.normal)} $out/alignments/${tnpair.normal.name}.cram |