Clustering metagenomic reads using a bidirectional Burrows-Wheeler transform