old_contents: from sequana.tools import bam_to_mapped_unmapped_fastq, reverse_complement, StatsBAM2Mapped
  from sequana import sequana_data
  from sequana.tools import bam_get_paired_distance, GZLineCounter, PairedFastQ
  
  
  def test_StatsBAM2Mapped():
      data = sequana_data("test.bam", "testing")
      res = StatsBAM2Mapped(data)
  
  
  def test_bam2fastq():
      data = sequana_data("test.bam", "testing")
      res = bam_to_mapped_unmapped_fastq(data)
  
  
  
  def test_reverse_complement():
      assert reverse_complement("AACCGGTTA") == 'TAACCGGTT'
  
  
  def test_reverse():
      from sequana.tools import reverse
      assert reverse("AACCGG") == 'GGCCAA'
  
  
  def test_distance():
      data = sequana_data("test.bam", "testing")
      distances = bam_get_paired_distance(data)
  
  
  def test_gc_content():
      from sequana.tools import gc_content
      data = sequana_data('test.fasta', "testing")
      gc_content(data, 10)['seq1']
  
  def test_gzlinecounter():
      assert len(GZLineCounter(sequana_data("test.fastq.gz"))) == 1000
  
  
  def test_paired_file():
      f1 = sequana_data("test.fastq.gz")
      f2 = sequana_data("test.fastq.gz")
  
      assert PairedFastQ(f1,f2).is_synchronised()
  
new_contents: from sequana.tools import bam_to_mapped_unmapped_fastq, reverse_complement, StatsBAM2Mapped
  from sequana import sequana_data
  from sequana.tools import bam_get_paired_distance, GZLineCounter, PairedFastQ
  from sequana.tools import genbank_features_parser
  
  def test_StatsBAM2Mapped():
      data = sequana_data("test.bam", "testing")
      res = StatsBAM2Mapped(data)
      res.to_html()
  
  
  def test_bam2fastq():
      data = sequana_data("test.bam", "testing")
      res = bam_to_mapped_unmapped_fastq(data)
  
  
  def test_reverse_complement():
      assert reverse_complement("AACCGGTTA") == 'TAACCGGTT'
  
  
  def test_reverse():
      from sequana.tools import reverse
      assert reverse("AACCGG") == 'GGCCAA'
  
  
  def test_distance():
      data = sequana_data("test.bam", "testing")
      distances = bam_get_paired_distance(data)
  
  
  def test_gc_content():
      from sequana.tools import gc_content
      data = sequana_data('test.fasta', "testing")
      gc_content(data, 10)['seq1']
      gc_content(data, 101, circular=True)['seq1']
  
  
  def test_genbank_features_parser():
      data = sequana_data("JB409847.gbk")
      genbank_features_parser(data)
  
  
  def test_gzlinecounter():
      assert len(GZLineCounter(sequana_data("test.fastq.gz"))) == 1000
  
  
  def test_paired_file():
      f1 = sequana_data("test.fastq.gz")
      f2 = sequana_data("test.fastq.gz")
      assert PairedFastQ(f1,f2).is_synchronised()
  
subject: add more tests