Parsing HMMER Reads in Biopython

[maxbio]

Hello all,

I have obtained HMMER reads on 60 metatranscriptomic sequences, against Pfam, dbCAN2, and MEROPS databases. I ran the protein sequences through hmmsearch and now want to parse them for analysis in R. I want to use Biopythons SearchIO to do this since it has a parser for HMMER reads (HMMERIO) but the instructions on the Biopython wiki are very unclear as to how to do this. Does anyone have experience parsing HMMER outputs in Biopython?

Thank you very much.