The cookbook example of opening and reading a Genbank file can easily be translated into Jython:
#!/usr/bin/env jython
import sys
from java.io import *
from java.util import *
from org.biojava.bio import *
from org.biojava.bio.seq.db import *
from org.biojava.bio.seq.io import *
from org.biojava.bio.symbol import *
if __name__ == "__main__":
br = BufferedReader( FileReader( sys.argv[1] ) )
sequences = SeqIOTools.readGenbank(br)
while sequences.hasNext():
seq = sequences.nextSequence()
print seq.seqString()
No comments:
Post a Comment