[Open-bio-l] FASTQ support in Biopython, BioPerl, and EMBOSS
Chris Fields
cjfields at illinois.edu
Fri Jul 31 01:19:56 UTC 2009
I do tend to agree, and I don't think any savings from a performance
hit will be worth the headache of having to repeatedly explain why
it's (silently) doing so, when a simple warning or error message
('value X out of range for fastq format y') would suffice.
chris
On Jul 30, 2009, at 6:52 PM, Aaron Mackey wrote:
> I would strongly warn against truncation, for any reason. Use the
> formulas you have for quality-encoding conversions, but do not
> assume that you know more than I do about what my data contains, or
> that you are in any way helping me by altering my data, silently or
> otherwise. Said another way, feel free to warn me that my data may
> contain garbage, and utterly fail to convert it for me, but do not
> try to fix it for me.
>
> -Aaron
>
> On Thu, Jul 30, 2009 at 5:50 PM, Peter <biopython at maubp.freeserve.co.uk
> > wrote:
> On Thu, Jul 30, 2009 at 9:08 PM, Chris Fields<cjfields at illinois.edu>
> wrote:
> >
> > I do think if it affects performance to a significant enough
> degree we
> > can do this silently, we just need to ensure this is well-
> documented.
>
> Agreed.
>
> > My opinions is this use will prove to be a edge case anyway (most
> will
> > want conversion to Sanger vs. Illumina/Solexa).
>
> Absolutely.
>
> Going from Solexa/Illumina to Sanger FASTQ will be more common
> (and there are no truncation issues). Going from Sanger FASTQ to
> Solexa or Illumina FASTQ will be rarer, and while a truncation is
> possible it requires very high scores (above PHRED 62) which are
> likely only to be possible from a consensus alignment or such like.
> i.e. Yes, it should be an edge case.
>
> I guess this expected usage supports the argument about issuing a
> warning on truncation, even with a modest performance overhead
> (because it only slows down the rarer expected usage).
>
> But let's get some benchmarks done to help settle this...
>
> Peter
> _______________________________________________
> Open-Bio-l mailing list
> Open-Bio-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/open-bio-l
>
More information about the Open-Bio-l
mailing list