Skip to content
New issue

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

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

Already on GitHub? Sign in to your account

Is metaviral mode's Length threshold flexible? #1437

Open
handibles opened this issue Jan 13, 2025 · 0 comments
Open

Is metaviral mode's Length threshold flexible? #1437

handibles opened this issue Jan 13, 2025 · 0 comments

Comments

@handibles
Copy link

Is your feature request related to a problem? Please describe. For generic questions use Q&A section in the Discussions forum above.

Thanks for all the work devs.

I'm comparing metaviral SPAdes with standard SPAdes+virsorter2. Usually we end up placing a 10Kbp threshold on sequence length post-hoc, but I'm wondering how increasing the Length threshold (as mentioned in the paper) to a higher value would affect the graph assembly? From the paper:

We classify such a path as long if its length exceeds a threshold Length (the default value 1000 nucleotides) and high coverage if its coverage exceeds a threshold Coverage (the default value 5×). Long high-coverage paths represent putative sequences of linear viruses.

I've found similar looking thresholds in the hybrid_aligner defn in config.info, but that might be for hybrid assembly specifically:

long_read_threshold 1000

Describe the solution you'd like

  • should the length parameter be considered flexible? Can it be referenced from the CLI?
  • Would modifying the value above change the programme as intended, or are those parameters for long/short hybrid assembly? If not, is there another location this can be increased from 1,000 to (e.g.) 10,000?

Describe alternatives you've considered

Currently, holding off on brutalising the lovely code.

Additional context

Thanks for the help.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant