[Bioperl-l] Bio/SeqIO/swiss.pm parsing error -

Erik er at xs4all.nl
Fri Nov 3 20:56:47 UTC 2006


I should have mentioned two more changes needed
in Bio/SeqIO/swiss.pm for the new Swissprot
format:

In Bio/SeqIO/swiss.pm, next_seq():
===============
  $params{'-namespace'} = ($2 eq 'Reviewed'   ) ? 'Swiss-Prot' :
                          ($2 eq 'Unreviewed' ) ? 'TrEMBL'     :
                       #  ($2 eq 'STANDARD'   ) ? 'Swiss-Prot' :
                       #  ($2 eq 'PRELIMINARY') ? 'TrEMBL'     :
                           $2;
===============


and in Bio/SeqIO/swiss.pm, write_seq():
===============
   $div = ($ns eq 'Swiss-Prot') ? 'Reviewed'    :
          ($ns eq 'TrEMBL')     ? 'Unreviewed'  :
       #  ($ns eq 'Swiss-Prot') ? 'STANDARD'    :
       #  ($ns eq 'TrEMBL')     ? 'PRELIMINARY' :
           $ns;
===============

above, the old lines are shown as outcommented.


again, hth,

Erik





More information about the Bioperl-l mailing list