这是indexloc提供的服务,不要输入任何密码
Skip to content
Draft
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
55 changes: 55 additions & 0 deletions sniffles2.wdl
Original file line number Diff line number Diff line change
@@ -0,0 +1,55 @@
version 1.0

task Sniffles {
input {
File bam
File bamIndex
File? referenceFasta
File? tandemRepeats
File? genotypeVcfs

String outputVcf = "output.vcf"
String options = " "

Boolean mosaic = false
Boolean outputRnames = false

Int threads = 4
String memory = "32GiB" # 16 was fine without --no-qc
Int timeMinutes = 1 + ceil(size(bam, "GiB") * 2)
String dockerImage = "quay.io/biocontainers/sniffles:2.6.2--pyhdfd78af_0"
}

meta {
citation: {
doi: "10.1038/s41587-023-02024-y"
}
}

command {
set -e
mkdir -p "$(dirname ~{outputVcf})"

sniffles \
~{"--reference " + referenceFasta} \
~{"--tandem-repeats " + tandemRepeats} \
~{"--genotype-vcf " + genotypeVcfs} \
--threads ~{threads} \
~{true="--output-rnames" false="" outputRnames} \
~{true="--mosaic" false="" mosaic} \
--input ~{bam} \
~{options} \
--vcf ~{outputVcf}
}

output {
File vcf = outputVcf
}

runtime {
cpu: threads
memory: memory
time_minutes: timeMinutes
docker: dockerImage
}
}
Loading