Description
D. melanogaster protein sequences from FlyBase were aligned against each
scaffold in the D. yakuba (DyakCAF1) assembly and the predicted gene models were constructed
using GeMoMa.
Methods
D. melanogaster protein sequences were aligned against the D. yakuba (DyakCAF1) genome
assembly using NCBI TBLASTN with the following parameters:
-evalue 1e-5
-max_intron_length 100000
-matrix BLOSUM80
-gapopen 13
-gapextend 2
-soft_masking true
-db_soft_mask 30
-best_hit_overhang 0.1
-best_hit_score_edge 0.1
The TBLASTN results were used by GeMoMa to produce an initial set of gene predictions.
The GeMoMa predictions are then filtered by the GAF
module in GeMoMa to
produce the final set of gene predictions.
References
Keilwagen J, Wenk M, Erickson JL, Schattat MH, Grau J, Hartung F.
Using intron position conservation for homology-based gene prediction.
Nucleic Acids Res. 2016 May 19;44(9):e89.