[Bioperl-l] How to check mRNA moltype?
Wed, 17 Oct 2001 11:37:32 +0100
moltype() is the correct method. It works on EMBL sequences so it should on
GenBank, too. I had a look at the code and RNA option was completely
I rewrote the lines doing it and committed them to bioperl-live. I hope you
are using CVS and can see the changes immediately.
This fix should go into 07 branch, too, but first I'd like to get feedback
that the fix really works on larger data sets. I tried it only on t/data
Does anyone have a system handy to try this? The potential problem is that
it is DNA/mRNA/... keyword is allowed to be missing completely from the
Incidentally, the LOCUS line is going to change slightly in GenBank
distributions in December. I tried the current parser on two examples
provided and it seems to keep on working properly. See
Hong Lu wrote:
> I want to differentiate mRNA sequence from other Sequences. I found $seqobj->moltype will always give me 'dna' value even if the sequence is mRNA type on the first line of genbank record.
> How to get the mRNA value fo error checking purpose?
> Henry Lu
> Bioinformatics Scientist
> Epigenesis Pharmaceuticals
> Find the best deals on the web at AltaVista Shopping!
> Bioperl-l mailing list
______ _/ _/_____________________________________________________
_/ _/ http://www.ebi.ac.uk/mutations/
_/ _/ _/ Heikki Lehvaslaiho email@example.com
_/_/_/_/_/ EMBL Outstation, European Bioinformatics Institute
_/ _/ _/ Wellcome Trust Genome Campus, Hinxton
_/ _/ _/ Cambs. CB10 1SD, United Kingdom
_/ Phone: +44 (0)1223 494 644 FAX: +44 (0)1223 494 468