2022-05-23 08:01:55 +00: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 18:16:45 +00:00
|
|
|
idcutoff = 0.985
|
2022-05-23 12:02:10 +00:00
|
|
|
outoption = "xcfert" // Nonsense text to check default case.
|
|
|
|
columns = ""
|
2022-05-30 11:01:06 +00:00
|
|
|
VSEARCH_USEARCHGLOBAL ( [[id:'test'], query], db, idcutoff, outoption, columns )
|
2022-05-23 12:02:10 +00: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 18:16:45 +00:00
|
|
|
idcutoff = 0.985
|
2022-05-23 12:02:10 +00:00
|
|
|
outoption = "userout"
|
|
|
|
columns = "query+target+id"
|
2022-05-30 11:01:06 +00:00
|
|
|
VSEARCH_USEARCHGLOBAL ( [[id:'test'], query], db, idcutoff, outoption, columns )
|
2022-05-23 08:01:55 +00:00
|
|
|
}
|