Conditional reciprocal blast confusion
0
0
Entering edit mode
6.2 years ago
Biogeek ▴ 470

Dear all,

I have 2 assemblies. A de novo assembly and a draft genome assembly of an organism. My aim is to confirm that specific genes of the de novo do in fact appear in the draft genome. I will then design primers around the best match obtained form the CRB analysis.

I am using crb-blast: https://github.com/cboursnell/crb-blast.

If I use the transcriptome as the query, several of the transcriptome genes will obtain a positive hit with the same gene within the genome (for instance, several of my transcripts provide a best CRBB hit to the same genome gene).

If I use the genome as the query, several of the genome genes will obtain a positive hit with the same transcript within the transcriptome.

Can someone please help me understand CRB. I was under the idea that you would end up with a 1 for 1 match, instead several matches are coming up for genes/ transcripts which are queried. I am aware that several transcripts are produced by a gene; however, is the idea of CRB not to ultimately provide a 1 for 1 best match?

If this is the case, how do you choose the best sequence to design the primers around? Do i keep the de novo as the query and the target as the genome? Thanks.

crb-blast • 1.2k views
ADD COMMENT

Login before adding your answer.

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