tools/samtools-queryRegion.nix


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22

{ bionix
, flags ? null
, regions ? []
}:

input:

with bionix;
with lib;
with types;

assert (matchFiletype "samtools-queryRegion" { fa = _: true; } input);

stage {

  name = "samtools-faidx";
  buildInputs = with pkgs; [ samtools ];
  buildCommand = ''
    ln -s ${input} input.fasta
    samtools faidx ${optionalString (flags != null) flags} input.fasta ${concatStringsSep " " regions} > $out
  '';
}