diff options
Diffstat (limited to 'tools/mutect-call.nix')
-rw-r--r-- | tools/mutect-call.nix | 42 |
1 files changed, 42 insertions, 0 deletions
diff --git a/tools/mutect-call.nix b/tools/mutect-call.nix new file mode 100644 index 0000000..9fd009b --- /dev/null +++ b/tools/mutect-call.nix @@ -0,0 +1,42 @@ +{bionix +, nixpkgs +, cosmic +, dbsnp}: + +with nixpkgs; +with lib; + +let + inherit (bionix.types) matchFiletype; + getVCFref = matchFiletype "mutect-call" {vcf = {ref}: ref;}; + getBAMref = matchFiletype "mutect-call" {bam = {ref, ...}: ref;}; + refs = map getVCFref [ cosmic dbsnp ]; + ref = head refs; +in + +assert (length (unique refs) == 1); + +{normal, tumour}: + +assert (ref == getBAMref normal && ref == getBAMref tumour); + +stdenv.mkDerivation { + name = "mutect"; + buildInputs = [ bionix.mutect.app ]; + buildCommand = '' + ln -s ${normal} normal.bam + ln -s ${tumour} tumour.bam + ln -s ${dbsnp} dbsnp.vcf + ln -s ${cosmic} cosmic.vcf + ln -s ${ref} ref.fa + ln -s ${bionix.samtools.faidx {} ref} ref.fa.fai + ln -s ${bionix.samtools.dict {} ref} ref.dict + mutect --analysis_type MuTect \ + --reference_sequence ref.fa \ + --cosmic cosmic.vcf \ + --dbsnp dbsnp.vcf \ + --input_file:normal normal.bam \ + --input_file:tumour tumour.bam \ + --out $out + ''; +} |