I have the fallowing code(see below) and it is giving me this error

Use of /c modifier is meaningless without /g at hwk5.pl line 31. Barew +ord found where operator expected at hwk5.pl line 31, near "-database + => '/class" (Might be a runaway multi-line // string starting on line 4)(Missing + operator before lass?) Unquoted string "lass" may clash with future reserved word at hwk5.pl +line 31. syntax error at hwk5.pl line 31, near "-database=> '/class" Execution of hwk5.pl aborted due to compilation errors.

Here is the code and I will mark line 31 by ********

use Bio::SeqIO; use Bio::SearchIO; use Bio::Tools::Run::StandAloneBlast; use strict print "Write the file name from where you are getting the sequence form\n\n"; $file=<STDIN>; chomp $file; print "Write put the name of the output file\n\n"; $output=<STDIN>; chomp $output; #declare the required object my $seqio= ''; my $seqobj= ''; my $factory = ''; my $blast_report= ''; my $out_file= ''; $seqio= Bio::SeqIO->new(-file => $file, -format=> 'fasta'); while ($seqobj = $seqio->next_seq()) { $factory = Bio::Tools::Run::StandAloneBlast->new( -program =>'blastn', ***************-database => '/class/bi617a/share/blast/data/rat.rna.fn +a', -e=> 1e-6, _READMETHOD=> "Blast" ); $out_file = $seqobj->display_id . $outfile; $factory->outfile($out_file); $blast_report = $factory->blastall($seqobj); } #the blast output name my $filename= $outfile; #create an object tot he blast output file my $input= Bio::SearchIO->(-format => 'blastn', -file => $filename); #Print out the query, hit names and percent indenty, percent query ali +gned print "QueryName\tQueryLenght\tLengthMatch\tAccesionNumber\tDescription\tRan +k\tScore\%ID\tStartPostQ\tEndPosQ\tStartPosH\tEndPosH\n"; #Go thru the BLAST reports one by one while (my $result = $report->next_result) { while($hit=$result->next_hit) {#Go through eadch HSP for this sequence while ($hsp= $hit->next_hsp) { #Print some tab-delimited data about this HSP print $result->query_name, "\t"; print $result->query_length, "\t"; print $hsp->hsp_length, "\t"; print $result->query_accession, "\t"; print $hit->description, "\t"; print $hsp->rank, "\t"; print $hsp->score, "\t"; print $hsp->percent_identity, "\t"; print $hsp->start('query'), "\t"; print $hsp->end('query'), "\t"; print $hsp->start('hit'), "\t"; print $hsp->end('hit'), "\t"; } } }

In reply to Bioperl FASTA parsing by rolandomantilla

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.