Recycler
Recycler copied to clipboard
Shamir-Lab
No output produced
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.
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!
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?