[Bioperl-l] Write segments with Bio::Tools::GFF
daniel.standage at gmail.com
Wed Nov 3 16:14:02 EDT 2010
I'm using a tool called AnnotationSketch (part of the GenomeTools package)
to create graphics from my GFF3 files. It complains when sequence-region
lines are not included, even when the reference is printed as a feature.
That's what got me interested in the problem in the first place.
Thanks for the info, I'll figure something out.
On Wed, Nov 3, 2010 at 3:09 PM, Scott Cain <scott at scottcain.net> wrote:
> Hi Daniel,
> Why do you need the sequence-region line? It is mostly just
> informational. If you want to define the reference sequence for the
> features in the GFF file, I would suggest printing out a full GFF line
> for the reference sequence instead anyway. That way, it can have all
> the information that a GFF line can encode, like a source and type, as
> well as information in the ninth column.
> That said, Bio::Tools::GFF probably does lack that method. The code
> in it is probably a little outdated.
> On Wed, Nov 3, 2010 at 2:42 PM, Daniel Standage
> <daniel.standage at gmail.com> wrote:
> > Hi
> > I've got a simple script that is filtering some GFF3 data. The
> > Bio::Tools::GFF class has methods for reading and writing features and
> > reading regions (segments), but I cannot find a method for writing
> > It's not hard to just print out the sequence-region line from the
> > LocatableSeq object, but then these lines are printed out before the
> > gff-version line (which is a no-no). Any suggestions about how to handle
> > this?
> > Thanks.
> > --
> > Daniel S. Standage
> > Graduate Research Assistant
> > Bioinformatics and Computational Biology Program
> > Iowa State University
> > _______________________________________________
> > Bioperl-l mailing list
> > Bioperl-l at lists.open-bio.org
> > http://lists.open-bio.org/mailman/listinfo/bioperl-l
> Scott Cain, Ph. D. scott at scottcain dot
> GMOD Coordinator (http://gmod.org/) 216-392-3087
> Ontario Institute for Cancer Research
Daniel S. Standage
Graduate Research Assistant
Bioinformatics and Computational Biology Program
Iowa State University
More information about the Bioperl-l