blob: 62a89f60c4c64815cc85b0992e65957405fbfcbe (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
|
{ bionix
, bwaIndexAttrs ? { }
, faidxAttrs ? { }
, indexAttrs ? { }
, assemblyAttrs ? { }
, collectMetricsAttrs ? { }
, softClipsToSplitReadsAttrs ? { flags = "REALIGN_ENTIRE_READ=true"; }
, config ? null
, heapSize ? "4g"
, shards ? 10
}:
with bionix;
with lib;
with types;
with gridss;
inputs:
let
getref = matchFiletype "gridss-variants" { bam = x: x.ref; };
ref = getref (head inputs);
sorted = matchFileSorting "gridss-variants" { coord = _: true; };
homoRef = length (unique (map getref inputs)) == 1;
linkInput = f: attrs: input: ''
BASENAME=$(basename ${input})
WRKDIR="''${BASENAME}.gridss.working"
if [[ ! -e $WRKDIR ]] ; then
mkdir $WRKDIR
fi
for f in ${f attrs input}/* ; do
ln -s $f $WRKDIR/$BASENAME.''${f##*.}
done
'';
linkSV = input: ''
BASENAME=$(basename ${input})
WRKDIR="''${BASENAME}.gridss.working"
if [[ ! -e $WRKDIR ]] ; then
mkdir $WRKDIR
fi
ln -s ${input} $WRKDIR/$BASENAME.sv.bam
ln -s ${bionix.samtools.index indexAttrs input} $WRKDIR/$BASENAME.sv.bai
'';
assembly = bionix.samtools.sort { }
(softClipsToSplitReads softClipsToSplitReadsAttrs
(bionix.gridss.shardedAssemble shards assemblyAttrs inputs));
mkLinks = ''
ln -s ${ref} ref.fa
ln -s ${bionix.samtools.faidx faidxAttrs ref} ref.fa.fai
for f in ${bionix.bwa.index bwaIndexAttrs ref}/*; do
ln -s $f
done
${concatMapStringsSep "\n" linkSV inputs}
${linkSV assembly}
${concatMapStringsSep "\n" (linkInput collectMetrics collectMetricsAttrs)
inputs}
${linkInput collectMetrics collectMetricsAttrs assembly}
'';
in
assert (all sorted inputs);
assert homoRef;
rec {
identify = stage rec {
name = "gridss-identifyVariants";
buildInputs = with pkgs; [ jre samtools ];
buildCommand = mkLinks + ''
java -Xmx${heapSize} -Dsamjdk.create_index=true \
-cp ${jar} gridss.IdentifyVariants \
VERBOSITY=WARNING \
REFERENCE_SEQUENCE=ref.fa \
${concatMapStringsSep " " (i: "INPUT='${i}'") inputs} \
ASSEMBLY=${assembly} \
WORKER_THREADS=$NIX_BUILD_CORES \
OUTPUT_VCF=out.vcf \
${
optionalString (config != null)
("OPTIONS_FILE=" + bionix.gridss.gridssConfig config)
} \
WORKING_DIR=$TMPDIR/ \
TMP_DIR=$TMPDIR/
mv out.vcf $out
'';
passthru = {
filetype = filetype.vcf { inherit ref; };
gridss.assembly = assembly;
multicore = true;
};
};
annotate = stage rec {
name = "gridss-annotateVariants";
buildInputs = with pkgs; [ jre ];
buildCommand = mkLinks + ''
ln -s ${
bionix.gridss.identifyVariants {
inherit bwaIndexAttrs faidxAttrs indexAttrs assemblyAttrs
collectMetricsAttrs softClipsToSplitReadsAttrs config;
} inputs
} input.vcf
java -Xmx${heapSize} -Dsamjdk.create_index=true \
-cp ${jar} gridss.AnnotateVariants \
VERBOSITY=WARNING \
REFERENCE_SEQUENCE=ref.fa \
${concatMapStringsSep " " (i: "INPUT='${i}'") inputs} \
ASSEMBLY=${assembly} \
INPUT_VCF=input.vcf \
OUTPUT_VCF=out.vcf \
WORKER_THREADS=$NIX_BUILD_CORES \
WORKING_DIR=$TMPDIR/ \
${
optionalString (config != null)
("OPTIONS_FILE=" + bionix.gridss.gridssConfig config)
} \
TMP_DIR=$TMPDIR/
mv out.vcf $out
'';
passthru = {
filetype = filetype.vcf { inherit ref; };
gridss.assembly = assembly;
multicore = true;
};
};
annotateAndAssemble = stage rec {
name = "gridss-annotateVariants-assemble";
buildCommand = ''
mkdir $out
ln -s ${annotate} $out/gridss.vcf
ln -s ${assembly} $out/gridss.bam
ln -s ${bionix.samtools.index indexAttrs assembly} $out/gridss.bam.bai
'';
};
}
|
