[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