[Bioperl-l] bug #900: Bio::SeqIO::swiss.pm - gene names

Paul Gordon gordonp@niji.imb.nrc.ca
Thu, 1 Mar 2001 07:24:24 -0400 (AST)


> Done. It would be nice if someone could jump in and dig up a
> real-world test case, i.e., a swissprot entry that indeed has
> multiple gene names. Ideally we can add a test to SeqIO.t ... (I
> relied on / OR / being the delimiter, as given in the bug report.)

Go nuts! :-)

% time perl -ne 'BEGIN{$/="//\n"}print "$1\n" if /\nGN.*?\nGN/s && /^ID\s+(\S+)/' sprot.dat
BRC1_DROME
BRC4_DROME
CALM_HUMAN
EF12_SCHPO
ERF2_YEAST
GLS1_YEAST
H3_CAEEL
H4_SCHPO
MAT1_YEAST
MAT2_YEAST
MS31_CAEEL
OR2A_DROME
R10A_YEAST
RDPO_SCHPO
RL11_YEAST
RL12_YEAST
RL18_YEAST
RL20_YEAST
RL23_YEAST
RL31_YEAST
RL35_YEAST
RL41_YEAST
RL44_YEAST
RLK5_SCHPO
RS11_YEAST
RS21_YEAST
RS22_YEAST
RS23_YEAST
RS25_YEAST
RS28_YEAST
RS6_YEAST
THTR_MYCTU
TIS5_ECOLI
TRA9_MYCTU
TUP1_YEAST
UBIQ_YEAST
VPS4_YEAST
YI11_HALN1
YIA4_MYCTU
ZDS1_YEAST
11.96u 3.87s 0:16.25 97.4%