Skip to content

Commit

Permalink
updated customemodule for more multiqc output
Browse files Browse the repository at this point in the history
  • Loading branch information
U13bs1125 committed Jul 9, 2024
1 parent 8051466 commit 20b5797
Showing 1 changed file with 18 additions and 0 deletions.
18 changes: 18 additions & 0 deletions modules/local/custommodule.nf
Original file line number Diff line number Diff line change
Expand Up @@ -20,6 +20,24 @@ process CUSTOMMODULE {
script:
def args = task.ext.args ?: ''
"""
echo "# id: 'repeat summary'" > repeatmqc_mqc.tsv
echo "# section_name: 'repeat masking summary statistics'" >> repeatmqc_mqc.tsv
echo "# format: 'tsv'" >> repeatmqc_mqc.tsv
echo "# plot_type: 'bargraph'" >> repeatmqc_mqc.tsv
echo "# description: 'This plot shows a brief summary of each genomes whose repeats has been masked'" >> repeatmqc_mqc.tsv
echo "# pconfig:" >> repeatmqc_mqc.tsv
echo "# id: 'repeat summary'" >> repeatmqc_mqc.tsv
echo "# title: 'repeat summary'" >> repeatmqc_mqc.tsv
echo "# ylab: ''" >> repeatmqc_mqc.tsv
echo "id\tTotal scaffold length\tTotal contig length\ttantan masked bases\twindowmasker masked bases\trmodeler masked bases" >> repeatmqc_mqc.tsv
printf "\$(basename $assemt .assembly_summary)\t" >> repeatmqc_mqc.tsv
grep 'Total scaffold length' $assemt | tail -n 1 | awk '{print \$4}' | tr '\n' '\t' >> repeatmqc_mqc.tsv
grep 'Total contig length' $assemt | tail -n 1 | awk '{print \$4}' | tr '\n' '\t' >> repeatmqc_mqc.tsv
grep 'soft-masked bases' $assemt | tail -n 1 | awk '{print \$4}' | tr '\n' '\t' >> repeatmqc_mqc.tsv
grep 'soft-masked bases' $assemw | tail -n 1 | awk '{print \$4}' | tr '\n' '\t' >> repeatmqc_mqc.tsv
grep 'soft-masked bases' $assemr | tail -n 1 | awk '{print \$4}' >> repeatmqc_mqc.tsv
echo "# id: 'tantan repeat summary'" > tantanmqc_mqc.tsv
echo "# section_name: 'tantan repeat masking summary statistics'" >> tantanmqc_mqc.tsv
echo "# format: 'tsv'" >> tantanmqc_mqc.tsv
Expand Down

0 comments on commit 20b5797

Please sign in to comment.