Extraer secuencia
Publicado: 2012-07-05 09:41 @445
Hola, a partir de un fichero con el siguiente formato:
Estoy intentando extraer, en otro fichero (out.txt), el primer supercont. Para ello, en teoría, (según lo que yo pretendía) bastaría con teclear 1.2 y dejar que ese valor se guardara en $a. Al hacer esto, en vez de obtener un fichero con sólo la primera secuencia, obtengo uno con todas las secuencias.
Seguramente sea algo muy tonto pero no tengo mucha experiencia programando en Perl , ¿alguien puede ayudarme?
Este es el código del programa:
Using text Syntax Highlighting
>supercont1.1 of E.coli M1.001
TTATACTTAACTTATTATTTAAATTAATATAGCTTATTATATAGTAACCTTATAGTTCTA
ATTATATTATAAATAATTTATTTATTATACTTATTTAGTAGTTATTAAATATCTTTATTT
>supercont1.2 of E.coli M1.001
ATATAATAAATAAGAAGTATATATAAAGTAATAATTTTTTATTTATTTATAAGAATTATT
AAAGTTAGTTAATTTATAGTTATTATAACTATTAAAGATTAAAAGTTAATAAGTACTTTT
(...)
TTATACTTAACTTATTATTTAAATTAATATAGCTTATTATATAGTAACCTTATAGTTCTA
ATTATATTATAAATAATTTATTTATTATACTTATTTAGTAGTTATTAAATATCTTTATTT
>supercont1.2 of E.coli M1.001
ATATAATAAATAAGAAGTATATATAAAGTAATAATTTTTTATTTATTTATAAGAATTATT
AAAGTTAGTTAATTTATAGTTATTATAACTATTAAAGATTAAAAGTTAATAAGTACTTTT
(...)
Coloreado en 0.000 segundos, usando GeSHi 1.0.8.4
Estoy intentando extraer, en otro fichero (out.txt), el primer supercont. Para ello, en teoría, (según lo que yo pretendía) bastaría con teclear 1.2 y dejar que ese valor se guardara en $a. Al hacer esto, en vez de obtener un fichero con sólo la primera secuencia, obtengo uno con todas las secuencias.
Seguramente sea algo muy tonto pero no tengo mucha experiencia programando en Perl , ¿alguien puede ayudarme?
Este es el código del programa:
Using perl Syntax Highlighting
- #!usr/local/bin/perl
- use strict;
- use warnings;
- print "write the number of the supercontig from out which you want to start erasing";
- my $a=<STDIN>;
- chomp $a;
- my $file="E_coli_m1.001_1_supercontigs.fasta";
- my $out_file="out.txt";
- my $line;
- open (IN, "<$file") or die "can't open $file :$!\n";
- open (OUT, ">$out_file") or die "can't open $out_file:$!\n";
- while ($line= <IN>) {
- if ($line =~m/\d$a/) {
- exit;}
- else {print OUT $line;
- }
- }
- close IN;
- close OUT;
- exit;
Coloreado en 0.002 segundos, usando GeSHi 1.0.8.4