This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Biopython SeqIO.parse not functioning for all entries in a seq

I'm working with some genbank seq files and have the following code:

for seq_record in SeqIO.parse("datafile_location, "genbank"):

And while it can run through most of the seqs in the seq file (which contains multiple seqs) I get the following error. Any thoughts about how to fix this?

Maybe delete the offending seq? It gets to record 92126 of 93145 and then throws the error.

I have tried re-downloading the seq file, but that doesn't fix the problem.

File "C:\python38\lib\site-packages\Bio\GenBank\Scanner.py", line 516, in parse_records record = self.parse(handle, do_features) File

C:\python38\lib\site-packages\Bio\GenBank\Scanner.py", line 499, in parse if self.feed(handle, consumer, do_features): File

"C:\python38\lib\site-packages\Bio\GenBank\Scanner.py", line 466, in feed self._feed_header_lines(consumer, self.parse_header()) File

"C:\python38\lib\site-packages\Bio\GenBank\Scanner.py", line 1801, in feed_header_lines previous_value_line = structured_comment_dict[ KeyError: 'Assembly-Data'

seqio biopython seqio.parse genbank files
• 2,198 views • link
updated by Biostar Biostar • written by mvhanson • mvhanson

Can you post your code here? There is a " missing in the for seq_record in SeqIO.parse("datafile_location, "genbank"): loop but I guess that's a typo in the post as opposed to the code?

• 0 views • link
written by Sej Modha Sej Modha

0 answers

No answers yet.

Log in to answer this question.

More posts like this