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 polished contigs despite ~50X coverage #207

Open
gettl008 opened this issue Jan 31, 2022 · 5 comments
Open

No polished contigs despite ~50X coverage #207

gettl008 opened this issue Jan 31, 2022 · 5 comments

Comments

@gettl008
Copy link

I've used Racon on > 200 assemblies and have only had this error arise with one sample...

The reads in this case are ONT reads that I mapped to my assembly using Minimap2. The coverage distribution (removing the high coverage tail) looks like this...
image

When I run Racon, everything looks normal, there are no errors, and it ends with
[racon::Polisher::] total = 24.029181 s
The output file, however, is empty. When I run with '-u', the output file gets properly filled with all the unpolished contigs.

I'm thinking this has something to do with the input reads as I have two different assemblies I've made for the same sample and both are having the same issue. Any ideas?

@rvaser
Copy link
Collaborator

rvaser commented Jan 31, 2022

Hello,
can you please paste the whole racon log?

Best regards,
Robert

@gettl008
Copy link
Author

gettl008 commented Jan 31, 2022

Of course...

[racon::Polisher::initialize] loaded target sequences 0.072242 s
[racon::Polisher::initialize] loaded sequences 6.214936 s
[racon::Polisher::initialize] loaded overlaps 5.285117 s
[racon::Polisher::initialize] aligning overlaps [====================] 4.301526 s
[racon::Polisher::initialize] transformed data into windows 0.978941 s
[racon::Polisher::polish] generating consensus [====================] 0.563739 s
[racon::Polisher::] total = 17.828767 s

@rvaser
Copy link
Collaborator

rvaser commented Jan 31, 2022

Do you have a fastQ file? If so, what is the average quality score?

@gettl008
Copy link
Author

gettl008 commented Jan 31, 2022

This is a published read set I pulled from NCBI and appears to be missing all quality info. No idea why...

@rvaser
Copy link
Collaborator

rvaser commented Feb 1, 2022

Could you provide the data for this example?

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