1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
#!/usr/bin/env python
"""Display the SeqFeatures produced by the parser.
This produces a ton of output so it is possible to hand check what is
produced by the parser with the original GenBank file to make sure
everything is being parsed and output properly.
Usage:
./check_output.py <GenBank file to parse>
"""
# standard library
from __future__ import print_function
import sys
# GenBank stuff to test
from Bio import GenBank
if len(sys.argv) != 2:
print("Usage ./check_output.py <GenBank file to parse>")
sys.exit()
parser = GenBank.FeatureParser(debug_level=2)
with open(sys.argv[1], 'r') as handle:
iterator = GenBank.Iterator(handle, parser)
while True:
cur_record = next(iterator)
if not cur_record:
break
print("***Record")
print("Seq: %s" % cur_record.seq)
print("Id: %s" % cur_record.id)
print("Name: %s" % cur_record.name)
print("Description: %s" % cur_record.description)
print("Annotations****")
for annotation_key in cur_record.annotations:
if annotation_key != 'references':
print("Key: %s" % annotation_key)
print("Value: %s" % cur_record.annotations[annotation_key])
else:
print("References*")
for reference in cur_record.annotations[annotation_key]:
print(str(reference))
print("Feaures")
for feature in cur_record.features:
print(feature)
|
