code: def num_species(self):
          
          all_headers = reduce(lambda x, y: set(x) | set(y),
                               (rec.get_names() for rec in self.records))
          return len(all_headers)
docs: Returns the number of species found over all records