[Bioperl-l] problem parsing MEME file

Stefan Kirov skirov at utk.edu
Tue Sep 21 07:08:06 EDT 2004


I had a similar problem in the past and never got an answer from the guy 
who wrote the parser. There is also Bio::Matrix::PSM::IO::meme, but it 
works only on DNA.
James Thompson  has offered to write a psi-blast parser and will create 
objects that will also be useful for a protein meme parser. Then it will 
be easy to add this functionality to Bio::Matrix::PSM::IO::meme.
Stefan

Diego Riano wrote:

>Hi all
>
>I have written a small script to parse MEME results, but I am getting
>and error.  I am using Bioperl 1.4 and MEME 3.0.10, it should work but
>if does not.  could anyone help me?
>
>Here is the script:
>
>#!/usr/bin/perl
>
>use strict;
>use warnings;
>
>use Bio::AlignIO;
>
>my $inputf = "$ARGV[0]";
>my $in  = Bio::AlignIO->new(-file => $inputf,
>                            -format => 'meme');
>
>
>while ( my $aln = $in->next_aln() ) {
>        print "X\n";
>    }
>
>This is the error i got:
>
>-------------------- WARNING ---------------------
>MSG: Unrecognized format:
>6_04024                     163 5.87e-128 DYSCGGDDSG
>TEEEQAAFMKELERFYREKLMEFKPPKFYGEGLNCLKLWRQVTGLGGYDQVTSCKLWRQVGESFKPPKTCTTVSWTFRNFYEKALLEYEKHKIETGEFQV AASNLTERIT
> 
>---------------------------------------------------
>
>thanks
>
>Diego Riano
>  
>



More information about the Bioperl-l mailing list