[BioPython] Downloading CDS sequences
Animesh Agrawal
animesh.agrawal at anu.edu.au
Wed Jan 21 07:14:42 UTC 2009
>just return a dummy SeqRecord with an empty
>sequence for the "coded_by" cases, rather than going online wasting
>the NCBI server time.
Ok. So that's the reason and I was wondering why you want to return empty
sequence.
-----Original Message-----
From: p.j.a.cock at googlemail.com [mailto:p.j.a.cock at googlemail.com] On Behalf
Of Peter
Sent: Tuesday, 20 January 2009 8:44 PM
To: Animesh Agrawal
Cc: BioPython Mailing List
Subject: Re: [BioPython] Downloading CDS sequences
On Tue, Jan 20, 2009 at 8:31 AM, Animesh Agrawal
<animesh.agrawal at anu.edu.au> wrote:
>
> Peter,
> Thanks a lot. ...
> I tested your new script for Downloading CDS sequences. It was working
fine
> for records like P24673 but couldn't get it working for precise
"coded_by"
> string situation unless I comment(#return SeqRecord(Seq(""))) statement in
> get_nuc_record() function. I don't understand why?
That was a deliberate mistake to test your understanding (joke). I'm
pleased you worked out what was wrong!
I put that line in to speed up my testing - and forgot to remove it.
Basically that line was to just return a dummy SeqRecord with an empty
sequence for the "coded_by" cases, rather than going online wasting
the NCBI server time. With hindsight I could have edited my example
GenBank file to focus on the cases of interest. Does that make sense?
Sorry for the confusion,
Peter
More information about the Biopython
mailing list