How to align genomic sequence with corresponding amino acid sequence

Question

Does anyone know of a program that can align a genomic sequence with introns with the corresponding amino acid sequence?

I have both the genomic sequence and the correct amino acid sequence but no information on the genemodel in e.g. gff or genbank format.

This is not what one would call a 'normal' alignment, but I would like it to look like this:

Hypothetical example:

nucleotide sequence with intron in italic:

ATGCACGATACGACTGACGTACGTACGTACGTACGTACGTACGTACGGAGACGTAGACTC

corresponding amino acid sequence:

MHDTTDVRTYVRRRRL

The intron does not encode for amino acid which should thus result in a gap in the 'alignment'. Desired 'alignment':

ATGCACGATACGACTGACGTACGTACGTACGTACGTACGTACGTACGGAGACGTAGACTC
  M  H  D  T  T  D  V  R              T  Y  V  R  R  R  R L

Answer 1

Have a look at exonerate's protein2genome! From the documentation:

Aligning a protein to genomic sequence:

Similarly, it is possible to align a protein sequence to the genome, (similar to GeneWise, but with heuristics).

exonerate --model protein2genome query.fasta target.fasta

This model will allow introns in the alignment, but also allow frameshifts, and exon phase changes when a codon is split by an intron.

Answer 2

The classic tools for this are genewise and the already mentioned exonerate. Both are capable fo aligning protein sequences to genomic DNA (or cDNA) while modelling splice sites to give you a gene model.

For example, this is the output¹ of genewise on your sequences:

EMBOSS_001         1 MHDTTDVRTYVR----RRRL                              
                     MHDTTDVRTYVR    RRRL                              
                     MHDTTDVRTYVRTYVRRRRL                              
EMBOSS_001         1 acgaaggcatgcatgcacac                              
                     taaccatgcatgcatggggt                              
                     gctgtcatgcatgcagatac                              

The output also includes a GFF of your sequence:

EMBOSS_001  GeneWise    match   1   60  22.88   +   .   EMBOSS_001-genewise-prediction-1
EMBOSS_001  GeneWise    cds 1   60  0.00    +   0   EMBOSS_001-genewise-prediction-1    

Back when I was using these tools daily to align proteins to genomes, I found that exonerate was much faster but sometimes slightly less accurate. This was a few years ago now though, so that might have changed. I would always suggest trying exonerate first and only going for genewise if it fails.

---

¹_{This link will only be valid for a short while, but I just used the sequences from the OP as input for the GeneWise webserver.}