[Biopython] gff3 writting

Brad Chapman chapmanb at 50mail.com
Tue Jul 2 09:26:17 UTC 2013


Mic;
Thanks for the feedback, comments below.

> I found here (
> http://biopython.org/wiki/GFF_Parsing#Writing_GFF3_from_scratch ) an
> example how to write GFF3 from scratch.
>
> I modified it in order to add one more features and sub_features, but the
> second sub_features are not visible:
[...]
> with the following code:
[...]
> top_feature2.sub_features2 = [SeqFeature(FeatureLocation(30, 35),
> type="exon", strand=1,
>                                        qualifiers=sub_qualifiers2),
>                             SeqFeature(FeatureLocation(37, 40),
> type="exon", strand=1,
>                                        qualifiers=sub_qualifiers2)]

You want to specify these as the `sub_features` attributes (not
`sub_features2`). Hope this helps sort it out,
Brad



More information about the Biopython mailing list