Question: SNAP error in MAKER exon:out_of_bounds
0
gravatar for preethibagopi5
3 days ago by
preethibagopi50 wrote:

Hi, I am running a genome annotation in MAKER2 in which I completed the evidence round of maker and got a GFF3 file as output. I used the GFF3 file from the first maker run to train the SNAP in which it is showing me the error as this.

MODEL20592 1 1 3 + errors(3): exon-1:out_of_bounds exon-2:out_of_bounds exon-3:out_of_bounds
MODEL20593 1 1 3 - errors(3): exon-3:out_of_bounds exon-2:out_of_bounds exon-1:out_of_bounds
MODEL20594 1 1 3 - errors(3): exon-3:out_of_bounds exon-2:out_of_bounds exon-1:out_of_bounds
MODEL58264 1 1 5 + errors(5): exon-1:out_of_bounds exon-2:out_of_bounds exon-3:out_of_bounds exon-4:out_of_bounds exon-5:out_of_bounds
MODEL58265 1 1 3 - errors(3): exon-3:out_of_bounds exon-2:out_of_bounds exon-1:out_of_bounds
MODEL58266 1 1 2 - errors(2): exon-2:out_of_bounds exon-1:out_of_bounds
MODEL58267 1 1 4 - errors(4): exon-4:out_of_bounds exon-3:out_of_bounds exon-2:out_of_bounds exon-1:out_of_bounds
MODEL58268 1 1 4 - errors(4): exon-4:out_of_bounds exon-3:out_of_bounds exon-2:out_of_bounds exon-1:out_of_bounds
MODEL58269 1 1 1 - errors(1): exon-1:out_of_bounds
MODEL39568 1 1 1 - errors(1): exon-1:out_of_bounds
***10652 genes, 0 OK, 0 warnings, 10652 errors***

All the gene models are showing errors. I checked with NCBI provided GFF3 file also to train SNAP but I ended up in same error telling that all gene models as error.

**Command Used*****

maker2zff genome.all.gff
fathom -validate geneome.ann genome.dna > snap_validate.output.txt

Can anyone please help me to solve this problem... I don't know whether this is the right place to ask? But I rised my question in maker forum also i didn't get any solution Thanks in advance.

maker2 snap annotation genome • 59 views
ADD COMMENTlink modified 2 days ago • written 3 days ago by preethibagopi50

Are you sure you're using the same genome assembly fasta with the same gffs? This looks like you have exons with coordinates longer than your assembly contigs/chromosomes, which normally only happens if it's the wrong one

ADD REPLYlink written 2 days ago by Philipp Bayer6.1k

Yes I am using the same GFF file from same genome assembly file..... A part of my GFF file:

Scaffold_94%3BHRSCAF%3D190 . contig 1 59085 . . .
Scaffold_94%3BHRSCAF%3D190 maker exon 4576 4872 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4576 4938 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4576 5007 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4519 4552 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4628 4938 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4501 4516 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4604 4938 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4501 4515 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4591 4938 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4479 4549 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4635 4938 . - .
Scaffold_94%3BHRSCAF%3D190 maker exon 4519 4573 . - .

ADD REPLYlink modified 2 days ago • written 2 days ago by preethibagopi50

out of curiosity: what are the %-signs doing in the GFF file?

ADD REPLYlink written 2 days ago by lieven.sterck5.1k

Name of the scaffold is having the % symbol in it. Is that will be the problem.....?

ADD REPLYlink written 2 days ago by preethibagopi50

Not sure but I have personally never encountered that before

ADD REPLYlink written 2 days ago by lieven.sterck5.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 562 users visited in the last hour