1. 将read中的OT,OB的所有C化成T,然后将CTOT,CTOB中所有G转化为A
2. 将reference的T,B所有C转化成T, 然后再将CTT,CTB的所有G转化为A, 生成一个bisulfide基因组
问题
为什么不能不用CT的reads, 这个毫无意义啊!反正和O是绝对的互补的
可能nondirectional library会有更多的信息
3.将reads和bisulfide基因组比对,找到其在bisulfide基因组上最适位点
可能的比对的结果
比对可能产生unique best match红色
4.再将原来的reads和该位置的正常基因组比对,找到甲基化的位点