Hello!

I have gene sequences (or parts of genes) of a few unknown species from 2 bacteria communities. I would like to BLAST on JGI with these 2 communities and choose the ones with E value match close to 0.0 and get the scaffold number and phylogeny.

This is an example of what it looks like when I do it on JGI when I do this manually:

How can I automate this process using Python or Java?

Thank you!