/mothur

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Already on GitHub? Sign in to your account

seq.error with align=F #333

Open
m-bakker opened this Issue May 12, 2017 · 0 comments

Comments

Assignees

@pschloss pschloss

@mothur-westcott mothur-westcott

Labels
None yet
Projects
None yet
Milestone
No milestone
3 participants
@m-bakker @pschloss @mothur-westcott
@m-bakker

m-bakker commented May 12, 2017

I'm finding that seq.error with align=F does not find the correct match in a reference database. For instance, it often uses a different reference for a given read than is chosen by classify.seqs against the same database. It is quick to call error=0 for incorrect pairings that generate pairwise alignments like:
xxxxxxxxxxxxxxxxxxxxxxxxxxxyyyyyyy----------------------------
-------------------------------yyyyyyyzzzzzzzzzzzzzzzzzzzzzzzzz

This is with fungal ITS reads.

@pschloss pschloss assigned pschloss and mothur-westcott May 18, 2017

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