old_contents: from sequana.pacbio import BAMPacbio
  from sequana import sequana_data
  from easydev import TempFile
  
  def test_pacbio():
      b = BAMPacbio(sequana_data("test_pacbio_subreads.bam"))
      assert len(b) == 8
      b.df
      assert b.nb_pass[1] == 8
  
      with TempFile() as fh:
          b.stride(fh.name, stride=2)
  
      with TempFile() as fh:
          b.filter_length(fh.name, threshold_min=500)
  
      b.hist_snr()
      b.hist_ZMW_subreads()
  
      b.hist_len()
      b.hist_GC()
      b.plot_GC_read_len()
      
  
new_contents: from sequana.pacbio import BAMPacbio
  from sequana import sequana_data
  from easydev import TempFile
  
  def test_pacbio():
      b = BAMPacbio(sequana_data("test_pacbio_subreads.bam"))
      assert len(b) == 8
      b.df
      assert b.nb_pass[1] == 8
  
  
      with TempFile() as fh:
          b.filter_length(fh.name, threshold_min=500)
  
      print(b)   #  check length
  
      b.stats
  
  
      # test nb_pass from scratch
      b = BAMPacbio(sequana_data("test_pacbio_subreads.bam"))
      b.nb_pass
  
      # test hist_snr from scratch
      b._df = None
      b.hist_snr()
  
      # test hist_len from scratch
      b._df = None
      b.hist_len()
  
      # test from scratch
      b._df = None
      b.hist_GC()
  
      # test from scratch
      b._df = None
      b.plot_GC_read_len()
  
      # test from scratch
      b._df = None
      b._nb_pass = None
      b.hist_ZMW_subreads()
  
  def test_pacbio_stride():
      b = BAMPacbio(sequana_data("test_pacbio_subreads.bam"))
      with TempFile() as fh:
          b.stride(fh.name, stride=2)
      with TempFile() as fh:
          b.stride(fh.name, stride=2, random=True)
  
subject: update and add tests