<div dir="ltr"><div class="markdown-here-wrapper" id="markdown-here-wrapper-205511" style><p style="margin:1.2em 0px!important">Hi, Carson. I’m annotating a 6 Mbp plant mitochondrial genome using GenBank coding nucleotide and protein sequences from related species. I’m seeing 50 genes annotated using MAKER 2.31, and 37 genes annotated using MAKER 2.31.3. The missing genes look good based on the evidence. I see protein_match evidence in the 2.31.3 GFF file, but no resulting gene and mRNA.</p>


<p style="margin:1.2em 0px!important">Is there a ChangeLog indicating the changes from 2.31 to 2.31.3? Do you know of a change that might cause this? What information can I give you that would help debug this? My maker_opts.ctl file follows.</p>


<p style="margin:1.2em 0px!important">Cheers,<br>Shaun</p>
<pre style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;font-size:1em;line-height:1.2em;margin:1.2em 0px"><code class="language-sh" style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;margin:0px 0.15em;padding:0px 0.3em;white-space:pre-wrap;border:1px solid rgb(234,234,234);background-color:rgb(248,248,248);border-top-left-radius:3px;border-top-right-radius:3px;border-bottom-right-radius:3px;border-bottom-left-radius:3px;display:inline;white-space:pre;overflow:auto;border-top-left-radius:3px;border-top-right-radius:3px;border-bottom-right-radius:3px;border-bottom-left-radius:3px;border:1px solid rgb(204,204,204);padding:0.5em 0.7em;display:block!important;display:block;padding:0.5em;color:rgb(51,51,51);background-color:rgb(248,248,255);background-repeat:initial initial">#-----Genome (these are always required)
genome=pg29mt-concat.fa #genome sequence (fasta file or fasta embeded in GFF3 file)
organism_type=eukaryotic #eukaryotic or prokaryotic. Default is eukaryotic

#-----EST Evidence (for best results provide a file for at least one)
est=cds_na.fa #set of ESTs or assembled mRNA-seq in fasta format

#-----Protein Homology Evidence (for best results provide a file for at least one)
protein=cds_aa.fa #protein sequence file in fasta format (i.e. from mutiple oransisms)

#-----Repeat Masking (leave values blank to skip repeat masking)
model_org=picea #select a model organism for RepBase masking in RepeatMasker
rmlib=rmlib.fa #provide an organism specific repeat library in fasta format for RepeatMasker
repeat_protein=/usr/local/opt/maker/libexec/data/te_proteins.fasta #provide a fasta file of transposable element proteins for RepeatRunner

#-----Gene Prediction
est2genome=1 #infer gene predictions directly from ESTs, 1 = yes, 0 = no

protein2genome=1 #infer predictions from protein homology, 1 = yes, 0 = no
trna=1 #find tRNAs with tRNAscan, 1 = yes, 0 = no

#-----External Application Behavior Options
cpus=4 #max number of cpus to use in BLAST and RepeatMasker (not for MPI, leave 1 when using MPI)

#-----MAKER Behavior Options
est_forward=1 #map names and attributes forward from EST evidence, 1 = yes, 0 = no
single_exon=1 #consider single exon EST evidence when generating annotations, 1 = yes, 0 = no
</code></pre>
</div></div>