-
Notifications
You must be signed in to change notification settings - Fork 16
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
An error when looking for 200 proteins in a genome #46
Comments
Could you share me with the proteins and the reference genome? |
Yes, thanks. |
I downloaded the genome and could get the results:
Based on the total number of bases, we are using the same reference. What version are you using? |
How do I know which version I have? |
|
0.11-r235-dirty |
Hi,
I get the following error when running miniprot looking for 200 proteins at a 73MB genome:
[M::mp_ntseq_read@0.3751.01] read 68163509 bases in 69053 contigs
[M::mp_idx_build@0.3791.01] 608218 blocks
[M::mp_idx_build@0.9872.83] collected syncmers
[M::mp_idx_build@13.9181.13] 23643989 kmer-block pairs
[M::mp_idx_print_stat] 1694743 distinct k-mers; mean occ of infrequent k-mers: 13.95; 0 frequent k-mers accounting for 0 occurrences
##gff-version 3
BUG! 1929 == 1929? 621 == 622? 40M1I57M3D83M2D98M3I8M1I59M1D26M15U117M1D46M2D30M5I27M1I5M63U12M
miniprot_: align.c:195: mp_extra_cal: Assertion `al == r->qe - r->qs' failed.
Aborted
I tried to run it looking for 20 proteins in the same genome and it worked fine, and also looked for the 200 proteins in a smaller genome, and it works fine. So I don't think it is a file format issue.
Any suggestions?
Thanks,
Gil
The text was updated successfully, but these errors were encountered: