How to calculate adjusted P-value (FDR, BH) from a table of P-values and select below a cut off (q value <0.25)

Hi friends,
I have a table of P-values like:

Feature X type LDA P-value
Prevotella_copri 5.17933651645 test 4.60356342912 0.0922867588261
Butyrivibrio_crossotus 4.69683479993 control 4.07248710826 0.735615882525
Bacteroides_dorei 4.60167732782 test 3.96437649555 0.843672782494
Akkermansia_muciniphila 4.28836733193 control 3.90356518613 2.98756798706E-05
Bacteroides_cellulosilyticus 4.27495706056 control 3.89732714054 0.389612546234
Bacteroides_intestinalis 4.19463599386 control 3.8851611951 0.0312193961172
Bacteroides_eggerthii 4.27603800675 test 3.87752454514 0.086688300821
Alistipes_putredinis 4.6767463495 control 3.86667535783 0.0419820771958
Ruminococcus_bromii 4.46456513989 test 3.808952528 0.225163911724
Bacteroides_uniformis 4.85724578241 control 3.80355808508 0.274845078066
Eubacterium_sp_CAG_180 4.28515167506 test 3.73437627148 0.0137754984604
Bacteroides_vulgatus 4.80795432151 control 3.72498343619 0.978424756919
Paraprevotella_xylaniphila 4.1622177234 test 3.71112597761 0.239742113249
Eubacterium_rectale 4.41796600305 test 3.70791598112 0.306888374573

I want to calculate the adjusted P-value (False Discovery Rate, Benjamini-Hochberg) and make an additional column in the output. I also expect one additional file where only q-value <0.25 will be featured.
Can anyone please tell me how can I do that?

Thanks

Check the function p.adjust(), especially its method parameter.