Gff3 - moving an annotation from scaffold level to pseudomolecule level
2
2
Entering edit mode
2.3 years ago

I have an annotation for a genome assembly in scaffold format in gff3/fasta format - someone else has ordered these scaffolds into pseudomolecules, so now I want to move the 'old' annotation in gff3 format from their scaffold positions into the pseudomolecules' position.

I could write a script to do this but I'm sure I'll introduce off-by-one errors everywhere. Is there any package out there which takes two gff3 files (gene positions on scaffolds, and scaffold positions in pseudo-molecules) and generates a new annotation gff3 with the gene positions on pseudomolecules?

Assembly annotation • 653 views
ADD COMMENT
4
Entering edit mode
2.3 years ago
kfletcher ▴ 40

Would recommend flo https://github.com/wurmlab/flo

It doesn't do exactly as you request. From memory, flo instead aligns the two assemblies and produces a new GFF file on the second assembly.

One thing, the gene lines have to be removed from the GFF file. They provide a script to do so. I don't think they provide a script to put the gene lines back in.

ADD COMMENT
1
Entering edit mode
2.3 years ago
Carambakaracho ★ 2.5k

Hi Philipp, might liftOver do what you want?

Though my impression is that data prep might just take as long as writing the script.

ADD COMMENT

Login before adding your answer.

Traffic: 2818 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6