Skip to content

Improve plot-bamstats qual-hm#2143

Merged
whitwham merged 2 commits into
samtools:developfrom
jkbonfield:plot-bamstats-colours
Dec 11, 2024
Merged

Improve plot-bamstats qual-hm#2143
whitwham merged 2 commits into
samtools:developfrom
jkbonfield:plot-bamstats-colours

Conversation

@jkbonfield

Copy link
Copy Markdown
Contributor

This combines @jgrg's #2116 with my own additions to cope with long-read technologies and to fix a bug introduced with single end libraries.

Before and after screenshots for the combined commits.

Novaseq

image
image

Revio

image
image

Ultima

image
image

Fixes warning "warning: y2 axis range undefined or overflow, resetting to [0:0]"
which resulted in missing y-ticks.

Fixed quals-hm layout to be 2 distinct plots.

Adjusts the colour of the heatmap to be easier to read, avoiding
multiple dark zones in the colour wheel.

Fixes broken links in the HTML table for characters that need URI escaping.

Reduce width of plot so that "Number of bases" heatmap label is not clipped
Comment thread misc/plot-bamstats Outdated
$args = get_defaults($opts,"$$opts{prefix}quals2.png",wh=>$is_paired ? '700,500' : '600,400');
$fh = $$args{fh};
my $pos_size = $is_paired ? " set rmargin 0; set lmargin 0; set tmargin 0; set bmargin 0; set origin 0.1,0.1; set size 0.4,0.8" : '';
# my $pos_size = $is_paired ? " set rmargin 0; set lmargin 0; set tmargin 0; set bmargin 0; set origin 0.1,0.1; set size 0.4,0.8" : '';

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Does this line still need to be in here?

Fix single-ended run plots.

Reduce rate of tick labels for long-reads.

Adjust colour wheel to bring out more fidelity with very low qualities
so we have better distinction between all zero and some low qual.

Replaces / fixes samtools#2116
@jkbonfield jkbonfield force-pushed the plot-bamstats-colours branch from 1d89c57 to f0964c3 Compare December 11, 2024 09:50
@whitwham whitwham merged commit aa8cb59 into samtools:develop Dec 11, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants