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
結(jié)果：
生成deduplicated.bam和deduplication_report.txt文件
2夺艰、將bam文件轉(zhuǎn)化為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
結(jié)果：
結(jié)果
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
結(jié)果：
結(jié)果
