This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Incorrect scaffold names GeneMark-ES

Hey everyone,

I've recently used GeneMark to de novo predict genes in an assembly I'm working on. I installed GeneMark-ES / ET v.4.33 .

The run generated a GTF, however the first column is using dufferent scaffold names to my own. My scaffolds are simply labelled "scaffold1, scaffold2, scaffold3 etc.". An example of the GeneMark-ES GTF output is below:

1_dna GeneMark.hmm exon 137 283 0 + . gene_id "1_g"; transcript_id "1_t"; 1_dna GeneMark.hmm CDS 137 283 . + 1 gene_id "1_g"; transcript_id "1_t"; 1_dna GeneMark.hmm exon 307 344 0 + . gene_id "1_g"; transcript_id "1_t"; 1_dna GeneMark.hmm CDS 307 344 . + 1 gene_id "1_g"; transcript_id "1_t"; 1_dna GeneMark.hmm exon 371 543 0 + . gene_id "1_g"; transcript_id "1_t";

When Ideally I'd want something like (note I don't know if 1_dna maps to scaffold 1 - this is just an example):

scaffold1 GeneMark.hmm exon 137 283 0 + . gene_id "1_g"; transcript_id "1_t"; scaffold1 GeneMark.hmm CDS 137 283 . + 1 gene_id "1_g"; transcript_id "1_t"; scaffold1 GeneMark.hmm exon 307 344 0 + . gene_id "1_g"; transcript_id "1_t"; scaffold1 GeneMark.hmm CDS 307 344 . + 1 gene_id "1_g"; transcript_id "1_t"; scaffold1 GeneMark.hmm exon 371 543 0 + . gene_id "1_g"; transcript_id "1_t";

My command was simply:

gmes_petap.pl --ES --cores 6 --sequence assembly.fasta

Any tips/insights into what might be happening would be appreciated.

Cheers

annotation ab initio prediction
• 1,555 views • link
written by alastair.j.lud • alastair.j.lud

1 answer

Digging around I found file dna.trace in the GeneMark generated info directory. This contains mapping information between the new scaffold names (i.e 1_dna) to the original scaffolds.

Not sure why the output didn't map them back to the originals, but the information is there.

• 0 views • link
written by alastair.j.lud • alastair.j.lud

Log in to answer this question.

More posts like this
  • Gene Specific coverage from WGS data
    written by Adyasha • Adyasha
    Hello everyone, I have annotation file like this ``` less -S Sars_cov_2.ASM985889v3.101.gtf | head -20 #!genome-build ASM985889v3 #!genome-version ASM985889v3 #!genome-date 2020-01 #!genome-build-accession NCBI:GCA_009858895.3 MN908947.3 ensembl …
  • GeneMark-ES output with error about perl module
    written by katak • katak
    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, …
  • Extracting transcript sequences with gene name (gffread)
    written by cgias • cgias
    Hello, I've extracted the transcripts using gffread from Cufflinks package using gffread -w transcripts.fa -g /path/to/genome.fa transcripts.gtf It gave me the sequences with the transcript_id …
  • STAR index with NCBI files
    written by cgias • cgias
    Hello, So when I finally made all the other threads I posted here work, I was told to change my genome & annotation file.... as …