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

Doesn't quite align to end and circ fails #23

Closed
tseemann opened this issue Apr 5, 2018 · 3 comments
Closed

Doesn't quite align to end and circ fails #23

tseemann opened this issue Apr 5, 2018 · 3 comments
Assignees
Labels

Comments

@tseemann
Copy link
Owner

tseemann commented Apr 5, 2018

tig00000001     dna     4736005

 Score = 55308 bits (29950),  Expect = 0.0
 Identities = 29995/30013 (99%), Gaps = 17/30013 (0%)
 Strand=Plus/Plus

Query  1        CGCTGTCGGCAAGAATATAGCGGCTTGATGCCAAAG-CGCCT-GGTCATTTCGACAAAAA  58
                |||||||||||||||||||||||||||||||||||| ||||| |||||||||||||||||
Sbjct  4704147  CGCTGTCGGCAAGAATATAGCGGCTTGATGCCAAAGGCGCCTGGGTCATTTCGACAAAAA  4704206

<snip>

Query  29988    ACGGTTTTTCAGT  30000
                |||||||||||||
Sbjct  4734143  ACGGTTTTTCAGT  4734155

@tseemann tseemann added the bug label Apr 5, 2018
@tseemann tseemann self-assigned this Apr 5, 2018
@tseemann
Copy link
Owner Author

tseemann commented Apr 5, 2018

I think this is just because the default --readlen if 30,000 is shorter than the overhang!

@tseemann
Copy link
Owner Author

tseemann commented Apr 5, 2018

Still can happen:

tig00000001     dna     4736706

Query  31657    CAATGAAACCGCAGCTGCATTAAATAAAAATACCAACTCATCCTTTCTTTCCGTAATCAA  31716
                ||||||||||||||| ||||||||||||||||||||||||||| |||||| |||||||||
Sbjct  4736624  CAATGAAACCGCAGC-GCATTAAATAAAAATACCAACTCATCC-TTCTTTTCGTAATCAA  4736681

Query  31717    AAATCCAGTTTT-GCTTCAGCTTT  31739
                |||||||||||| ||||||| |||
Sbjct  4736682  AAATCCAGTTTTTGCTTCAGATTT  4736705

This one was only 1 bp short of the end. Need to apply the same rule to the "1" starting position as we do to the end position.

@tseemann
Copy link
Owner Author

tseemann commented Apr 5, 2018

Added new --fuzz parameter.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

1 participant