Hi, i m trying to align a fasta file(sequence.fasta) that contains 10 proteins.I want to use MUSCLE in order to do the alignment.I 've installed MUSCLE 3.6 and now im trying to have the right code.I m really new to all this stuff and my onle experience with perl is during the last 3 months so i need your help. MY Code: (or basically the code from bioperl.wiki)

#!/ usr/bin/ perl use strict; use warnings; use Bio::Tools::Run::Alignment::Muscle; # Build a muscle alignment factory my $factory =Bio::Tools::Run::Alignment::Muscle -> new(@params ); # Pass the factory a list of sequences to be aligned . my $inputfilename = 'sequence.fasta '; # $aln is a SimpleAlign object . my $aln = $factory -> align ( $inputfilename ); # Perform a profile alignment on a MSA to include more seqs my $alnfilename = 'sequence.fasta '; my $seqsfilename = 'outfile '; $aln = $factory -> profile ( $alnfilename , $seqsfilename );
After trying to execute it, i get the message " Global symbol "@params" requires explicit package name". How do i define @params? I mean wthats the code thats missing? My task is to make a phylogentic tree from this fasta file, but the first step is to make the alignhment.The howtos don't really help me...

In reply to Aligning protein sequences using MUSCLE through Perl by DimitriC

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":



  • Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
  • Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
  • Read Where should I post X? if you're not absolutely sure you're posting in the right place.
  • Please read these before you post! —
  • Posts may use any of the Perl Monks Approved HTML tags:
    a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
  • You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
            For:     Use:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.