[Bioperl-l] Next-gen modules
Chris Fields
cjfields at illinois.edu
Fri Jul 24 18:27:32 UTC 2009
On Jul 24, 2009, at 8:00 AM, Peter wrote:
> Hi all,
>
> On Fri, Jul 24, 2009 at 1:19 PM, Chris Fields<cjfields at illinois.edu>
> wrote:
>>>
>>> Have you guys (BioPerl) have also gone for "fastq-sanger" instead of
>>> just "fastq" for the Sanger Standard version of FASTQ (like EMBOSS)?
>>> Does BioPerl use just "fastq" to mean anything?
>>
>> Short answer: yes, and yes.
>>
>> Slightly longer answer: I've set up SeqIO so it converts "new(-
>> format =>
>> 'foo-bar')" to new(-format => 'foo, -variant => 'bar'). In the fastq
>> constructor, if the variant is expected but isn't defined (i.e. for
>> 'fastq')
>> it defaults to sanger. Makes it a bit easier maintenance-wise if a
>> new
>> variant pops up.
>
> Right, so BioPerl understands "fastq" and "fastq-sanger" to mean the
> Sanger standard FASTQ files.
Yes.
> I've just updated Biopython to also allow "fastq-sanger" as an alias
> for
> "fastq", so we are consistent here:
> http://lists.open-bio.org/pipermail/biopython-dev/2009-July/
> 006466.html
>
> Biopython, BioPerl and EMBOSS now all agree on the format names:
> * "fastq-sanger" - PHRED scores offset 33
> * "fastq-solexa" - Solexa scores offset 64
> * "fastq-illumina" - PHRED scores offset 64
>
> And Biopython and BioPerl also agree on the meaning of "fastq" as
> an alias for "fastq-sanger". Unfortunately EMBOSS differs here, see:
> http://lists.open-bio.org/pipermail/emboss-dev/2009-July/000599.html
>
> Does BioJava or BioRuby have a SeqIO equivalent where they need
> to give different sequence formats unique names? If so, we should
> talk to them soon...
>
> Peter
Not sure, but it would be nice to have consistency there, yes.
chris
More information about the Bioperl-l
mailing list