This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Genome Annotation - Finding Genes Program

Hi,

May be my question is very simple, I want to find some genes in a eukaryotic genome, I have my genes (mRNA - sheep, cattle, goat, human) and a genome (alpaca) in multifasta, I tried with homology-based programs, AAT and GeneSeqer, but they don't give me a GFF output, ATT gave me a alignment and a statistics files for each genomic sequence, so It gave me to many files, and GeneSeqer gave me a statistics and alignment in one file, I would to have only one statistics file, with the start end position, and a GFF output. Any ideas?

annotation gene
• 5,285 views • link
updated by Biostar Biostar • written by Daniel Fernández • Daniel Fernández

2 answers

Use MAKER.

• 0 views • link
updated by Ram Ram • written by Yannick Wurm Yannick Wurm

Well I was trying to install the MAKER this days, the program is very complete it use many prograsm and the Exonarate too, but for me is a little difficult to install it, this need some prerequisites, DBI, DBD::SQLite, Proc::ProcessTable, threads, IO::All and IO::Prompt; I installed all except DBD::SQLite and IO::All, in its README file did not say how to, maybe is like the others, (perl Makefile.PL; make; make test; make install) but I don't sure, I don't want to spoil the installation, can you tell how I can installed that, please?.

• 0 views • link
updated by Ram Ram • written by Daniel Fernández • Daniel Fernández

Exonerate will do the job. It can align genes, mRNA, proteins, EST, or genome to genome.

It outputs GFF (--showtargetgff yes).

The set of parameters I use (fungal genomes):

--showvulgar yes
--softmaskquery yes
--softmasktarget yes
--minintron 20
--maxintron 3000
--ryo ">%qi length=%ql alnlen=%qaln>%ti length=%tl alnlen=%taln"
--showalignment no
--showtargetgff yes
--bestn 5
--percent 70
• 1 views • link
updated by Ram Ram • written by Leszek Leszek

Log in to answer this question.

More posts like this