Substitution Error in pal2nal perl script
0
0
Entering edit mode
7.2 years ago
apulunuj ▴ 30

I need some help in understanding how to work, the pal2nal.pl script. The program can be obtained from the website, http://www.bork.embl.de/pal2nal/.

Whenever I run the script I get an error message at the line, (the line of interest has a comment beside it). THe error message says 'Substitution loop at line 338, which is the line of interest'. I don't see why the script should work.

The part of the script with the line:

undef(@nucid);

undef(%id2nucseq);

$nseq = -1;

foreach $i (0..$#nucfiles) {

   open(NUCFILE, "< $nucfiles[$i]") || die "Can't open $nucfiles[$i]"; ##Line of Interest

    $nucfiledata = <NUCFILE>;

    close(NUCFILE);

    $nucfiledata =~ s/\x0D\x0A|\x0D|\x0A/\n/g;

    foreach (split(/\n/, $nucfiledata)) {

        if (!/^#/ && /\S+/) {

            if (/^>(\S+)/) {

                ++$nseq;

                $tmpnucid = $1;

                push(@nucid, $tmpnucid);

            } else {

                s/[^a-zA-Z]//g;

                $id2nucseq{$tmpnucid} .= $_;

            }

        }

    }

}
alignment nucleotide seq perl script • 1.7k views
ADD COMMENT
0
Entering edit mode

It seems you are trying pal2nal on a massive file - see here. What is it? What is your version of perl?

ADD REPLY
0
Entering edit mode

perl 5.10.1 is what I am using/

ADD REPLY

Login before adding your answer.

Traffic: 2723 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6