Foro - Perl en Español

por **abraham03** » 2016-09-12 14:32 @647

Hola, he creado este script para extraer múltiples secuencias contenidas en una carpeta en formato fastq y pasarlas a formato fasta y que mantenga el mismo nombre de cada archivo fastq:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/perl
 
use strict;
use warnings;
use Getopt::Long;
 
my ($dir, $files, $file_name, $file_format, $extension);
 
    GetOptions (
                'dr=s'  =>\$dir
                );
 
foreach $files (glob("$dir/*.fastq")) {
        open LINE, "<", $files or die "can't read open $files";
            while(
                    defined(my $head    = <LINE>)       &&
                    defined(my $seq     = <LINE>)       &&
                    defined(my $qhead   = <LINE>)       &&
                    defined(my $quality = <LINE>)           
            ){
                        substr($head, 0, 1, '>'); 
                        
                        $extension = "\\.fastq";
                        $files =~ s/(.*)$extension(.*)/$1/;
                        $file_name =$files;
                        open OUTFILE, '>>', $file_name.".fasta";
                        print OUTFILE $head, $seq;
            }
}
close LINE;
close OUTFILE;
exit;
 
Coloreado en 0.003 segundos,  usando GeSHi 1.0.8.4

El código funciona bien, solo que quiero agregarle más opciones. Una de ellas es convertir un solo archivo si así se desea, para lo cual quisiera saber cómo organizar en una subrutina la siguiente parte del código:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

 while(
                    defined(my $head    = <LINE>)       &&
                    defined(my $seq     = <LINE>)       &&
                    defined(my $qhead   = <LINE>)       &&
                    defined(my $quality = <LINE>)           
            ){
                        substr($head, 0, 1, '>'); 
                        
                        $extension = "\\.fastq";
                        $files =~ s/(.*)$extension(.*)/$1/;
                        $file_name =$files;
                        open OUTFILE, '>>', $file_name.".fasta";
                        print OUTFILE $head, $seq;
            }
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

Para usarlo en las dos opciones (carpeta con múltiples fastq y un solo archivo fastq). Actualmente tengo que usar dos scripts para las diferentes opciones y quisiera tenerlo en uno solo y creo que usar una subrutina me evitará tener código repetido.

¡¡Muchas gracias!!

#!/usr/bin/env perl
use strict;
use warnings;
use Getopt::Long;
 
my ($dir, $file_name, $file_format);
my $extension = "\\.fastq";
 
GetOptions (
    'dr=s' => \$dir
);
 
if (-f $dir) {
    procesa($dir);
}
else {
    foreach my $file (glob("$dir/*.fastq")) {
        procesa($file);
    }
}
 
 
sub procesa {
    my $file = shift;
 
    ($file_name = $file) =~ s/(.*)$extension.*/$1/;
 
    open LINE,    '<',   $file             or die "can't read open $file\n";
    open OUTFILE, '>>', "$file_name.fasta" or die "can't write $file_name\n";
 
    while (
        defined(my $head    = <LINE>)       &&
        defined(my $seq     = <LINE>)       &&
        defined(my $qhead   = <LINE>)       &&
        defined(my $quality = <LINE>)
    ) {
        print OUTFILE ">$head", $seq;
    }
 
    close OUTFILE;
    close LINE;
}
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

por **abraham03** » 2016-09-14 15:00 @667

Como siempre, mil gracias, no solo por contestar, ¡¡¡ si no también por hacer las cosas sencillas !!!

Foro - Perl en Español

De fastaq a fasta en subrutina

De fastaq a fasta en subrutina

Publicidad

Re: De fastaq a fasta en subrutina

Re: De fastaq a fasta en subrutina

¿Quién está conectado?