[Bioperl-l] sh: 1: Syntax error: "(" unexpected
Antony03
antony.vincent.1 at ulaval.ca
Fri Mar 8 22:56:59 EST 2013
Hi!
I'm trying to use emboss program (water or needle) with a bioperl script. I
got this error when I try the script: sh: 1: Syntax error: "(" unexpected
Here is the script:
#!/usr/bin/perl
# get an EMBOSS factory
use Bio::Factory::EMBOSS;
use Bio::AlignIO;
use Bio::SeqIO;
$f = Bio::Factory::EMBOSS -> new();
# get an EMBOSS application object from the factory
$water = $f->program('water');
$inputfilename = "krt1.fasta";
$seq_to_test = Bio::SeqIO->new(-file => $inputfilename ,
'-format' => 'fasta');
$inputfilename2 = "krt2.fasta";
@seqs_to_check = Bio::SeqIO->new(-file => $inputfilename2 ,
'-format' => 'fasta');
# here is an example of running the application -
# water can compare 1 sequence against 1 or more sequences
# in a database using Smith-Waterman
#this would be a list of seqs to compare
# (could be just 1)
my $wateroutfile = 'out.water';
$water->run({-sequences => $seq_to_test,
-seqall => \@seqs_to_check,
-gapopen => '10.0',
-gapextend => '0.5',
-outfile => $wateroutfile});
# now you might want to get the alignment
my $alnin = Bio::AlignIO->new(-format => 'emboss',
-file => $wateroutfile);
while ( my $aln = $alnin->next_aln ) {
# process the alignment -- these will be Bio::SimpleAlign objects
}
Have you an idea of what can do this error?
Thanks!
--
View this message in context: http://old.nabble.com/sh%3A-1%3A-Syntax-error%3A-%22%28%22-unexpected-tp35154292p35154292.html
Sent from the Perl - Bioperl-L mailing list archive at Nabble.com.
More information about the Bioperl-l
mailing list