This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Search if a single primer could bind to multiple sequences with Python

Hi all,

I have one primer and a list of sequences. I want to check to which of the sequences (if any) the primer could bind to. It's like an in-silico PCR (so mismatches should be more penalized at the 3' end of the primer than at the 5' one etc.), but with only one primer.

I can't seem to find a tool that accepts only one primer instead of a pair. Anyone has any suggestion? Possibly, I would like to integrate this tool in a Python script.

python primer primer extension alignment
• 3,208 views • link
updated by Camp • Camp • written by gprezza • gprezza

Not a python solution, But you can try seqkit for locating primer sequence in sequence(s) of interest. without penalty, you can use seqkit for locating pcr primer on forward and reverse strand of your sequences:

$ seqkit locate -p <primer_sequence> <target_sequence(s) in fasta format>
• 0 views • link
written by cpad0112 cpad0112

I need a way to check if my primer could be used as such in the target sequence. For this reason, I need a tool that weighs mismatches based on their location in the primer (i.e. almost ignore the ones in the 5' end), exactly like in-silico PCR tools do. If i get it correctly, seqkit doesn't do that.

• 0 views • link
written by gprezza • gprezza

Hey, did u find a solution to your problem? I am in a very similar position. I have to run a pair of primers against a bunch of genomes I've been told to try primersearch on Emboss by python script but im finding it quite complicated

• 0 views • link
written by Camp • Camp

1 answer

You could use PrimerBLAST if you have a second primer that is supposed to work with this primer you have. You can then provide PrimerBLAST with the target sequence plus both fwd and rev primer and a background genome, so PrimerBLAST will check if (against that genome) the primer pair is specific to generate the sequence you want. PrimerBlAST is essentially Primer3 (so its scores and metrics) plus BLAST, so pretty much exactly what you need.

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this