>>> from Bio import SeqIO
>>> files = ["gbvrli.seq"]
>>> gb_vrl = SeqIO.index_db("gbvrl.idx", files, "genbank")
>>> print "%i sequences indexed" % len(gb_vrl)
958086 sequences indexed

>>> print gb_vrl["GQ333173.1"].descriptionHIV-1 isolate F12279A1 from Uganda gag protein (gag) gene, partial cds.

>>> print gb_vrl["GQ333173.1"].description
HIV-1 isolate F12279A1 from Uganda gag protein (gag) gene, partial cds.

>>> print gb_vrl.get_raw("GQ333173.1")
LOCUS GQ333173 459 bp DNA linear VRL 21-OCT-2009
DEFINITION HIV-1 isolate F12279A1 from Uganda gag protein (gag) gene, partial
cds.
ACCESSION GQ333173
...
//