[Bioperl-l] assistance with bioperl
Mark A. Jensen
maj at fortinbras.us
Wed Sep 30 18:24:59 UTC 2009
Antonina--
Try the following:
Make sure that cysprot1a.msf and cysprot1b.fa are in the current directory,
or use full path names for the files.
MAJ
----- Original Message -----
From: "Antonina Iagovitina" <antonina.iagovitina at epfl.ch>
To: <bioperl-l at lists.open-bio.org>
Sent: Wednesday, September 30, 2009 2:09 PM
Subject: [Bioperl-l] assistance with bioperl
> Here is the error message I get when I try to align a sequence to an existing
> alignment. Please help
> I am using Windows XP and Clustalw version1.83
>
> MSG:
> ERROR: Could not open sequence file (-profile)
> No. of seqs. read = -1. No alignment!
>
> use Bio::AlignIO;
> use Bio::SeqIO;
> use Bio::Seq;
> use Bio::Tools::Run::Alignment::Clustalw;
>
> my @params = ('ktuple' => 2, 'matrix' => 'BLOSUM');
> my $factory = Bio::Tools::Run::Alignment::Clustalw->new(@params);
> $str = Bio::AlignIO->new(-file=> 'cysprot1a.msf');
> $aln = $str->next_aln();
> $str1 = Bio::SeqIO->new(-file=> 'cysprot1b.fa');
> $seq = $str1->next_seq();
> $aln = $factory->profile_align($aln,$seq);
> end
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
>
More information about the Bioperl-l
mailing list