You can get some good information from the BBS site on routes, etc. but you may need to convert it to a more friendly format (i.e. a CSV file). Once you have your dataframes linked up, you can spatially subset based on some of the data that has already been assigned per BBS route…
BBS <- BBS[BBS$BCR %in% c(1,2), ]