Discussion :

[crochepatte]

Bonjour,
Je souhaiterai dans un sequence, remplacer les motifs correspondant à une expression regulière par le meme nombre d'un autre caractere (N par exemple)....

Exemple: avec une expression reguliere

Code :

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

1
2
3
 
 
s/([ATGC]{1})\1{4,}/ ... /g

Je souhaiterai que

"ATGCAAAAAAATGCTCG" devienne "ATGCNNNNNNTGCTCG"

Je pense que c'est assez facile, mais je n'ai pas trouvé l'astuce

Merci

[MarneusCalgarXP]

Code :

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

s/A/N/g

[GLDavid]

Salut ami bioinformaticien !

Voici ma méthode :

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
 
use strict;
use warnings;
 
my $s = "ATGCAAAAAAATGCTCG";
print $s."\n";
if($s=~/(A{4,}|C{4,}|G{4,}|T{4,})/){
        my $motif = $1;
        my $substitute;
        for(0..(length($motif)-1)){
                $substitute.='N';
        }
        $s=~s/$motif/$substitute/;
}
print $s."\n";

Apparemment, cela marche.
Et en plus, c'est facilement adaptable en fonction

En espérant t'avoir aidé.

@++

[Leishmaniose]

Salut,

Tout d'abord je me demande pourquoi ne pas avoir posté sur le sous-forum Perl > bioinformatique qui serait plus adapté pour ta requête. Ensuite, d'après ce que j'ai compris de ta question, il s'agit de masquer les répétition dans une séquence ADN... Pour cela, il existe déjà plusieurs programmes disponibles comme par exemple : repeatmasker, etc. sans compter les filtres utilisés par des programmes tel que blast et qui permettent de masquer les répétitions tout en cherchant les similitudes par alignement de séquences, notes que ces filtres sont utilisables seuls comme c'est le cas de dust. Mais bon si ça te dit de traduire ça en Perl...

Bonne continuation et @+ sur Forum des développeurs > Autres langages et outils > Perl > Bioinformatique

[crochepatte]

Merci beaucoup,

MarneusCalgarXP -> Pas si simple, je cherche pas à remplacer seulement des motifs d'une paire de base...

GLDavid-> J'avais deja une solution similaire mais je me demandais si la fonction substitute ne prenait pas en compte ce genre de requetes...Tant pis, je reste à ma bonne vieille méthode...

Leishmaniose-> Tout depends de la notion de "repeat" (on va gardé ce terme puisque tu l'as employé). La notion même de repeat au sens où je la considere dans mon application n'est pas du tout identique à celle considérée par les principaux outils disponibles via le web...Et combien même aprés avoir trouver un outil disponible sur le Web, celui doit etre disponible gratuitement, et integrable dans mon logiciel...Par conséquent, etant donné toutes ses contraintes, je trouve que "developper" une fonction de 10 lignes au maximum, est bien moins couteux et forcement parfaitement adapté.
Concernant mon post dans cette section, j'ai bien posté ici consciemment, car à mon sens, ma question ne se limite pas au domaine de la bioinformatique, mais concerne l'utilisation générale de la fonction "substitute"... J'espere avoir été clair

Donc à bientot sur le Forum des développeurs > Autres langages et outils > Perl > je choisis la section en fonction de mes demandes...

[Leishmaniose]

Dans ce cas est-il possible d'avoir une idée sur ta notion de repeat...

et pour ce qui est des utilitaires que j'ai cité dans mon message, je te confirme que ce sont des solutions totalement free et qui s'adaptent bien au développement d'autres logiciels... Par exemple le dust est écrit en C, son code source est disponible ainsi que plusieurs binary adaptés à une multitude de plate-forme...

Pour finir, je respecte parfaitement ta liberté de choisr le sous forum qui te parait le plus adapté pour ton post, et dans ce cas précis je comprend que ton post concerne l'utilisation globale de la fonction "substitute" d'où le choix du sous-forum langage. Néanmoins, je persiste à dire que la solution à ton problème ne peut être que bioinformatique... Je m'explique : As tu pensé à faire un tour sur le module BioPerl...! je suis sûr que tu trouveras une fonction toute prête pour ce genre d'application...

Bonne continuation.