| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112 |
- #######################
- # globe
- #######################
- PEAKSNAKE_HOME: /home/zhxd/workspace/20231225_ChIP_ATAC/PeakSnake
- genome: ref/genome.fa
- gtf: ref/genes_longest.gtf
- tss_tes_bed: ref/genes_longest.bed6
- tss_tes_shuffle_bed: ref/genes_shuffle.bed6
- seq_type: SE # 或者 "PE"
- peak_type: narrow # narrow | broad
- # 有效基因组大小, 参考 https://deeptools.readthedocs.io/en/develop/content/feature/effectiveGenomeSize.html#effective-genome-size
- gsize: 119481543
- # idr: 只支持 2 重复,适合 narrowPeak) ; intersect (2 + 重复 ) ; norep: 没有重复
- peak_selection: idr
- #######################
- # samples
- #######################
- samples:
- ZAT6_ABA_rep1: Col0_mock_rep1
- ZAT6_ABA_rep2: Col0_mock_rep2
- H3K4me3_WT_rep1: Input_WT_rep1
- H3K4me3_WT_rep2: Input_WT_rep1
- contrasts:
- - ZAT6_ABA_vs_H3K4me3_WT
- #######################
- # fastp
- #######################
- fastp: "--adapter_sequence AGATCGGAAGAGCACACGTCTGAACTCCAGTCA --adapter_sequence_r2 AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT --qualified_quality_phred 20 --length_required 50"
- #######################
- # bowtie2
- # PE 建议设置 --no-mixed --no-discordant, 只保留成对比对
- # ATAC-seq 必须设置 --maxins 1000 ~ 2000
- #######################
- bowtie2_index: ref/genome # 不设置则从genome.fa 构建
- bowtie2:
- #######################
- # alignmentSieve
- # blackregion: 黑名单区域 线粒体 叶绿体 大肠杆菌等
- # ATAC-seq 在 extra 中添加 --ATACshift
- #######################
- alignmentSieve:
- minMappingQuality: 25
- blacklist: ref/blacklist.bed
- extra:
- #######################
- # bam to bw
- #######################
- bamCoverage: "--binSize 50 --normalizeUsing RPKM"
- #######################
- # MACS3
- # fold: 用于过滤 peaks
- # qscore: 用于过滤 peaks
- # macs3 callpeak ChIP-seq 一般不需要设置特殊参数
- # ATAC-seq 推荐: --nomodel --shift -100 --extsize 200 或 --nomodel --shift -75 --extsize 150
- #######################
- macs3:
- qscore: 3
- extra:
- #######################
- # replicate intersect
- #######################
- intersect:
- min_overlap: 0.5
- #######################
- # IDR
- #######################
- idr:
- idr_threshold: 0.05
- #######################
- # 提取序列用于 motif 分析
- #######################
- motif:
- top_n: 1000
- summit_flank: 50
- #######################
- # Peank 注释
- # 只关注 TSS 附近 用 ChIPseeker
- # 否则UROPA
- #######################
- uropa:
- do: yes
- feature: ['gene']
- feature_anchor: ['start'] # ['start', 'center', 'end']
- relative_location: ['Upstream', 'OverlapStart'] # ['Upstream', 'Downstream', 'OverlapStart', 'OverlapEnd','PeakInsideFeature', 'FeatureInsidePeak']
- distance: [5000, 5000]
- #######################
- # diff analysis
- #######################
- diff_peaks:
- do: yes
- method: edgeR # edgeR | DESeq2
- padj: 0.05
- #######################
- # cross correlation
- # ChIP-seq 需要进行, ATAC-seq 不需要
- #######################
- cross_correlation:
- do: yes
- spp: "-s=0:5:500"
|
