diff options
author | Justin Bedo <cu@cua0.org> | 2019-09-25 15:26:34 +1000 |
---|---|---|
committer | Justin Bedo <cu@cua0.org> | 2019-09-25 15:26:34 +1000 |
commit | b2b9a3be4e8544f569cf6155b8a34f405d4a2591 (patch) | |
tree | 24fdb1580832b693fc7078ca78a86ef0360932cf /tools | |
parent | f5161aa32ad0123755e1a33ebec2e54c287f6ac0 (diff) |
facets: assume sorted vcf
Diffstat (limited to 'tools')
-rw-r--r-- | tools/facets-call.nix | 7 |
1 files changed, 1 insertions, 6 deletions
diff --git a/tools/facets-call.nix b/tools/facets-call.nix index 021cacc..a736627 100644 --- a/tools/facets-call.nix +++ b/tools/facets-call.nix @@ -14,11 +14,6 @@ stage { name = "facets"; buildInputs = [ facets.app ]; buildCommand = '' - # Facets requires lexical sorting on the VCF files - grep '^#' ${vcf} > input.vcf - grep -v '^#' ${vcf} | LC_ALL=C sort -t $'\t' -k1,1 -k2,2n >> input.vcf || true - - # Now actually run facets - snp-pileup input.vcf $out ${concatStringsSep " " bams} + snp-pileup ${vcf} $out ${concatStringsSep " " bams} ''; } |