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

Building eukaryotic rrna operon database #63

Open
DeniRibicic opened this issue Mar 27, 2022 · 0 comments
Open

Building eukaryotic rrna operon database #63

DeniRibicic opened this issue Mar 27, 2022 · 0 comments

Comments

@DeniRibicic
Copy link

Hi,

I am in a process of trying to build up an eukaryotic full rrna operon database in order to be able to use it with long amplicon ONT generated sequences.
The idea is to use the ncbi reference eukaryotic representative genomes database (ref_euk_rep_genomes) and extract rrna parts using barrnap.
In order to do so, and to retrieve rrna operons of expected size, would the argument --reject do the work here? Say I am expecting size of ca. 4500 bp, and use --reject on 3000 bp, everything below this size length will be filtered out?

Also, I'd appreciate any comments or suggestions on how you would do this, or is this a valid approach at all?

Thanks

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

1 participant