[BioPython] Trying to transcribe RNA = error?

Peter biopython at maubp.freeserve.co.uk
Tue Jul 15 21:15:15 UTC 2008


Now that the latest release is out, I wanted to return to a point
raised by Sebastian Bassi last month:

Do people agree that the Bio.Seq.transcibe() method should reject RNA sequences?

We already reject Seq objects with a protein alphabet.  I want to
accept Seq objects with a generic or nucleotide alphabet, and
sequences with a DNA alphabet are of course fine, but I had been on
the fence about trying to (re)transcribe RNA.

Unless anyone speaks up, I will side with Sebastian and add a check to
the  Bio.Seq.transcibe() method to reject RNA sequences, and a similar
check on the Bio.Seq.back_transcibe() method to reject DNA sequences.

Peter



More information about the Biopython mailing list