Rows (gene IDs) in pipeline's Rds expression counts and rows in the recommended external geneCounts.tsv.gz are in different order
read 19 files Error: The rows (genes) of the count matrices to be merged are not the same. Execution halted Error in rule AberrantExpression_pipeline_Counting_mergeCounts_R:
I am using drop 1.5.0 with my own samples (bams obtained from rnaseq nf core pipeline), mapped to gencode29, together with your external samples from https://zenodo.org/records/6078397. I inspected the gene IDs (i.e. rownames) in external as well as in my data and values are exactly the same, just the order of rows is different (and that probably causes failing the pipeline). External rownames (gtex/v8/exported_counts/Whole_Blood--GRCh38--gencode29/geneCounts.tsv.gz) are sorted, but rownames from your RDS files (Output/processed_data/aberrant_expression/v29/counts/somename.Rds) are not.