docs: def write_output(output_report, sample_name, multi_positions, genus, percent_contam, contam_stddev, total_gene_length,
                   database_download_date, snp_cutoff=3, cgmlst=None):
      """
      Function that writes the output generated by ConFindr to a report file. Appends to a file that already exists,
      or creates the file if it doesn't already exist.
      :param output_report: Path to CSV output report file. Should have headers SampleName,Genus,NumContamSNVs,
      ContamStatus,PercentContam, and PercentContamStandardDeviation, in that order.
      :param sample_name: string - name of sample
      :param multi_positions: integer - number of positions that were found to have more than one base present.
      :param genus: string - The genus of your sample
      :param percent_contam: float - Estimated percentage contamination
      :param contam_stddev: float - Standard deviation of percentage contamination
      :param total_gene_length: integer - number of bases examined to make a contamination call.
      :param cgmlst: If None, means that rMLST database was used, so use rMLST snp cutoff. Otherwise, some sort of cgMLST
      database was used, so use a different cutoff.
      """
      # If the report file hasn't been created, make it, with appropriate header.
      if not os.path.isfile(output_report):
          with open(os.path.join(output_report), 'w') as f:
              f.write('Sample,Genus,NumContamSNVs,ContamStatus,PercentContam,PercentContamStandardDeviation,BasesExamined,DatabaseDownloadDate\n')
  
      if multi_positions >= snp_cutoff or len(genus.split(':')) > 1:
          contaminated = True
      else:
          contaminated = False
      with open(output_report, 'a+') as f:
          f.write('{samplename},{genus},{numcontamsnvs},'
                  '{contamstatus},{percent_contam},{contam_stddev},'
                  '{gene_length},{database_download_date}\n'.format(samplename=sample_name,
                                                                    genus=genus,
                                                                    numcontamsnvs=multi_positions,
                                                                    contamstatus=contaminated,
                                                                    percent_contam=percent_contam,
                                                                    contam_stddev=contam_stddev,
                                                                    gene_length=total_gene_length,
                                                                    database_download_date=database_download_date))