Skip to content

Commit

Permalink
modify strand setting to match with immuno.wdl. Note that the strand …
Browse files Browse the repository at this point in the history
…input will be first, second, unstranded, not RF, FR, XS.
  • Loading branch information
mhoang22 committed Jan 22, 2025
1 parent 7bfc583 commit 4ea6617
Showing 1 changed file with 8 additions and 2 deletions.
10 changes: 8 additions & 2 deletions definitions/tools/regtools.wdl
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ task regtools {
String? output_filename_vcf = "splice_variant.vcf"
String? output_filename_bed = "splice_variant.bed"

String strand # enum [XS, RF, FR]
String strand = "unstranded" # [first, second, unstranded]
Int? window_size
Int? max_distance_exon # max distance from exon/intron boundary to annotate a variant in exonic region as splicing variant
Int? max_distance_intron
Expand Down Expand Up @@ -34,12 +34,18 @@ task regtools {
disks: "local-disk ~{space_needed_gb} HDD"
}

Map[String, String] strandness = {
"first": "RF",
"second": "FR",
"unstranded": "XS"
}

command <<<
/regtools/build/regtools cis-splice-effects identify \
-o ~{output_filename_tsv} \
~{if defined(output_filename_vcf) then "-v ~{output_filename_vcf}" else ""} \
~{if defined(output_filename_bed) then "-j ~{output_filename_bed}" else ""} \
-s ~{strand} \
-s ~{strandness[strand]} \
~{if defined(window_size) then "-w ~{window_size}" else ""} \
~{if defined(max_distance_exon) then "-e ~{max_distance_exon}" else ""} \
~{if defined(max_distance_intron) then "-i ~{max_distance_intron}" else ""} \
Expand Down

0 comments on commit 4ea6617

Please sign in to comment.