code: def report(self, align_bam, ref_file, gtf_file, is_paired=False, rrna_file="null"):
          
  
          # collect duplication metrics
          dup_metrics = self._get_current_dup_metrics(align_bam)
          align_metrics = self._collect_align_metrics(align_bam, ref_file)
          insert_graph, insert_metrics = (None, None)
          if is_paired:
              insert_graph, insert_metrics = self._insert_sizes(align_bam)
  
          rnaseq_metrics = self._rnaseq_metrics(align_bam, gtf_file, rrna_file)
  
          summary_info = self._parser.get_summary_metrics(align_metrics,
                                                  dup_metrics,
                                                  insert_metrics=insert_metrics,
                                                  rnaseq_metrics=rnaseq_metrics)
          graphs = []
          if insert_graph and file_exists(insert_graph):
              graphs.append((insert_graph,
                             "Distribution of paired end insert sizes"))
          return summary_info, graphs
docs: Produce report metrics for a RNASeq experiment using Picard
          with a sorted aligned BAM file.