Boxplot

Genome Sequencing Accuracy

Read accuracy distribution across sequencing platforms

Output
Genome Sequencing Accuracy
Python
import matplotlib.pyplot as plt
import numpy as np

np.random.seed(42)
platforms = ['Illumina', 'PacBio', 'Nanopore', 'Ion Torrent', 'BGI']
data = [
    99.9 - np.random.exponential(0.05, 200),
    99.5 - np.random.exponential(0.3, 200),
    98 - np.random.exponential(1, 200),
    99.2 - np.random.exponential(0.2, 200),
    99.7 - np.random.exponential(0.1, 200)
]

fig, ax = plt.subplots(figsize=(10, 6), dpi=100)
ax.set_facecolor('#0d1117')
fig.patch.set_facecolor('#0d1117')

colors = ['#27D3F5', '#F5276C', '#6CF527', '#F5B027', '#4927F5']
bp = ax.boxplot(data, widths=0.5, patch_artist=True, showfliers=True, notch=True,
                flierprops=dict(marker='v', markerfacecolor='#9C2007', markersize=4, alpha=0.6),
                medianprops=dict(color='#F5D327', linewidth=2))

for patch, color in zip(bp['boxes'], colors):
    patch.set_facecolor(color)
    patch.set_alpha(0.8)
    patch.set_edgecolor('white')
for i, color in enumerate(colors):
    bp['whiskers'][i*2].set_color(color)
    bp['whiskers'][i*2+1].set_color(color)
    bp['caps'][i*2].set_color(color)
    bp['caps'][i*2+1].set_color(color)

ax.axhline(99, color='#F54927', linewidth=1, linestyle='--', alpha=0.7, label='99% Threshold')
ax.set_xticklabels(platforms, rotation=15, ha='right')
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.spines['left'].set_color('#333333')
ax.spines['bottom'].set_color('#333333')
ax.yaxis.grid(True, color='#1a1a3f', linewidth=0.5, zorder=0)
ax.set_axisbelow(True)
ax.tick_params(axis='both', colors='#888888', labelsize=9, length=0, pad=8)
ax.set_ylabel('Read Accuracy (%)', fontsize=11, color='white', fontweight='500')
ax.set_title('Genome Sequencing Platform Accuracy', fontsize=14, color='white', fontweight='bold', pad=15)
ax.legend(loc='lower right', facecolor='#0d1117', edgecolor='#333333', labelcolor='white')

plt.tight_layout()
plt.show()
Library

Matplotlib

Category

Statistical

Did this help you?

Support PyLucid to keep it free & growing

Support