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

Why skip 40 aligned_pairs? #29

Closed bing2022ca closed 1 year ago

bing2022ca commented 1 year ago

In aligned_bam_to_cpg_scores.py line 449: for query_pos, ref_pos in read.get_aligned_pairs(matches_only=True)[20:-20]:

Why skip the first and last 20 pairs here?

ctsa commented 1 year ago

Thanks @bing2022ca. Interesting question but I don't have the original design criteria.