I\'ve had this issue with the naming (contrast) of the deseq2 output, where the naming out the output was not reproducible. Normally the comparison would look like: variable_lev