I have tried to make rarefaction curves using iNEXT.
ggiNEXT(out, type=1, facet.var="order") +
facet_wrap(~order, scales="free") +
theme_bw(base_size = 18) +
theme(legend.position="right") +
theme_minimal() +
scale_y_continuous(name="ASVs") +
scale_x_continuous(name="Number of Reads")