[Biopython] GFF3Writer

Michal mictadlo at gmail.com
Sat Apr 9 04:02:13 EDT 2011


On 04/08/2011 10:06 PM, Brad Chapman wrote:
> Thanks for the feedback. I added an example of how to build up
> SeqRecord and SeqFeature objects and write them to GFF:
>
> http://www.biopython.org/wiki/GFF_Parsing#Writing_GFF3_from_scratch
>
> I also made a small update to the output code to support numbers and
> non-lists in the features, so to run the example code you should
> grab the latest version from:
>
> https://github.com/chapmanb/bcbb/tree/master/gff
>
> Hope this helps,
> Brad
Thank you it works. How is it possible to create ##sequence-region with 
the writer like described in http://www.sequenceontology.org/gff3.shtml:

##gff-version   3
##sequence-region   ctg123 1 1497228
ctg123 . gene            1000  9000  .  +  .  ID=gene00001;Name=EDEN
ctg123 . TF_binding_site 1000  1012  .  +  .  Parent=gene00001
ctg123 . mRNA            1050  9000  .  +  .  ID=mRNA00001;Parent=gene00001

Thank you in advance.

Michal


More information about the Biopython mailing list