Discussion :

[Jasmine80]

Bonjour,

Le problème est assez simple mais je ne sais pas quelle est la meilleure façon de procéder.
J'ai deux séquences $seq et $seq_ali et je cherche à obtenir $seq_ali_modif.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
use strict;
use warnings;
 
my $seq = 'ATgcTGtC';
my $seq_ali = 'ATG..CTGT.C';
 
 
my $seq_ali_modif = 'ATg..cTGt.C';

Pourriez-vous me conseiller?


Merci;

[Jasmine80]

Voici la solution que j'ai trouvé mais cela ne me parait pas terrible, je me demande si il n'y a pas moyen de mieux faire avec des expressions régulières.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
use strict;
use warnings;
 
my $seq = 'ATgcTGtC';
my $seq_ali = 'ATG..CTGT.C';
my $seq_ali_modif = "";
 
my @a_seq_ali = split (//, $seq_ali);
my @a_seq = split (//, $seq);
 
my $j = 0;
 
foreach my $nuc (@a_seq_ali){
    if($nuc !~ /\./){
        print "$nuc\t$a_seq[$j]\n";
        if($a_seq[$j] !~ /$nuc/i){
            print "ERREUR\n";
        }
        else{
            $seq_ali_modif .= $a_seq[$j];
        }
        $j++;
    }
    else{
        $seq_ali_modif .= '.';
    }
}
 
print $seq_ali_modif ."\n";
# ATg..cTGt.C

Avec une expression régulière dans ce genre

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
my ($sous_seq, $reste) = ($seq_ali =~ /^([a-z]+)\.(.*)$/i);
# on récupère la première partie avant le point
if ($seq =~ /($sous_seq)(.*)$/i){
    $seq_ali_modif = $1.$reste;
    # on place les minuscules
}

Merci pour vos conseils,