Error in rule AberrantSplicing_Overview_R when padjCutoff too low:

[noeldjitro]

Hi,

when padjCutoff is too low, the results_per_junction.tsv and results.tsv is header only.
This caused an error in my analysis because "type" is not found in results_per_junction.tsv.

I managed to replicate a similar error by changing the "padjCutoff" from 1 to 0.1 in the drop demo config.yaml:

• The header in results.tsv:

• However, the header in results_per_junction.tsv is only:

• this is causing error in AberrantSplicing_Overview_R because type is not found in results_per_junction.tsv header

The results_per_junction.tsv and results.tsv are generated from previous script: 08_extract_results_FraseR.R

Can you fix this problem, please?

Thank you,
Noel

[vyepez88]

Hi Noel, we'll fix it on our side.
You can overcome it by simply running the module independently: snakemake aberrantSplicing