1、去除冗余的reads

for i in control_P_R1 Insulin_hcg_P_R1  insulin_P_R1;do deduplicate_bismark --bam $i.clean_bismark_bt2_pe.bam --output_dir  ~/WGBS/dedup

結果：

生成deduplicated.bam和deduplication_report.txt文件

2、將bam文件轉化為sam文件

for i in control_P_R1 Insulin_hcg_P_R1 insulin_P_R1;do samtools view – h $i.clean_bismark_bt2_pe.deduplicated.bam > $i.clean_bismark_bt2_pe.deduplicated.sam;done

結果：

結果

3、計算每個胞嘧啶的甲基化值

bismark_methylation_extractor --bedGraph -p --comprehensive --no_overlap --cytosine_report --counts --report --genome_folder ~/WGBS/ref control_P_R1.clean_bismark_bt2_pe.deduplicated.sam -o ~/WGBS/extr1

結果：

結果