docs: def estimate_percent_contamination(contamination_report_file):
      """
      Estimates the percent contamination of a sample (and standard deviation).
      :param contamination_report_file: File created by read_contig,
      :return: Estimated percent contamination and standard deviation.
      """
      contam_levels = list()
      with open(contamination_report_file) as csvfile:
          reader = csv.DictReader(csvfile)
          for row in reader:
              lowest_count = 99999
              base_counts = row['Bases'].split(';')
              for count in base_counts:
                  num_bases = int(count.split(':')[1])
                  if num_bases < lowest_count:
                      lowest_count = num_bases
              total_coverage = int(row['Coverage'])
              contam_levels.append(lowest_count*100/total_coverage)
      return '%.2f' % (np.mean(contam_levels)), '%.2f' % np.std(contam_levels)