#!/usr/bin/perl -w

# Ejercicio 14. Trabajo con GenBank

 

# Abrir el fichero GenBank.gb para lectura

 

open (GEN, "GenBank.gb");

@todo = <GEN>;

close (GEN);

 

# Dentro de las líneas del fichero buscamos una que tenga al comienzo la palabra ORIGIN,

# esta me da la separación entre la información que hay que guardar en ambos ficheros.

 

for ($i = 0; $i < scalar @todo; $i++) {

      if ($todo[$i] =~ /ORIGIN/) {

            $linea_separacion = $i;

print "La linea de separacion es la $linea_separacion\n\n";

            last;

            }

      }

 

 

open (NOSEC, ">fich1.txt");

 

for ($i = 0; $i < $linea_separacion; $i++) {

      print NOSEC "$todo[$i]";

      }

close (NOSEC);

 

 

open (SEC, ">fich2.txt");

for ($i = $linea_separacion+1; $i < scalar @todo; $i++) {

      print SEC "$todo[$i]";

      if ($todo[$i] =~ /tca/) {

            print "La cadena tca aparece en la línea $i en la posicon $-[0]\n";

            }

      }

 

close (SEC);

 

# Salir del programa

 

exit;