[Biopython] GFF.writer

Brad Chapman chapmanb at 50mail.com
Fri May 3 06:51:37 EDT 2013


Mic;

> I found here ( http://www.biopython.org/wiki/GFF_Parsing#Writing_GFF3 ) an
> example how to write gff3 files.
>
> Is it possible not to write "##sequence-region ID1 1 20"?

It only writes that line if you supply a sequence for the SeqRecord:

https://github.com/chapmanb/bcbb/blob/master/gff/BCBio/GFF/GFFOutput.py#L105

so manually removing the sequence before writing will prevent it from
outputting that metadata. However, the GFF3 spec encourages writing
these metadata items:

http://www.sequenceontology.org/gff3.shtml

Brad


More information about the Biopython mailing list