[Bioperl-l] It seems a bug: Bio/SeqFeatureI.pm?

Sang Chul Choi goshng at naver.com
Fri Oct 15 00:28:45 EDT 2004


Dear Friends:

In the file Bio/SeqFeatureI.pm, there seems a small bug which should be fixed.
I tried to used the method of Seq Object called spliced_seq(). But, when I
have Bioperl parse the line like 

complement(join(2352..2759,2809..3015))

It just complements each fragment and then joins them together. And, I have
a strange translated protein. So, I changed the order of operation, and
it works.

            # This is dumb, subseq should work on locations...
            if( $loc->strand == 1 ) {
                $seqstr .= $called_seq->subseq($loc->start,$loc->end);
            } else {
                #$seqstr = $called_seq->trunc($loc->start,$loc->end)->revcom->seq() . $seqstr;
                $seqstr .= $called_seq->trunc($loc->start,$loc->end)->revcom->seq();
            }

Thanks,

Sang Chul



More information about the Bioperl-l mailing list