Página 1 de 1

Mostrar posición de codón

NotaPublicado: 2015-06-12 18:34 @815
por perl biología
Buenas. Dada una secuencia he buscado si posee el codón ATG de inicio con los códigos:

#!/usr/bin/perl

$ secuencia = AAAATGCGGTTATGCCCATGGGGTG;
$ s = 'ATG';

$posicon = 0;
$n = 0;
$n++ while $posicion = 1 + index($secuencia, $s, $posicion);

print "SE HAN ENCONTRADO $n VECES ESTE STRING ($s) EN LA SECUENCIA.\n";


exit;

Y me muestra el número de veces que tengo la secuencia ATG.

A continuación quisiera que me mostrase en el lugar que está esa secuencia. ¿Cómo podría realizarlo?

Muchas gracias.

Re: Mostrar posición de codón

NotaPublicado: 2015-06-12 18:44 @822
por explorer
Basta con agregar un print() dentro del bucle while():
Sintáxis: [ Descargar ] [ Ocultar ]
Using perl Syntax Highlighting
  1. #!/usr/bin/perl
  2.  
  3. $secuencia = 'AAAATGCGGTTATGCCCATGGGGTG';
  4. $s         = 'ATG';
  5.  
  6. $posicion = 0;
  7. $n = 0;
  8.  
  9. while ($posicion = 1 + index($secuencia, $s, $posicion)) {
  10.     $n++;
  11.     print "Se ha encontrado en la posición $posición\n";      # informa de la posición, basada en '1'
  12. }
  13. print "SE HAN ENCONTRADO $n VECES LA SECUENCIA ($s) EN LA SECUENCIA.\n";
Coloreado en 0.002 segundos, usando GeSHi 1.0.8.4