Skip to content

Commit 21152ea

Browse files
committed
keeping strand search mmseqs/easysearch
1 parent 9e3c7cf commit 21152ea

File tree

1 file changed

+6
-1
lines changed

1 file changed

+6
-1
lines changed

conf/modules.config

+6-1
Original file line numberDiff line numberDiff line change
@@ -12,6 +12,7 @@
1212

1313
variant_caller = params.variant_caller
1414
intermediate_variant_caller = params.intermediate_variant_caller
15+
search_mode = params.mmseqs_searchtype
1516
assemblers = params.assemblers ? params.assemblers.split(',').collect{ it.trim().toLowerCase() } : []
1617
id_threshold_cdhit = Math.max(params.identity_threshold > 1 ? params.identity_threshold * 0.01 : params.identity_threshold, 0.80)
1718
save_final_reads = params.save_final_reads ? (!params.skip_hostremoval ? 'host' : (!params.skip_complexity_filtering ? 'complexity' : (!params.skip_humid ? 'deduplication' : 'trimming'))) : 'nothing'
@@ -1652,7 +1653,11 @@ process {
16521653
}
16531654

16541655
withName: MMSEQS_EASYSEARCH {
1655-
ext.args= { "--search-type ${params.mmseqs_searchtype} --rescore-mode 3 --format-output query,target,theader,fident,qlen,tlen,alnlen,mismatch,gapopen,qstart,qend,tstart,tend,evalue,bits" } // don't change format-output
1656+
ext.args= [
1657+
mmseqs_searchtype == "3" ? "--search-type ${params.mmseqs_searchtype} --strand 2" : "--search-type ${params.mmseqs_searchtype}",
1658+
"--rescore-mode 3",
1659+
"--format-output query,target,theader,fident,qlen,tlen,alnlen,mismatch,gapopen,qstart,qend,tstart,tend,evalue,bits" // don't change format-output
1660+
]
16561661
publishDir = [
16571662
path: { "${params.outdir}/consensus/quality_control/mmseqs-search"},
16581663
mode: params.publish_dir_mode,

0 commit comments

Comments
 (0)