Question: Gff3 - moving an annotation from scaffold level to pseudomolecule level
0
gravatar for Philipp Bayer
4 months ago by
Philipp Bayer6.1k
Australia/Perth/UWA
Philipp Bayer6.1k wrote:

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?

annotation assembly • 245 views
ADD COMMENTlink modified 4 months ago by Carambakaracho1.2k • written 4 months ago by Philipp Bayer6.1k
1
gravatar for kfletcher
4 months ago by
kfletcher10
kfletcher10 wrote:

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 COMMENTlink written 4 months ago by kfletcher10
1
gravatar for Carambakaracho
4 months ago by
Carambakaracho1.2k
Switzerland/Basel
Carambakaracho1.2k wrote:

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 COMMENTlink written 4 months ago by Carambakaracho1.2k
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: 1613 users visited in the last hour