[Bioperl-l] warning: Bio::Index::Fastq;

Peter Cock p.j.a.cock at googlemail.com
Tue Mar 11 16:12:00 UTC 2014


On Tue, Mar 11, 2014 at 4:01 PM, EpiMAN <eh_mcb at yahoo.com> wrote:
> test <http://bioperl.996286.n3.nabble.com/file/n17374/test>
>
> grep -C 10 "^@@DDDDDFFFFFIIIDCBBBB"

Very helpful - it does look like it came from a valid FASTQ
(the -- lines are where grep cuts each snippet).

This does seem like an indexing bug where a quality line
starts with '@' signs, but like Chris I thought this had been
fixed once before?

https://github.com/bioperl/bioperl-live/blob/master/Bio/Index/Fastq.pm

Peter



More information about the Bioperl-l mailing list