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

overlapping prefix adapters #53

Closed
marcelm opened this Issue Mar 13, 2015 · 3 comments

Comments

Projects
None yet
1 participant
@marcelm
Owner

marcelm commented Mar 13, 2015

From marcel.m...@tu-dortmund.de on May 08, 2012 17:27:10

The new -g ^ADAPTER option isn’t enough. There has been a request to allow less strict anchoring, where the adapter overlaps the beginning of the read.

This is easily achieved by this change:
-PREFIX = align.STOP_WITHIN_SEQ2
+PREFIX = align.STOP_WITHIN_SEQ2 | align.START_WITHIN_SEQ1

The question is whether that is the desired behaviour or whether both versions should be possible.

Original issue: http://code.google.com/p/cutadapt/issues/detail?id=43

@marcelm

This comment has been minimized.

Owner

marcelm commented Mar 13, 2015

From seb.th...@gmail.com on November 07, 2014 08:01:56

Yes this behaviour could very usefull. I have tons of examples were my MID has been "eaten" at the beginning of the read. Combining -e and -O allows substitutions and that's too permissive.

And by the way, thanks for the trick.

@marcelm

This comment has been minimized.

Owner

marcelm commented Mar 13, 2015

From marcel.m...@tu-dortmund.de on November 10, 2014 05:35:23

Thanks for the feedback! I’ll try to come up with a way of allowing this. It may take a while, howere (although I’m aware that this report is already two years old).

@marcelm

This comment has been minimized.

Owner

marcelm commented Jan 27, 2017

This can be achieved with -g XXXXADAPTER. The X characters are interpreted as IUPAC wildcards that do not match any nucleotide. If the number of X is lower than the number of allowed errors, the match will never be found within the read.

@marcelm marcelm closed this in 40907f8 Apr 30, 2018

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment