from sequana import sequana_data
from sequana.compare import RNADiffCompare

c = RNADiffCompare(
    sequana_data("rnadiff_salmon.csv", "doc/rnadiff_compare"),
    sequana_data("rnadiff_bowtie.csv", "doc/rnadiff_compare")
)
c.plot_volcano()
indices = c.r1.df.query("log2FoldChange>1 or log2FoldChange<-1").index.values
indices  = [x for x in indices if x in c.r1.df.index and x in c.r2.df.index]
c.plot_geneset(indices, showlines=True)