This script performs taxonomic annotation, filtering, and alpha/beta diversity analysis of metagenomic sequencing data obtained from water column and sediment samples collected along a pollution gradient in the Virilla River, Costa Rica. Sampling was performed during the dry and rainy seasons.