Gene Sequencing Alignment Problem

2 posts / 0 new
Last post
anonya's picture
Gene Sequencing Alignment Problem


I'm trying to align an 800-amino-acid sequence coded by a gene from a
known genome to a large set of ordered contig sequences from a new
genome of a related strain whose sequence is not yet fully resolved.

Using tblastn, I've so far been unsuccessful in finding
positions/sequence sections that align well to the sequence
from the gene of interest and also do not align better to another
sequence from another gene.

Effectively, sequences that have the best alignment to a gene get
'taken' by that gene, leaving the gene of interest with no significant
alignments. Thus, I’m unable to map the position/sequence sections of
this gene among the contigs.

Is there a way I can decisively check whether or not this gene exists
in its entirety, or partially, in the contig set
of the newly sequenced genome, and if yes, determine its location?

Or is there something inherently flawed about my method/approach that
I could please get a recommendation as to how to correct?

Nikail Collins
Nikail Collins's picture

Perhaps, try the nucleitide level first. also, make sure sequence is going in the correct direction.  If you know the ORF sequence, you could search for that within your new sequence to make sure you are on the right track to start.  Does this answer your question?