Multiqc: Trimmomatic logs not parsed properly

Created on 26 Nov 2018  路  4Comments  路  Source: ewels/MultiQC

Description of bug:
I try to assemble a report with Trimmomatic logs (amongst others). Trimmomatic was run in Galaxy, all the logs therefore refer to a fastq_r1 input file used to name the sample. To circumvent this, I am using the instruction:

extra_fn_clean_exts:
    - type: regex_keep
      pattern: '[0-9]+_mATAC_t[0-9]{1,2}'

to extract sample names from file names like 180518_mATAC_t1_1_trimmomatic.err.

Strangely this instruction is ignored for trimmomatic logs while the regex perfectly works for other modules ; so I end up with a single sample (i.e. fastq_r1) in the report section.

MultiQC Error log:

[INFO   ]         multiqc : This is MultiQC v1.6
[INFO   ]         multiqc : Template    : default
[INFO   ]         multiqc : Searching 'qc/trimmomatic/'
[INFO   ]     trimmomatic : Found 1 logs
[INFO   ]         multiqc : Compressing plot data
[INFO   ]         multiqc : Report      : trimmo_test.html
[INFO   ]         multiqc : Data        : trimmo_test_data
[INFO   ]         multiqc : MultiQC complete

File that triggers the error:
Attached is a zip with 10 logs and the config file

trimmomatic.zip

MultiQC run details (please complete the following):

  • Command used to run MultiQC:
    multiqc -c multiqc_ATAC_config.yaml -n trimmo_test qc/trimmomatic/
  • MultiQC Version: MultiQC v1.6
  • Operating System: CentOS Linux release 7.4.1708 (Core)
  • Python Version: Python 2.7.15
  • Method of MultiQC installation: conda
bug

All 4 comments

Thanks for the detailed issue 馃憤 I'll look into it as soon as I can.

Phil

Hi @cgirardot,

The problem here is that Trimmomatic gets it's sample names from the command line used. As you pointed out, this is the same for every file in your case. The extra_fn_clean_exts only works with the sample name that has been collected, so doesn't help in this case as the filename is not used at all.

To get around this, I've added a new option to the code which allows you to set a config flag to tell the module to use the filename instead of the file contents. To use it, add the following config option:

trimmomatic:
    s_name_filenames: true

Once you have this, everything should work as excepted.

Let me know how you get on!

Phil

right on time fro Christmas !

Thx

this is working great , thanks again.

Was this page helpful?
0 / 5 - 0 ratings

Related issues

berguner picture berguner  路  5Comments

cyklee picture cyklee  路  6Comments

pellonen picture pellonen  路  8Comments

MattBashton picture MattBashton  路  7Comments

tluquez picture tluquez  路  6Comments