Route: atac

ATAC-seq analysis.

Segments:

• Align to the reference genome (Bowtie2).
• Remove duplicate reads (Sambamba).
• Generate genome browser tracks.
• Call peaks (MACS and HMMRATAC).

Usage

Set up a new analysis (common across all routes). If running for the first time, check the detailed usage instructions for an explanation of every step.

cd <project dir>
git clone --depth 1 https://github.com/igordot/sns
sns/generate-settings <genome>
sns/gather-fastqs <fastq dir>

Run atac route.

sns/run atac

Check for potential problems.

grep "ERROR:" logs-sbatch/*

Output

Primary results:

• BAM-DD: Deduplicated BAM files. Can be used for additional analysis.
• BIGWIG: BigWig files normalized to the total number of reads. Can be used for visual inspection of peaks.
• peaks-*: Peaks.

Run metrics:

• summary-combined.atac.csv: Summary table that includes the number of reads, alignment rate, and fraction of PCR duplicates.
• summary.qc-fragment-sizes.png: Distribution of fragment sizes.