deeptools 3.2.1 plotProfile / multiqc

[nservant]

Hi,

I think there is a small bug in the plotProfile output (v3.2.1) and MultiQC.

  computeMatrix reference-point \\                                                                                                                                                                                 
                --referencePoint TSS \\                                                                                                                                                                            
                -R $geneBed \\                                                                                                                                                                                     
                -S ${bigwig} \\                                                                                                                                                                                    
                -o ${prefix}_matrix.mat.gz \\                                                                                                                                                                      
                --outFileNameMatrix ${prefix}.computeMatrix.vals.mat.gz \\                                                                                                                                         
                -b 1000 -a 1000 --skipZeros \\                                                                                                                                                                     
                -p ${task.cpus}                                                                                                                                                                                    
                                                                                                                                                                                                                   
  plotProfile -m ${prefix}_matrix.mat.gz \\                                                                                                                                                                        
              -o ${prefix}_bams_profile.pdf \\                                                                                                                                                                     
              --outFileNameData ${prefix}.plotProfile.tab                                                                                                                                                          

MultiQC crashes because it is looking for the 'TSS' index, which is not in the output file.

In the output, I have

bin labels		-1.0Kb	tick 	1.0Kb

Replacing 'tick' by 'TSS' fixes the issue Best Nicolas

[dpryan79]

Hmm, yeah, I'll have to think about how to change multiQC to handle that. I think I added a header line in recent deepTools versions, so maybe that's what can be used for file detection in multiQC.