Skip to content

Commit

Permalink
fix: snap task
Browse files Browse the repository at this point in the history
Fix for wrong input in hmm-assembler.pl step. Also adds the ~{extra_params} so it works in default use case.
  • Loading branch information
MmasterT committed Jan 25, 2024
1 parent cab2268 commit 1654b81
Showing 1 changed file with 6 additions and 4 deletions.
10 changes: 6 additions & 4 deletions annotation/prediction_module/main.wdl
Original file line number Diff line number Diff line change
Expand Up @@ -48,7 +48,7 @@ workflow ei_prediction {

String? codingquarry_extra_params
String? glimmer_extra_params
String? snap_extra_params
String snap_extra_params = " -lcmask "
String? augustus_extra_params
String? evm_extra_params

Expand Down Expand Up @@ -293,6 +293,7 @@ workflow ei_prediction {
genome = def_reference_genome,
transcripts = select_first([def_training_models]),
pretrained_hmm = snap_training,
species = species,
extra_params = snap_extra_params
}
}
Expand Down Expand Up @@ -1060,7 +1061,8 @@ task SNAP {
IndexedReference genome
File transcripts
File? pretrained_hmm
String? extra_params
String? extra_params,
String species,
RuntimeAttr? resources
}

Expand Down Expand Up @@ -1099,8 +1101,8 @@ task SNAP {
cd params
forge ../export.ann ../export.dna
cd ..
hmm-assembler.pl ~{genome.fasta} params > ~{basename(genome.fasta)}.hmm
snap ~{basename(genome.fasta)}.hmm ~{genome.fasta} > snap.predictions.zff
hmm-assembler.pl ~{species} params > ~{basename(genome.fasta)}.hmm
snap ~{extra_params} ~{basename(genome.fasta)}.hmm ~{genome.fasta} > snap.predictions.zff
else
snap ~{extra_params} ~{pretrained_hmm} ~{genome.fasta} > snap.predictions.zff
fi
Expand Down

0 comments on commit 1654b81

Please sign in to comment.