[BioPython] Trying to transcribe RNA = error?

[Martin MOKREJŠ]

Martin MOKREJŠ wrote:

>>> What would be really meaningful if the object would be a genbank record
>>> with introns and transcribe() would rip out all the introns or even
>>> better allow me to choose which-ones to splice-out. Like to retain
>>> exons 1, 5, 7.
>>
>> I was thinking it would be nice to take a SeqRecord with SeqFeatures
>> (i.e. a parsed GenBank file) and ask Biopython to give you the
>> nucleotide sequence for the feature.  In the case of intron/exons
>> these are held using the subfeatures, so it would make sense to obey
>> those.  I see this as a variation on splicing/subsetting (in python
>> speak) a sequence.
> 
> Personally I find this the only reason to justify existence of such
> a function. It it can abstract me from the underlying format and give
--------------^^ If

> me a true, spliced transcript. Of course I would rather prefer to raise
> an exception if there would be more than one intron and I would have
> not specified which intron index-positions to use.
----------------------^^^^^^ exon

Sorry for the typos.