Foro - Perl en Español

Buenas tardes. He conseguido que mi programa lea de tres en tres una secuencia y otro programa que me encuentre la secuencia ATG y me diga dónde está.

Ahora estoy tratando de que dada una secuencia, me encuentre la secuencia ATG y cuente de tres en tres a partir de ATG.

#!/usr/bin/perl
$secuencia = 'AAAATGCGGTTATGCCCATGGGTAAGTG';
$indica_orf = 0; #Indicada cuando encuentra un orf
$orf = "";
my @var = (); #La usamos para guardar los orfn
my $count = 0; #Contador de orf
my $posi = 0; #Posición inicial
my $posf = 0; #Posición final
for (my $i = 0; $i < length $secuencia; $i=$i+3) {
my $str ="";
$a = substr ($secuencia, $i, 3); #Para que me vaya leyendo las letras de 3 en 3.
if (($a eq 'ATG') and ($indica_orf == 0)) { #Si $a es igual a ATG y el indicador esta a 0
$indica_orf = 1;
$posi = $i+1;
}
if ($indica_orf == 1) { #Si estamos dentro de un orf
$str = ($a);
$orf = $orf.$str; #Sustituye las 3 letras de a por su valor y va concatenando
}
if ( $a =~ /TA[AG]|TGA/i ) { #Cuando a sea igual a TAA TAG o TGA
if ($indica_orf == 1) {
$indica_orf = 0; #Sal del orf
}
}

$orf = "";

Llegado aquí me he liado un poco y me da error. ¿Me podríais ayudar?
¡¡¡Gracias!!!

No sé a qué te refieres, yo quiero que me lea la secuencia a partir de ATG para ver si está en fase hasta que encuentre la secuencia de parada.

Es que dada la secuencia, intento que me cuente de tres en tres a partir de ATG hasta la secuencia de parada.

Pero me da error, y ya no sé qué más cambiarle ya que lo estoy intentando de varias formas.

Lo que quiero es que me muestre si está en fase o no.

Disculpe por las molestias y ante todo muchas gracias.

He diseñado un programa que dada un secuencia, encuentra el triplete de inicio ATG y el de parada.

#!/usr/bin/perl
use strict;
use warnings;

# Leer el fichero

my $sequence = 'AAATTTTTTATGCCCTTTATGGGGTTTTAGTGACCCGGGTG';
my $i = 0;
my @orfs = $sequence =~ /(ATG(?:.{3})*?)(TAA|TAG|TGA|.{1,3}$)/g ;

foreach my $orf (@orfs) {

{
my @tripletes = ( $orf =~ m/\w{3}/g );
my $nueva_secuencia = join "-", @tripletes;
print $orf. "\n\n";
$i++;
}
}
print $i. "\n";

En este caso me sale por el terminal:


que sería la secuencia empezando por ATG hasta llegar al codón de parada TAG. Me muestra también el codón de parada TAG.

Pero tengo un problema que no entiendo: ¿Qué significa el 2 que sale por el terminal?

Cuando cambio la secuencia y pongo 'AAATTTTTTATGCGCCTTTATGGGGTTTTAGTGACCCGGGTG y no está en fase me saldría esto:



por lo que ahora me sale GTG que no es el codón de parada, y me muestra la secuencia desde ATG hasta el codón final.

Yo quería que sólo me muestre la secuencia de tres en tres a partir de ATG hasta el codón de parada. ¿Qué me está fallando?

Muchas gracias.

El dos sé que sale del último print, solo que lo que yo quería que me mostrase no me lo mostraba y no entendía bien por qué me salía el dos, y cuando no, me daba 0 al modificar la secuencia.

Pero ante la duda lo quité ya que no llegaba a mi fin.

Gracias por la corrección.