From 86b1b192501eecc8a4e566795327772c32b1633d Mon Sep 17 00:00:00 2001 From: Justin Bedo Date: Wed, 5 Dec 2018 10:48:05 +1100 Subject: bowtie: init --- test-tnpair.nix | 11 +++++++---- 1 file changed, 7 insertions(+), 4 deletions(-) (limited to 'test-tnpair.nix') 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 -- cgit v1.2.3