aboutsummaryrefslogtreecommitdiff
path: root/test-tnpair.nix
diff options
context:
space:
mode:
authorJustin Bedo <cu@cua0.org>2021-07-28 08:07:32 +1000
committerJustin Bedo <cu@cua0.org>2021-07-28 09:17:56 +1000
commited98174bbefc85f9895c870fdb7518e2a534a150 (patch)
tree913c1eb5cb384213d22af6cf4f4e2b8fd2a51276 /test-tnpair.nix
parent134f5e0bfcd5f22e3adee4b707e10693cf71e635 (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.nix4
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;