search

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

find multiple mapping region #35

Open jeongmeani opened 4 years ago

jeongmeani commented 4 years ago

Hi, these days, i analysis bisulfite sequencing using BSseeker2.

but some target reads are aligned to multiple location of genome. so i want to know where is this region.

what can i do?

Best Regards.

Jeongmin

guoweilong commented 4 years ago

Hi Jeongmin,

In BS-Seeker2, it only keeps the mapping results of uniquely mapped reads in final results.

Actually, the mapping results of all mappable reads were in temporary folder, and were removed at the end of running. You can try to edit the code for removing these files as comment, and see all the files left in temporary folder.

Best, Weilong