diff options
author | Justin Bedo <cu@cua0.org> | 2021-07-28 08:07:32 +1000 |
---|---|---|
committer | Justin Bedo <cu@cua0.org> | 2021-07-28 09:17:56 +1000 |
commit | ed98174bbefc85f9895c870fdb7518e2a534a150 (patch) | |
tree | 913c1eb5cb384213d22af6cf4f4e2b8fd2a51276 /test-tnpair.nix | |
parent | 134f5e0bfcd5f22e3adee4b707e10693cf71e635 (diff) |
bugfix: use unsorted type for alignment output
Most aligners are only pseudo-name sorted (i.e., pairs adjacant).
Diffstat (limited to 'test-tnpair.nix')
-rw-r--r-- | test-tnpair.nix | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/test-tnpair.nix b/test-tnpair.nix index 7cc5453..6eb2e23 100644 --- a/test-tnpair.nix +++ b/test-tnpair.nix @@ -11,7 +11,7 @@ let ref = fetchfa ./examples/ref.fa; - alignWithRG = rg: bwa.align { inherit ref; flags = "-R'@RG\\tID:${rg}\\tSM:${rg}'";}; + alignWithRG = rg: x: nameSort (bwa.align { inherit ref; flags = "-R'@RG\\tID:${rg}\\tSM:${rg}'";} x); sort = samtools.sort {}; nameSort = samtools.sort {nameSort = true;}; flagstat = samtools.flagstat {}; @@ -41,7 +41,7 @@ let glvariants = strelka.call {} (builtins.attrValues alignments); platypusVars = platypus.call {} (builtins.attrValues alignments); octopusVars = octopus.call {} (builtins.attrValues alignments); - shards = map (bwa.align {inherit ref;}) (shard.fastQPair 2 normal.files); + shards = map (x: nameSort (bwa.align {inherit ref;} x)) (shard.fastQPair 2 normal.files); }; tnpairResult = processPair tnpair; |