From e7a2705bb3818db7647c7a37ef675760c323bd26 Mon Sep 17 00:00:00 2001 From: Justin Bedo Date: Thu, 2 May 2019 07:45:02 +1000 Subject: samtools-merge: detect correct sorting --- test-tnpair.nix | 2 ++ 1 file changed, 2 insertions(+) (limited to 'test-tnpair.nix') diff --git a/test-tnpair.nix b/test-tnpair.nix index 1b5f347..f2e84ac 100644 --- a/test-tnpair.nix +++ b/test-tnpair.nix @@ -13,6 +13,7 @@ let alignWithRG = rg: bwa.align { inherit ref; flags = "-R'@RG\\tID:${rg}\\tSM:${rg}'";}; sort = samtools.sort {}; + nameSort = samtools.sort {nameSort = true;}; flagstat = samtools.flagstat {}; check-fastqc = fastqc.check {}; check-fastp = fastp.run {}; @@ -56,6 +57,7 @@ let (ln (gridss.call (with tnpairResult.alignments; [normal tumour])) "gridss2") (ln (gridss.callAndAssemble (with tnpairResult.alignments; [normal tumour])) "gridss3") (ln (samtools.merge {} [tnpairResult.alignments.tumour tnpairResult.alignments.normal]) "alignments/merged.bam") + (ln (samtools.merge {} [(nameSort tnpairResult.alignments.tumour) (nameSort tnpairResult.alignments.normal)]) "alignments/merged-namesorted.bam") (ln (samtools.view { outfmt = types.toCram; } (tnpairResult.alignments.tumour)) "alignments/${tnpair.tumour.name}.cram") #(ln (samtools.view { outfmt = types.toCram; } (tnpairResult.alignments.normal)) "alignments/${tnpair.normal.name}.cram") (ln (flagstat tnpairResult.alignments.tumour) "alignments/${tnpair.tumour.name}.flagstat") -- cgit v1.2.3