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 #!/usr/bin/perl -w use strict; use Bio::Restriction::Analysis; # get a DNA sequence from somewhere my $seq = Bio::PrimarySeq->new (-seq =>'CTAGCTTAATTCATTAGCTCTGACTGCAACGGGCAATATGTCTCAG', -is_circular => 1, # plasmide -primary_id => 'synopsis', -molecule => 'dna'); # now start an analysis. # this is using the default set of enzymes my $ra = Bio::Restriction::Analysis->new(-seq=>$seq, -is_circular => 1); # AluI is one them. Where does it cut? # This is will return an array of the sequence strings my $enz = 'AluI'; my @frags = $ra->fragments($enz); # how big are the fragments? print "AluI fragment lengths: ", join(' & ', map {length $_} @frags), "\n";
Le module Bio::Restriction::Analysis permet de voir le résultat d'une digestion avec une enzyme de restriction sur une séquence (qui est donc coupée). J'aimerais pouvoir spécifier que la séquence est circulaire.
Le module Bio::PrimarySeq (objet traité par le module précédent) crée des objets contenant la séquence et permet de définir si elle est linéaire ou circulaire mais apparemment cela ne suffit pas et Bio::Restriction::Analysis la traite comme linaire.
Le module Bio::Restriction::Analysis possède une 'méthodes interne' nommée _circular qui permet de traiter les séquences circulaires. Je ne m'y connais pas en OO, pourriez-vous me dire comment l'utiliser dans le script ci dessus?
Merci,
Partager