[Bioperl-l] bug in location->to_FTstring

Philip Lijnzaad p.lijnzaad@med.uu.nl
Wed, 04 Dec 2002 15:41:48 +0100


>> mmm ... this doesn't appear to be part of BioPerl 1.0.2. Is it on CVS trunk
>> or branch? Cheers,
> Yes - only CVS live code.

OK. Incentally, it does appear that your out-commented code _is_ needed, at
least sometimes: if I uncomment your comments (like I needed to in yesterday's
exercise) and then use SeqIO to do a read and write as

  my $in  = Bio::SeqIO->newFh(-fh => \*STDIN , '-format' => 'EMBL');
  my $out = Bio::SeqIO->newFh(-fh=> \*STDOUT, '-format' => 'EMBL');
  print $out $_ while <$in>;

I get differences like:

< FT   CDS             join(complement(173191..173432), ...

> FT   CDS             complement(join(complement(173191..173432), ... 

If I reintroduce your comments, things work fine again. So something is fishy
... bug #1074 is still at large. Cheers,


                                                                      Philip
-- 
Philip Lijnzaad
Genomics Laboratory
Dept. of Biomedical Genetics
University Medical Center (UMC), Utrecht
Stratenum, STR 3.223
room 3.126 (not in on Thursdays)
P.O. Box 85060, 3508 AB Utrecht
(Universiteitsweg 100, 3584 CG Utrecht)
The Netherlands
tel: +31 (0)30 253 8464
fax: +31 (0)30 253 9035