bio-chipseq-visualization
Visualize ChIP-seq data using deepTools, Gviz, and ChIPseeker. Create heatmaps, profile plots, and genome browser tracks. Visualize signal a…
Maintainer FreedomIntelligence · Last updated April 1, 2026
CRISPR screen: guide counting → MAGeCK → hit calling → visualization.
Original source
https://github.com/FreedomIntelligence/OpenClaw-Medical-Skills/tree/main/skills/bio-workflows-crispr-screen-pipeline
Skill Snapshot
Source Doc
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
counts = pd.read_csv('experiment.count.txt', sep='\t', index_col=0)
counts_numeric = counts.iloc[:, 1:]
qc_stats = {}
for col in counts_numeric.columns:
total = counts_numeric[col].sum()
zeros = (counts_numeric[col] == 0).sum()
gini = calculate_gini(counts_numeric[col].values)
qc_stats[col] = {'total_reads': total, 'zero_count_guides': zeros, 'gini': gini}
qc_df = pd.DataFrame(qc_stats).T
print('QC Summary:')
print(qc_df)
## QC thresholds
assert qc_df['zero_count_guides'].max() < len(counts) * 0.2, 'Too many zero-count guides'
assert qc_df['gini'].max() < 0.4, 'Gini index too high (uneven distribution)'
print('QC passed!')
mageck test
-k experiment.count.txt
-t Day14_Rep1,Day14_Rep2,Day14_Rep3
-c Day0
-n negative_screen
--pdf-report
--gene-lfc-method alphamedian
Related skills
Visualize ChIP-seq data using deepTools, Gviz, and ChIPseeker. Create heatmaps, profile plots, and genome browser tracks. Visualize signal a…
Generate consensus FASTA sequences by applying VCF variants to a reference using bcftools consensus. Use when creating sample-specific refer…
Visualize copy number profiles, segments, and compare across samples. Create publication-quality plots of CNV data from CNVkit, GATK, or oth…
Circular genome visualization with Circos or pycirclize.