[maker-devel] Alternative start codons

Fields, Christopher J cjfields at illinois.edu
Tue Mar 19 13:04:18 MDT 2013


We had a user notice that MAKER is not observing alternative start codons for bacterial genomes.  For instance, this predicted transcript:

>Xf_Mul_000007-RA transcript Name:"Protein of unknown function" offset:79 AED:0.42 eAED:1.00 QI:79|-1|0|1|-1|1|1|20|24
GTGGGATACAGGCCGCTGATCGCTGATGGCGCGTACCTGAAACTGCTGCTGGACTACTAC
GTTACAGTGCAGCCTTTGCATGCCGATTGGAAAGATCTATATATCATCGCTTGCGCTATT
ACAGCGGCTAAAAAGAGTCTTCAATTTGGCGTAATTCAGTCATTGGCGGGGTAG

Yields this protein sequence.

>Xf_Mul_000007-RA protein AED:0.42 eAED:1.00 QI:79|-1|0|1|-1|1|1|20|24
MPIGKIYISSLALLQRLKRVFNLA

I'm pretty sure I know what is going on, namely that MAKER is treating the 5' end as UTR and looking for the first ATG (there is one in the sequence above).  Is there any way to change this behavior, though?  For instance, allow alternative start codons like GTG/TTG?

chris



More information about the maker-devel mailing list