[Bioperl-l] Fasta sequence width question
Brian Osborne
bosborne11 at verizon.net
Wed Feb 22 19:28:12 UTC 2012
Francis,
You can set this yourself.
http://search.cpan.org/~cjfields/BioPerl-1.6.901/Bio/SeqIO/fasta.pm#width
Brian O.
On Feb 22, 2012, at 1:55 PM, Francisco J. Ossandón wrote:
> Hello,
> I have a question about the width used by Bioperl in the Fasta format. This
> format recommends that lines of text be shorter than 80 characters, but
> there is no really fixed length for the sequence lines.
> http://en.wikipedia.org/wiki/FASTA_format
> http://blast.ncbi.nlm.nih.gov/blastcgihelp.shtml
>
> I usually download sequences from NCBI and their Fasta sequences always use
> a sequence line width of 70 characters (like
> http://www.ncbi.nlm.nih.gov/protein/50842445?report=fasta ), while Bioperl
> uses a default of 60 to write in Fasta format using Bio::SeqIO::Fasta
> ("BEGIN { $WIDTH = 60}"). Is there a particular reason to set the width at
> 60 characters?? Or maybe it could be changed to 70 to match NCBI??
>
> Cheers,
>
> Francisco J. Ossandon
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
More information about the Bioperl-l
mailing list