julio_514 has asked for the wisdom of the Perl Monks concerning the following question:
#c:\perl\perl.exe use strict; use warnings; use Bio::Seq; use Bio::Tools::Run::Alignment::Clustalw; BEGIN { $ENV{CLUSTALDIR} = 'c:/CLUSTALW2/' } my $dna_seq = "NACCNANCGGGGGNNAAANTNNNACACTNCNGGGGGGNTNNTTANTGNGTNACAC +ANNNGNCTNNNGGNNNCNANANTCNACAG +GACTTAGAAGNCTNCANNGGANCATNCCCTGCTANACN +ANGNNANCACTCTNCAGNNCNNANCNTNGGGCNNCCTNNTCNNTNNNCCATNNGNGCTNC +ANNANCN +TANANNGATNANTANAGNANNTTNTGTTATGNNGNGTG"; my$tn5_seq = "GGGGNCCCGCTTACGATACTTCGTTATAGCATACATTATACGAAGTTATCAGATCC +CCCTGGATGGAAAACGGGAAAGGTTCCG +TCCAGGACGCTACTTGTGTATAAGAGTCAGGTT"; $dna_seq =~ s/\s+//g; $tn5_seq =~ s/\s+//g; my $id_1 = 'test1_DNA'; my $id_2 = 'test2_DNA'; my $alphabet = 'DNA'; my $seq_obj = Bio::Seq->new(-seq => $dna_seq, -display_id => $id_1, -alphabet => $alphabet); my $tn5_obj = Bio::Seq->new(-seq => $tn5_seq, -display_id => $id_2, -alphabet => $alphabet); my @params = ('ktuple' => 2, 'matrix' => 'BLOSUM'); my $factory = Bio::Tools::Run::Alignment::Clustalw->new(@params); my @seq_array =(); push @seq_array, $seq_obj, $tn5_obj; my $seq_array_ref = \@seq_array; my $aln = $factory->align($seq_array_ref);
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re: bioperl and clustalw help
by Anonymous Monk on Oct 23, 2008 at 09:41 UTC | |
by julio_514 (Acolyte) on Oct 23, 2008 at 13:02 UTC | |
by Anonymous Monk on Feb 09, 2015 at 12:49 UTC | |
by Corion (Patriarch) on Feb 09, 2015 at 12:51 UTC |