• Publicidad

Mostrar posición de codón

Perl aplicado a la bioinformática

Mostrar posición de codón

Notapor perl biología » 2015-06-12 18:34 @815

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.
perl biología
Perlero nuevo
Perlero nuevo
 
Mensajes: 12
Registrado: 2015-06-03 09:53 @453

Publicidad

Re: Mostrar posición de codón

Notapor explorer » 2015-06-12 18:44 @822

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
JF^D Perl programming & Raku programming. Grupo en Telegram: https://t.me/Perl_ES
Avatar de Usuario
explorer
Administrador
Administrador
 
Mensajes: 14486
Registrado: 2005-07-24 18:12 @800
Ubicación: Valladolid, España


Volver a Bioinformática

¿Quién está conectado?

Usuarios navegando por este Foro: No hay usuarios registrados visitando el Foro y 0 invitados

cron