ksahlin / NGSpeciesID

Reference-free clustering and consensus forming of long-read amplicon sequencing
GNU General Public License v3.0
49 stars 14 forks source link

Allow IUPAC code in primer sequences #7

Closed jdalino closed 3 years ago

jdalino commented 3 years ago

Add IUPAC equivalence map in find_barcode_locations so edlib can understand IUPAC code as correct match for primer research.

jdalino commented 3 years ago

Hello @ksahlin, I worked with NGSpeciesID recently (very nice tool btw) and found that IUPAC code in primer sequences are not recognized. I added a equivalence map in the find_barcode_locations for edlib to consider them as match. Feel free to edit or refuse if you don't want this to be the standard behaviour.

ksahlin commented 3 years ago

Super thanks for this pull request, very useful! I wanted to wait with the merge until I had fixed the travis automatic build, which is now done.