<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; color: rgb(0, 0, 0); font-size: 14px; font-family: Calibri, sans-serif;"><div>Nothing in the scoring or gene selection has changed.</div><div><br></div><div>Changes are:</div><div>Fix trnascan naming so codon is included in name</div><div>Fix fgenesh parsing when used with correct_est_fusion</div><div>Fix final ID bug when '/' character used in GFF3 input ID.</div><div>Fix a start codon issue that could come up under when the right set of parameters were used (primarily correct_est_fusion and protein2genome).</div><div><br></div><div><br></div><div>If you can provide both gff3 outputs form comparison, I could probably tell you why. Set up both runs to make sure that settings are indeed identical.</div><div><br></div><div>--Carson</div><div><br></div><span id="OLK_SRC_BODY_SECTION"><div style="font-family:Calibri; font-size:11pt; text-align:left; color:black; BORDER-BOTTOM: medium none; BORDER-LEFT: medium none; PADDING-BOTTOM: 0in; PADDING-LEFT: 0in; PADDING-RIGHT: 0in; BORDER-TOP: #b5c4df 1pt solid; BORDER-RIGHT: medium none; PADDING-TOP: 3pt"><span style="font-weight:bold">From: </span> Shaun Jackman <<a href="mailto:sjackman@gmail.com">sjackman@gmail.com</a>><br><span style="font-weight:bold">Reply-To: </span> Shaun Jackman <<a href="mailto:sjackman@gmail.com">sjackman@gmail.com</a>><br><span style="font-weight:bold">Date: </span> Monday, May 5, 2014 at 6:09 PM<br><span style="font-weight:bold">To: </span> "<a href="mailto:maker-devel@yandell-lab.org">maker-devel@yandell-lab.org</a>" <<a href="mailto:maker-devel@yandell-lab.org">maker-devel@yandell-lab.org</a>><br><span style="font-weight:bold">Subject: </span> [maker-devel] Fewer genes in MAKER 2.31.3<br></div><div><br></div><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>
_______________________________________________
maker-devel mailing list
<a href="mailto:maker-devel@box290.bluehost.com">maker-devel@box290.bluehost.com</a>
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a>
</span></body></html>