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

Dealing with primer modifications #8

Open
sscoster opened this issue Dec 5, 2022 · 3 comments
Open

Dealing with primer modifications #8

sscoster opened this issue Dec 5, 2022 · 3 comments

Comments

@sscoster
Copy link

sscoster commented Dec 5, 2022

Hi,

I'd love to use minibar to help me demultiplex custom barcodes and primers with ONT sequencing. I saw a note on the documentation that IUPAC codes Y and R are recognized, but I also have W. Is that recognized?

And then in addition, my primer code has internal deoxyInosine and I'm trying to figure out how to deal with this.

Do you have suggestions for me to try for these issues?

Many Thanks!

@jbh-cas
Copy link
Collaborator

jbh-cas commented Dec 6, 2022

Hello,

The other authors of our paper specified the YR usage for minibar but I see in papers such as "Molecular Tools to Identify and Characterize Malignant Catarrhal Fever Viruses (MCFV) of Ruminants and Captive Artiodactyla" the use of W M and K as well as deoxyInosine which I believe is best matched as N, i.e, any of ACGT.

We can look at adding those to minibar's equality map. Are you familiar with any other primer extensions that find common usage.

best,

jh

@sscoster
Copy link
Author

sscoster commented Dec 7, 2022 via email

@jbh-cas
Copy link
Collaborator

jbh-cas commented Dec 11, 2022

Version 0.24 adds N and the other 2 base equivalent letters available in the primer sequence.

From 0.24 README:

Primer sequences can contain, in addition to the nucleotide identifiers ACGT, the IUPAC codes Y, R, W, S, M, K or N. In the case of Y, a C or T counts as a match. And for R, A or G counts as a match; likewise W: AT, S: CG, M: AC, K: GT, and N any of ACGT. Both primers are assumed to be the same length; consequently, there is one edit distance value used for the primers. Barcodes are also assumed to be the same length.

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

No branches or pull requests

2 participants