Re: Conversion of Phred 33 -> Phred 64 quality
"Fields, Christopher J" <[email protected]>
| Newsgroups | gmane.comp.lang.perl.bio.general |
|---|---|
| Message-ID | <[email protected]> |
This would be the so-called fastq-illumina’ or ‘fastq-solexa’ variants. FASTQ does handle these but the output for modern Illumina pipelines is to generate Sanger-based FASTQ.
Bio::SeqIO::fastq can interconvert these but it’s painfully slow last time I checked (mainly in part to the underlying sequence quality classes). They’re not really designed for tons of data and needs refactoring to deal with large sequence data sets. Lately, frankly, I’ve been using non* methods for doing this, namely seqtk:
https://github.com/lh3/seqtk
chris
> On Jan 26, 2015, at 8:23 AM, Mark A. Jensen <[email protected]> wrote:
>
> Hi folks,
> I know I could RTFM, but maybe someone knows off the top of their head: I understand that Illumina at one time made a switch in the constant added to quality scores to generate the FASTQ that comes off their instruments. This leads to a certain incomparability of data before and after that switch. This is about all I know of the issue; does anyone here have experience with this? Are there any BP modules that do this translation?
> much appreciated-
> MAJ
> _______________________________________________
> Bioperl-l mailing list
> [email protected]
> http://mailman.open-bio.org/mailman/listinfo/bioperl-l
_______________________________________________
Bioperl-l mailing list
[email protected]
http://mailman.open-bio.org/mailman/listinfo/bioperl-l