Sample swap in modkit dmr compare-tx-sites

[Kilian-Winiger]

Hi!

The modkit dmr compare-tx-sites seems to swap samples. It reports the counts of sample a for sample b and vice versa. This reverses also effect size and log2 fold change.

May it be that these lines cause that?

        let cond_a_counts = res.isoform_counts.pop().unwrap();
        let cond_b_counts = res.isoform_counts.pop().unwrap();

GeneDmrScore

[ArtRand]

Hello @Kilian-Winiger,

Thanks again for the catch. Here is a build with the change as you suggested. I've branched this from your previous issue, so you should be able to use this build instead of the one attached here. You can find the code to this build on this branch.

Thanks.

modkit_v0.6.3_2273723.zip