search

PacificBiosciences / pb-CpG-tools

Collection of tools for the analysis of CpG data
BSD 3-Clause Clear License
62 stars 5 forks source link

strand specific extraction #37

Closed xiaoyunguo closed 1 year ago

xiaoyunguo commented 1 year ago

Hi,

I want to do strand specific methylation extraction, I wonder if there is a way to do it? say setting the haplotype tag in some way?

Thanks

Xiaoyun

ctsa commented 1 year ago

Hi Xiaoyun, Do you mean the read's alignment orientation? Although it is not supported, you are right that you could have the output arbitrarily partitioned using the haplotype tag, just add HP aux tag values of 1 and 2 two the two sets of reads for which you would like separate tracks.

ctsa commented 1 year ago

Closing as answered, please reopen if needed.