Página 1 de 1

bioinfo: Leyendo FASTQ con CPAN

NotaPublicado: 2017-10-08 18:05 @795
por explorer
«Buenas. Hace un par de años describía en otra entrada cómo leer de manera eficiente ficheros FASTQ con ayuda de kseq.h Para ello definí una clase C++ que llamábamos desde Perl5 con Inline::CPP. La entrada de hoy es para contar que se puede hacer lo mismo instalando el módulo Bio::DB::HTS::Kseq desde CPAN, previa instalación en tu sistema de la librería htslib».

Sintáxis: [ Descargar ] [ Ocultar ]
Using perl Syntax Highlighting
  1. use strict;
  2. use Bio::DB::HTS::Kseq;
  3.  
  4. my ($length,$header,$sequence,$quality);
  5.  
  6. my $kseq = Bio::DB::HTS::Kseq->new("sample18MB.fq.gz");
  7. my $iter = $kseq->iterator();
  8. while(my $r = $iter->next_seq()) {
  9.  
  10.   ($header,$sequence,$quality) = ($r->name,$r->seq,$r->qual);
  11.   print ">$header\n$sequence\n$quality\n";
  12. }
Coloreado en 0.002 segundos, usando GeSHi 1.0.8.4

Artículo