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

No output produced #41

Open
sum732 opened this issue Dec 25, 2019 · 3 comments
Open

No output produced #41

sum732 opened this issue Dec 25, 2019 · 3 comments

Comments

@sum732
Copy link

sum732 commented Dec 25, 2019

I have got the recent pull of Recycler installed. The contigs/scaffolds were produced by SPAdes-3.11.1. The kmer size etc were all under 55 etc. Recycler is NOT generating any outputs. The BAM file was produced as per the tutorial.

Any ideas what could be the issues? I am have got about 7 contigs from the assembly ( SPAdes) >500bp ( largest is about 1500bp). Simply mapping these contigs to a reference via BLAST, I have got a really good score and %identity. Also the edges do overlaps. However, no output from Recycler.

@dpellow
Copy link
Collaborator

dpellow commented Dec 26, 2019

Thanks for your interest in using Recycler.
Does the Recycler run complete without producing output or is it getting stuck during the run? If there are many short contigs besides the 7 long ones you mentioned, Recycler could just be taking a long time to complete running.
Recycler will only assemble plasmids from overlapping contigs that form a circular path and have mostly uniform coverage. If this is not the case in your data then Recycler will not assemble plasmids.
Otherwise, we can help you to debug if you can send us more details of the dataset.
Hope that helps!

@sum732
Copy link
Author

sum732 commented Dec 26, 2019

Thank you for the quick response.
Recycler does complete the run and finishes with no output. From "SPAdes" there are short contigs as well.
Here is what is shown on STDERR:

(7.5625, 430.5257300917254, 1336.0)                                                             
================== path, coverage levels when added ====================                        
(4936, ' nodes remain in component\n')                                                          
==================final_paths identities after updates: ================  

And then empty output is produced.
For this plasmid in total there are about 70K PE reads. About the uniform coverage, that needs to be checked. Can you recommend parameters that I should try? I have changed the kmer etc with no affects. Can you recommend any other tools that you think may be useful if coverage is the issue?

@dpellow
Copy link
Collaborator

dpellow commented Dec 26, 2019

I can try to help you with this - send me an email at dpellow@post.tau.ac.il and I'll see what I can do.

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

2 participants