[Biopython] Adaptor trimmer and dimers

Peter biopython at maubp.freeserve.co.uk
Thu Oct 15 12:20:47 EDT 2009


On Thu, Oct 15, 2009 at 5:00 PM, natassa <natassa_g_2000 at yahoo.com> wrote:
> Hallo Biopythoners,
> I followed a recent thread conversation about adaptor trimming,
> which I intend to do on Illumina runs, and I am not sure I know
> where exactly in github I could find Brad Chapman's code for
> trimming AFTER modifications that he has done based on the
> thread conversation. ...

I guess you mean Brad's August Blog Post:
http://bcbio.wordpress.com/2009/08/09/trimming-adaptors-from-short-read-sequences/
and the following mailing list thread which included some tips on
speeding up the Biopython side of things:
http://lists.open-bio.org/pipermail/biopython/2009-August/005417.html

For anyone else interested, there are some simple examples in the
tutorial (using SeqRecord slicing - elegant and simple, but a bit slow):
http://biopython.org/DIST/docs/tutorial/Tutorial.html#sec:FASTQ-slicing-off-primer
http://biopython.org/DIST/docs/tutorial/Tutorial.html#sec:FASTQ-slicing-off-adaptor

And I did a blog post about low level FASTQ handling for speed
at the cost of flexibility and simplicity (using some of the same
ideas from the August mailing list discussion):
http://news.open-bio.org/news/2009/09/biopython-fast-fastq/

Peter


More information about the Biopython mailing list