code: def estimate_genome_size(self):
          
          logging.info('Estimating genome size using kmercountexact')
          for sample in self.metadata:
              # Initialise the name of the output file
              sample[self.analysistype].peaksfile = os.path.join(sample[self.analysistype].outputdir, 'peaks.txt')
              # Run the kmer counting command
              out, err, cmd = bbtools.kmercountexact(forward_in=sorted(sample.general.fastqfiles)[0],
                                                     peaks=sample[self.analysistype].peaksfile,
                                                     returncmd=True,
                                                     threads=self.cpus)
              # Set the command in the object
              sample[self.analysistype].kmercountexactcmd = cmd
              # Extract the genome size from the peaks file
              sample[self.analysistype].genomesize = bbtools.genome_size(sample[self.analysistype].peaksfile)
              write_to_logfile(out, err, self.logfile, sample.general.logout, sample.general.logerr, None, None)
docs: Use kmercountexact from the bbmap suite of tools to estimate the size of the genome