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
Inconsistent gene calls between Pyrodigal-2.0.2 and Prodigal-2.6.3+31b300a #26
Comments
Hi @aaronmussig , thanks for using Pyrodigal! I'll have a look at it. I think the bug is coming from region masking ( |
Okay, I think found the discrepancy between Pyrodigal and Prodigal: Prodigal only masks regions of 50 |
$ python compare.py --prodigal prodigal --genome GCA_900004415.fna -c -m
Genomes closed=True masked=True
Hits genome=GCA_900004415: prodigal=6693, pyrodigal=6693, equal=True Looks like that was it, I'll make a new release. |
Fixed in |
You're quick! Thanks for looking into this |
I'd rather not leave this kind of bugs around before going on holidays, ahah 😉 |
Hello,
I just wanted to say I appreciate the port of Prodigal, it's great work!
There is a case where running GCA_900004415.1 gives different results when running
Pyrodigal-2.0.2
compared toProdigal-2.6.3
andProdigal-2.6.3+31b300a
.Notably, there are quite a few differences where the
gc_cont
value differs slightly, there are also a few cases where genes differ (content and quantity).Below are the commands that I ran (output attached):
The text was updated successfully, but these errors were encountered: