2022-05-23 04:01:55 -04:00
|
|
|
#!/usr/bin/env nextflow
|
|
|
|
|
|
|
|
nextflow.enable.dsl = 2
|
|
|
|
|
|
|
|
include { VSEARCH_USEARCHGLOBAL } from '../../../../modules/vsearch/usearchglobal/main.nf'
|
|
|
|
|
|
|
|
workflow test_vsearch_usearchglobal {
|
|
|
|
|
|
|
|
query = file(params.test_data['sarscov2']['genome']['transcriptome_fasta'], checkIfExists: true)
|
|
|
|
db = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
|
2022-05-26 14:16:45 -04:00
|
|
|
idcutoff = 0.985
|
2022-05-23 08:02:10 -04:00
|
|
|
outoption = "xcfert" // Nonsense text to check default case.
|
|
|
|
columns = ""
|
2022-05-30 07:01:06 -04:00
|
|
|
VSEARCH_USEARCHGLOBAL ( [[id:'test'], query], db, idcutoff, outoption, columns )
|
2022-05-23 08:02:10 -04:00
|
|
|
}
|
|
|
|
|
|
|
|
workflow test_vsearch_usearchglobal_userout {
|
|
|
|
|
|
|
|
query = file(params.test_data['sarscov2']['genome']['transcriptome_fasta'], checkIfExists: true)
|
|
|
|
db = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
|
2022-05-26 14:16:45 -04:00
|
|
|
idcutoff = 0.985
|
2022-05-23 08:02:10 -04:00
|
|
|
outoption = "userout"
|
|
|
|
columns = "query+target+id"
|
2022-05-30 07:01:06 -04:00
|
|
|
VSEARCH_USEARCHGLOBAL ( [[id:'test'], query], db, idcutoff, outoption, columns )
|
2022-05-23 04:01:55 -04:00
|
|
|
}
|