from sequana import Contigs, sequana_data
filename = sequana_data("test_contigs_ex1.fasta")
c = Contigs(filename)
c.plot_scatter_contig_length_vs_nreads_cov()