Is There A Package To Map Primer Pair On Genome In R/Bioconductor?
2
0
Entering edit mode
11.1 years ago

Hi, everyone!

I have several primer pairs and want to map them onto genome, then calculate the coverage.

My plan now is as follows:

  1. using BSgenome to load the genome library;

  2. using function matchPattern in BSgenome to get the start and end information of the primer pairs;

  3. using function getseq to get sequences;

  4. alignment these sequences with genome sequence and calculate the coverage.

I'm wondering is there a packages which can map the primer pair on genome and get the sequences directly?

Thanks in advance!

bioconductor coverage primer • 3.9k views
ADD COMMENT
1
Entering edit mode

Hmmm. I've got a wrapper for invoking NCBI's e-pcr and importing a table of results...... What your after...?

ADD REPLY
0
Entering edit mode
10.5 years ago
Dan ▴ 520

Check this question: How to find a 28bp 'primer' sequence in a genome?

ADD COMMENT
0
Entering edit mode
10.5 years ago
Ryan Thompson ★ 3.6k

You could probably use Rsubread to map the primer pair as a paired-end dataset with only one pair in it. Make sure to set the appropriate options for multi-mapping and such.

ADD COMMENT

Login before adding your answer.

Traffic: 1272 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