BSSeeker / BSseeker2

A versatile aligning pipeline for bisulfite sequencing data
http://pellegrini.mcdb.ucla.edu/BS_Seeker2/
MIT License
61 stars 25 forks source link

Use .bam files generated with bwa-mem or re-make them? #49

Open desmodus1984 opened 9 months ago

desmodus1984 commented 9 months ago

Hi,

I am interested in detecting methylation regions and I have been using MethylKit. The method is not very powerful and I want to use more powerful approaches to detect differential methylation. One way is with HOME (https://github.com/ListerLab/HOME). It mentions using input from BSseeker2. I would like to know if I can use the bam files generated with bwa-mem, and read-deduplicated and bam-converted with samtools, or if I should do the entire process from the beginning - read trimming/ alignment.

Thank you very much;

Juan Pablo