We conducted metagenomics NGS sequencing and then used the args_oap pipeline to detect antibiotic resistance genes in both control and trait groups. How should we interpret the output from args_oap for comparing these groups? And besides args_oap, are there any recommended alternative or improved approaches (bioinformatics tools) for this type of analysis?