This is a test version of Biostars. For the public version, visit https://www.biostars.org.
GeneMark-ES output with error about perl module

Hi. I am running GeneMark-ES for ab initio gene annotation of insect genomes.

After GeneMark-ES completed, I get the seemingly final output file named genemark.gtf, but the following standard output error is also generated.

(in cleanup) (in cleanup) at /path-to/lib/site_perl/5.26.2/Object/InsideOut.pm line 1953 during global destruction.

The actual standard output error is the repeat of the above.

The genemark.gtf is like this:

$ tail genemark.gtf
scaffold_209549 GeneMark.hmm    CDS 241 306 .   +   2   gene_id "188161_g"; transcript_id "188161_t"; cds_type "Internal"; count "4_7";
scaffold_209549 GeneMark.hmm    intron  307 328 0   +   1   gene_id "188161_g"; transcript_id "188161_t"; count "4_6";
scaffold_209549 GeneMark.hmm    exon    329 361 0   +   .   gene_id "188161_g"; transcript_id "188161_t";
scaffold_209549 GeneMark.hmm    CDS 329 361 .   +   2   gene_id "188161_g"; transcript_id "188161_t"; cds_type "Internal"; count "5_7";
scaffold_209549 GeneMark.hmm    intron  362 395 0   +   1   gene_id "188161_g"; transcript_id "188161_t"; count "5_6";
scaffold_209549 GeneMark.hmm    exon    396 428 0   +   .   gene_id "188161_g"; transcript_id "188161_t";
scaffold_209549 GeneMark.hmm    CDS 396 428 .   +   2   gene_id "188161_g"; transcript_id "188161_t"; cds_type "Internal"; count "6_7";
scaffold_209549 GeneMark.hmm    intron  429 457 0   +   1   gene_id "188161_g"; transcript_id "188161_t"; count "6_6";
scaffold_209549 GeneMark.hmm    exon    458 517 0   +   .   gene_id "188161_g"; transcript_id "188161_t";
scaffold_209549 GeneMark.hmm    CDS 458 517 .   +   2   gene_id "188161_g"; transcript_id "188161_t"; cds_type "Internal"; count "7_7";

The point which I think strange is that 3rd line in the output file above doesn't contain "mRNA" or "gene", unlike the standard gtf file.

gmes.log ends with: /path-to/gmes_petap.pl : [Sun Jul 19 22:22:38 2020] mv genemark.gtf.tmp genemark.gtf

Is the output file I have got all right? Thank you very much in advance.

Best,

genome annotation ab initio genemark

0 answers

No answers yet.

Log in to answer this question.