May want to put this function within either helper_functions.R or dada2_pipeline.R. I'm leaning towards it being in the helper_functions.R file as this step is extra and not part of the main pipeline for DADA2.

http://benjjneb.github.io/dada2/tutorial.html#bonus-handoff-to-phyloseq