[Bioperl-l] Empty FASTA files with Bio::SeqIO

Aaron J Mackey ajm6q@virginia.edu
Thu, 28 Dec 2000 21:44:59 -0500 (EST)


On Wed, 20 Dec 2000, Hilmar Lapp wrote:

> The correctly FASTA-formatted empty seq ought to have an empty line after
> the '>'-line. I think we should check for that, just to be sure we're not
> misinterpreting something.

The fasta programs accept:

>sequence 1
>sequence 2
ATCGCGCA
>sequence 3
>
>sequence 5

>sequence 6
GATTACA

Note the "sequence 4" that has no description at all.

Of course these examples represent garbage input, but I just wanted to
clarify *all* the possibilities.

-Aaron

-- 
 o ~   ~   ~   ~   ~   ~  o
/ Aaron J Mackey           \
\  Dr. Pearson Laboratory  /
 \ University of Virginia  \
 /  (804) 924-2821          \
 \  amackey@virginia.edu    /
  o ~   ~   ~   ~   ~   ~  o