[maker-devel] Names/IDs from pred_gff are included in final gff
Ole Kristian Tørresen
o.k.torresen at ibv.uio.no
Wed Jan 17 10:52:13 MST 2018
Ok, but I have an entry in the final gff like this:
ID=maker-GmG20150304_scaffold_2371-pred_gff_maker-gene-0.125-mRNA-1;Parent=maker-GmG20150304_scaffold_2371-pred_gff_maker-gene-0.125;Name=ENSGMOT00000000668.1;_AED=0.00;_eAED=0.00;_QI=819|1|1|1|1|1|4|112|726;score=89.75616
(The name is derived from a pred_gff entry which is the results of mapping an old annotation to the new assembly).
This is then called ENSGMOT00000000668.1 protein AED:0.00 eAED:0.00 QI:819|1|1|1|1|1|4|112|726 in the proteins.fasta file. Which is unfortunate, because it apparently mapped 12 places in the assembly.
I have set map_forward=0, but keep_preds=1 (filtering on domain presence and AED score later).
This and another file (result of genemark_gtf2gff3), is not input as match/match_part to MAKER, but with gene/exon/CDS/mRNA. Could that be the issue?
Ole
> On 16 Jan 2018, at 19:15, Carson Holt <carsonhh at gmail.com> wrote:
>
> pred_gff will maintain it’s name in the match/match_part feature as the information is pulled directly from the input GFF3. But any pred_gff feature that becomes a final model will be renamed to something like “scaffold_1517-pred_gff_GeneMark.hmm-gene-0.6-mRNA-1” unless you specify map_forward=1 to maintain the original name.
>
> —Carson
>
>
>
>> On Jan 4, 2018, at 6:21 AM, Ole Kristian Tørresen <o.k.torresen at ibv.uio.no> wrote:
>>
>> Hi,
>> as far as I can see, names or IDs of features in gffs given to pred_gff is included in the final output as the name of the feature. As far as I can understand, this is not expected behaviour (it is for model_gff). This is with MAKER 2.31.9.
>>
>> I have these settings:
>> map_forward=0
>> keep_preds=1
>>
>> I thought that map_forward had to be 1 to get the names for the old GFFs. Can you replicate this?
>>
>> Thank you.
>>
>> Sincerely,
>> Ole K. Tørresen
>>
>> _______________________________________________
>> maker-devel mailing list
>> maker-devel at box290.bluehost.com
>> http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org
>
More information about the maker-devel
mailing list