I used the protein files provided by bacmet database to build a custom database of args_oap, extracted the GI_number of the protein database fasta file according to the tips of github, and then used the annotation file to build structure.txt.I have since changed the column names to level1, level2, level3. There was no problem in running stageone. The problem was stagetwo. How to solve it?