Skip to content
New issue

Have a question about this project? # for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “#”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? # to your account

Minimum frequency threshold is being ignored #382

Open
AndreaAguadoM opened this issue Jun 6, 2023 · 2 comments
Open

Minimum frequency threshold is being ignored #382

AndreaAguadoM opened this issue Jun 6, 2023 · 2 comments
Labels
bug Something isn't working

Comments

@AndreaAguadoM
Copy link

Description of the bug

Hello.
I am trying to use viralrecon as to detect variants with a minimun frequency of 0.1. Nevertheless, when using viralrecon with the instructions given below, I have found the following warning:

WARN: Found unexpected parameters:

  • --min_allele_freq: 0.1

the min_allele_freq is being ignored. What can I do?

Thanks in advance!

Command used and terminal output

$ nextflow run nf-core/viralrecon --input samplesheet.csv --genome 'MN908947.3' --min_allele_freq 0.1 -profile singularity --primer_bed nCoV-2019.artic.V4_1.bed --outdir viralrecon_output/ --platform illumina --protocol amplicon

Relevant files

No response

System information

No response

@AndreaAguadoM AndreaAguadoM added the bug Something isn't working label Jun 6, 2023
@maxulysse
Copy link
Member

There is no min_allele_freq params.
Where did you find reference to this params?

@TKMarkCheng
Copy link

I think this is in reference to --min-allele_freq in ivar_variants_to_vcf.py

"--allele_freq_threshold",

Which is a post-processing step downstream of ivar_variants.

nextflow.config itself actually doesn't take such an argument. In modules_illumina.config, the threshold is fixed at 0.25, which is likely the issue.

withName: 'IVAR_VARIANTS' {
ext.args = '-t 0.25 -q 20 -m 10'

# for free to join this conversation on GitHub. Already have an account? # to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

3 participants