aboutsummaryrefslogtreecommitdiff
path: root/scripts/variant-calling-bowtie-bwa-mem.scm
diff options
context:
space:
mode:
Diffstat (limited to 'scripts/variant-calling-bowtie-bwa-mem.scm')
-rw-r--r--scripts/variant-calling-bowtie-bwa-mem.scm14
1 files changed, 5 insertions, 9 deletions
diff --git a/scripts/variant-calling-bowtie-bwa-mem.scm b/scripts/variant-calling-bowtie-bwa-mem.scm
index 59c1337..598754b 100644
--- a/scripts/variant-calling-bowtie-bwa-mem.scm
+++ b/scripts/variant-calling-bowtie-bwa-mem.scm
@@ -3,25 +3,22 @@
(use-modules (ccwl ccwl))
(define fastq-1
- (input "fastq_1" #:type 'File))
+ (input "fastq_1"))
(define fastq-2
- (input "fastq_2" #:type 'File))
+ (input "fastq_2"))
(define bowtie2-index
(input "bowtie2_index"
- #:type 'File
#:other '((secondary-files . #(".1.bt2" ".2.bt2" ".3.bt2" ".4.bt2"
".rev.1.bt2" ".rev.2.bt2")))))
(define bwa-mem-index
(input "bwa_mem_index"
- #:type 'File
#:other '((secondary-files . #(".amb" ".ann" ".bwt" ".pac" ".sa")))))
(define reference
(input "reference"
- #:type 'File
#:other '((secondary-files . #(".fai")))))
(define (bam2vcf prefix)
@@ -29,7 +26,7 @@
(list (pipeline "samtools"
(list (command "samtools_view"
(list "samtools" "view" "-h" "-b" "-F" "4"
- (input (string-append prefix "_input") #:type 'File)))
+ (input (string-append prefix "_input"))))
(command "samtools_sort"
(list "samtools" "sort" "-n"))
(command "samtools_fixmate"
@@ -44,7 +41,7 @@
#:source "samtools_markdup/stdout")))
(command "samtools_index"
(list "samtools" "index"
- (input "dedup_bam" #:type 'File))
+ (input "dedup_bam"))
#:outputs (list (output "indexed_dedup_bam"
#:type 'File
#:binding '((glob . "$(inputs.dedup_bam.basename)"))
@@ -55,7 +52,6 @@
(list (command "freebayes"
(list "freebayes" "--ploidy" "1"
"--bam" (input "indexed_dedup_bam"
- #:type 'File
#:other '((secondary-files . #("bai"))))
"-f" reference))
(command "bcftools_view"
@@ -67,8 +63,8 @@
#:type 'File
#:source "bcftools_view/stdout")))
(command "tabix"
- (list "tabix" (input "vcf" #:type 'File))
#:outputs (list (output "indexed_vcf" #:type 'File
+ (list "tabix" (input "vcf"))
#:binding '((glob . "$(inputs.vcf.basename)"))
#:other '((secondary-files . #(".tbi")))))
#:other `((hints (Initial-work-dir-requirement