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

Implement solution for managing low-coverage samples #59

Open
georgiesamaha opened this issue Oct 10, 2024 · 0 comments
Open

Implement solution for managing low-coverage samples #59

georgiesamaha opened this issue Oct 10, 2024 · 0 comments
Labels
enhancement New feature or request

Comments

@georgiesamaha
Copy link
Member

Is your feature request related to a problem? Please describe.
When working with low-coverage and rare samples, a user will rerun a barcode in the same flowcell to increase the number of reads. Fastq files will then be merged into the same barcode directory and run through the pipeline. It is challenging to know if barcodes have enough reads to complete assemblies. Sometimes low coverage samples will cause Unicycler to throw an error re: miniasm.

Describe the solution you'd like
We should be doing some basic evaluation of read depth for barcodes at the input check stage of the workflow, these could throw a warning for samples that may not pass a specific threshold (don't know what this would be). Additionally, we could implement a no-miniasm implementation of unicycler for samples that may not meet this threshold to give it the best chance of success. Things to consider:

  • What would be a reasonable threshold for read depth?
  • What other qc checks should we do, as read depth is not the only important metric for genome assembly?
  • What impact does skipping miniasm have on the assembly and does this matter in our application?
  • How can we implement this in a way that is relatively straight forward and easy to understand for an end user?

Describe alternatives you've considered
None really, except skipping and recalling samples with low coverage.

# for free to join this conversation on GitHub. Already have an account? # to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

1 participant