I have RNAseq data taken from two different ENA project IDs (PrjA & PrjB). Both PrjA & PrjB contains test and control samples. PrjA was PE 50bp sequencing while PrjB was SE 100bp. Is it possible we can combine the data to identify DEG between test and control? Would batch effect removal by Combat be a better option or use the design = project+condition in DESeq2?