3
$\begingroup$

I want to convert my gff3 annotation files to genbank format for use in Mauve. I found the seqret tool here https://www.ebi.ac.uk/Tools/sfc/emboss_seqret/ which can perform this task, but my files (bacterial genomes) are too big.

Can anyone suggest a simple method to make this conversion?

$\endgroup$
2
  • 3
    $\begingroup$ Hi there. I don't know any simple method for what you ask, mainly because gff3 format does not contain sequences (only annotations) while fasta & genbank do. This thread on biostars might help tho : Gff3 + Fasta To Genbank. As the title suggest, you will need your genome gff3 and fasta reference as inputs. $\endgroup$ Jan 9 '20 at 16:56
  • 1
    $\begingroup$ THanks for taking the time to consider my question Paul. My gff3 files have the sequence in FASTA at the end of the file (following a ##FASTA indicator), but i take your point - will need some further leg work on this one! Thanks for the link also. $\endgroup$
    – Mark
    Jan 10 '20 at 14:26
2
$\begingroup$

I believe seqret is already the simplest approach. If you install the emboss command line utility, genome (file) sizes will not be a problem. You can even install it with conda, see https://anaconda.org/bioconda/emboss.

After that, you can use it as the example:

seqret -sequence {genome file} -feature -fformat gff -fopenfile {gff file} -osformat genbank -osname_outseq {output prefix} -ofdirectory_outseq gbk_file -auto

Hope it helps

$\endgroup$
1
$\begingroup$

You can use BioPerl to do this. Read in the GFF and make gene objects and write out those gene objects in Genbank format.

Perhaps this tool (doing the opposite) will get you started: https://metacpan.org/pod/bp_genbank2gff3.pl

Here is a script doing what you ask but it doesn't use any libraries (so is probably buggy): https://gist.github.com/avrilcoghlan/5386810

Hope that helps, Dan.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.