Replace usearch with an open source tool in AssemblePairs-reference

Issue #45 resolved
Jason Vander Heiden created an issue

vsearch doesn't provide local alignment methods, so we cannot replace usearch with vsearch in AssemblePairs-reference. We need an alternative.

Binaries for BWA and swipe and linux only, so that unfortunately invalidates them. Perhaps the SSW implementation in scikit-bio?

Comments (4)

  1. Jason Vander Heiden reporter

    Added blastn support in 9c90566. usearch and blastn do not yield the same results, apparently because blastn is being more strict returning hits for some unknown reason.

  2. Jason Vander Heiden reporter

    Syncing the word size and masking between ublast and blastn seems to have resolved most of the differences. Still some testing to do, but I'm going to go with the blastn solution for now.

    Dockerfiles have also been updated with blast+ installs.

  3. Jason Vander Heiden reporter

    Seems to work fine as of 5102751. There are minor differences between the usearch and blastn alignments, but this may be impractical to resolve (1 indel problem out of 100 alignments).

  4. Log in to comment