Getting Exact Number Of Genomic Copies For A Specific Sequence
2
3
Entering edit mode
13.2 years ago
Fucitol ▴ 140

Hi all,

I'm trying to retrieve a set of sequences that can be used for primer design. These sequences have some requirements, e.g. present in 4 different species, not present in a number of species on a blacklist, at least 200bp in length and preferably 50 or more copies of the sequence in all species.

So far I've succeeded and have a total of about 10 sequences that meet these requirements, however I'm unsure if I'm correct about the number of copies for each of the sequences.

For calculating the number of copies, I'm using plain BLAST with 95% identity and count the number of hits for all sequences (where alignment length >= 200bp and non-overlapping start/end positions in the target genome).

Is this a correct method of assessing the number of copies throughout a genome? What other methods are available that I can use for validation?

Thanks for any tips or solutions!

Edit: Thanks for the responses so far, however my goal here is not to design the primers myself, that will be taken care of by someone else. The primer designer gave those requirements as that is what their software requires. The eventual primers will be LAMP primers (process described at Eiken GENOME SITE. So for now I'm only interested in counting the number of copies for each of my candidate sequences.

blast primer • 2.9k views
ADD COMMENT
0
Entering edit mode
ADD REPLY
1
Entering edit mode
13.2 years ago
Rob Syme ▴ 540

A good free solution to this problem is e-PCR and re-PCR, by NCBI.

You can download the software directly from the NCBI ftp site.

These will pick up sites in your genomic fasta sequence where primers will bind to produce a product.

ADD COMMENT
0
Entering edit mode
13.2 years ago
Fabian Bull ★ 1.3k

You may have a look at Tallymer. It gives you the occurrence of any k-mer in a given sequence dataset.

ADD COMMENT

Login before adding your answer.

Traffic: 2183 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6