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

Circlator with short reads? #88

Closed
jcthrash opened this issue Apr 11, 2017 · 2 comments
Closed

Circlator with short reads? #88

jcthrash opened this issue Apr 11, 2017 · 2 comments

Comments

@jcthrash
Copy link

Hello,
I'm trying to verify a circularized genome. I have a single contig with what looks like overlapping ends, and I'd like to try and use circlator to finish the process, but all options appear to require long reads. I've tried doing the bwa mapping manually, but when I try to move on to the other steps, I get a series of errors. Before I go any further, can you tell me if there are any tricks to get circlator to run with short reads?

Thank you.

@martinghunt
Copy link
Collaborator

It was really written for long reads, so all the defaults are assuming pacbio/nanopore read length. I am guessing here so you would have to try it, but something like this:

--b2r_length_cutoff 1000 --b2r_min_read_length 30 --merge_min_length 50 --merge_min_length_merge 50 --merge_ref_end 100 --merge_reassemble_end 100

... and there may be other options that need tweaking. Essentially, make all the numbers a lot smaller!

@jcthrash
Copy link
Author

Ok, thank you!

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