[Bioperl-l] Sim4 output parsing using Bioperl

Arnaud Kerhornou axk at sanger.ac.uk
Tue Jul 1 11:45:37 EDT 2003


Hi

I'd like to use Bioperl to parse Sim4 output.
I run sim4 using this command, "sim4 genomicDNA.fa cDNAs.fa P=1 A=4"

I can't get the parsing going through the all sim4 output, it stops on 
this type of result:

------------------>
seq1 = DICTY6P4_0001.gfseq.masked, 145194 bp
seq2 = CONTIGS.fa (>AU275246), 258 bp

 >DICTY6P4_0001
 >AU275246  FC-IC Dictyostelium discoideum cDNA clone FC-IC1786, mRNA 
sequence

seq1 = DICTY6P4_0001.gfseq.masked, 145194 bp
<------------------

There is a "last" statement in the method "parse_next_alignment" in 
Bio::Tools::Sim4::Results.pm

 /^seq1/ && do {
       if($started) {
           $self->_pushback($_);
           last;
       }
...
}

Why this last statement ?

Thanks
Arnaud



More information about the Bioperl-l mailing list