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

exonerate_hits.py returning reverse-complement sequences #38

Open
mossmatters opened this issue Aug 8, 2018 · 0 comments
Open

exonerate_hits.py returning reverse-complement sequences #38

mossmatters opened this issue Aug 8, 2018 · 0 comments

Comments

@mossmatters
Copy link
Owner

mossmatters commented Aug 8, 2018

A user identified a seemingly random bug where half of the sequences were very different. Upon closer inspection, the sequences were reverse complement.

This is a result of addressing: #27 in January 2018. Sequences recovered before the update to 1.3 should not be affected.

This should not have been fixed. When there is one single hit, HybPiper should return the exonerate results without modification. Strandedness is only needed when building a "supercontig" from multiple exonerate hits.

Users who recovered sequences with HybPiper 1.3 can fix their sequences by updating HybPiper to 1.3.1 (for example with git pull) and re-running just the exonerate step. Add the following flags to read_first.py:

reads_first.py --prefix sampleName --no-blast --no-assemble --no-distribute

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