[Bioperl-l] Bio::SeqIO::fasta.pm

Hilmar Lapp hlapp@gnf.org
Fri, 23 Aug 2002 10:07:51 -0700


Apparently not.

Any votes from the core or from people on the list? I don't have a strong view here. I recall having seen somewhere that indeed the presence of the name is optional, but I'm not sure.

One could argue for tolerant reader, strict writer ...

	-hilmar

> -----Original Message-----
> From: Wiepert, Mathieu [mailto:Wiepert.Mathieu@mayo.edu]
> Sent: Friday, August 23, 2002 9:27 AM
> To: 'bioperl-l@bioperl.org'
> Subject: FW: [Bioperl-l] Bio::SeqIO::fasta.pm
> 
> 
> I was just testing something out, and it was accepting a 
> blank in the header.  If I actually wanted raw sequence, you 
> are correct of course.
> 
> I found an old thread that made me think a > with a space or 
> nothing after it was valid, which was why I asked.  I hadn't 
> found anything recent, though I could have missed it.  The 
> thread never said what was the agreed upon standard?
> 
> http://bioperl.org/pipermail/bioperl-guts-l/1999-November/001311.html
> 
> Mathieu Wiepert
> Medical Information Resources
> Mayo Foundation
> (507) 266-2317 Fax (507)-284-0360
> wiepert.mathieu@mayo.edu 
> 
> 
> > -----Original Message-----
> > From: Hilmar Lapp [mailto:hlapp@gnf.org]
> > Sent: Friday, August 23, 2002 11:23 AM
> > To: Wiepert, Mathieu
> > Cc: Bioperl
> > Subject: RE: [Bioperl-l] Bio::SeqIO::fasta.pm
> > 
> > 
> > Didn't someone post a FASTA definition document or link a 
> > while ago? AFAIK the ID is mandatory, you can't just have an 
> > empty line.
> > 
> > Have you tried reading your seqs as format raw? Not sure 
> > whether that one strips all non-seq characters, and whether 
> > it can handle multiple seqs, but generally speaking 'raw' 
> > format is what you would do if you really don't care about 
> > anything except the sequence itself.
> > 
> > 	-hilmar
> > 
> > > -----Original Message-----
> > > From: Wiepert, Mathieu [mailto:Wiepert.Mathieu@mayo.edu]
> > > Sent: Friday, August 23, 2002 7:50 AM
> > > Cc: Bioperl
> > > Subject: [Bioperl-l] Bio::SeqIO::fasta.pm
> > > 
> > > 
> > > Hi,
> > > 
> > > Is it valid to have a fasta file with no header info (other 
> > > than the >, with no spaces after).  I.e. something like
> > > >
> > > ACACACACA
> > > 
> > > 
> > > Would lead to blank primary_id, not sure what effect that 
> > > would have down the line.
> > > 
> > > I ask because fasta.pm throws the "Can't parse fasta header" 
> > > error when there *is* a space after the >, but goes merrily 
> > > on if there is a \n.  And then dies later.  Just wasn't sure 
> > > what the expected agreed upon behavior is.  
> > > 
> > > I ran into this because I was testing a few sequence blasts, 
> > > and don't need a header really.  I figure it's not much use 
> > > to have a fasta file with a bunch of "empty" header lines, 
> > > but fasta only has a comment
> > > 
> > > # FIX incase no space between > and name \AE
> > > 
> > > So I wasn't sure what the intent is.
> > > 
> > > -Mat
> > > _______________________________________________
> > > Bioperl-l mailing list
> > > Bioperl-l@bioperl.org
> > > http://bioperl.org/mailman/listinfo/bioperl-l
> > > 
> > 
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
>