@@ -156,12 +156,11 @@ eats_query(const std::uint32_t c) {
156156 return bam_cigar_type (bam_cigar_op (c)) & 1 ;
157157}
158158
159- /* The three functions below here should probably be moved into
160- bsutils.hpp. I am not sure if the DDG function is needed, but it
161- seems like if one considers strand, and the CHH is not symmetric,
162- then one needs this. Also, Qiang should be consulted on this
163- because he spent much time thinking about it in the context of
164- plants. */
159+ /* The three functions below here should probably be moved into bsutils.hpp. I
160+ am not sure if the DDG function is needed, but it seems like if one
161+ considers strand, and the CHH is not symmetric, then one needs this. Also,
162+ Qiang should be consulted on this because he spent much time thinking about
163+ it in the context of plants. */
165164[[nodiscard]] static inline bool
166165is_chh (const std::string &s, const std::size_t i) {
167166 return i + 2 < std::size (s) && is_cytosine (s[i]) && !is_guanine (s[i + 1 ]) &&
@@ -180,11 +179,10 @@ is_c_at_g(const std::string &s, const std::size_t i) {
180179 !is_guanine (s[i + 1 ]) && is_guanine (s[i + 2 ]);
181180}
182181
183- /* Right now the CountSet objects below are much larger than they need
184- to be, for the things we are computing. However, it's not clear
185- that the minimum information would really put the memory
186- requirement of the program into a more reasonable range, so keeping
187- all the information seems reasonable. */
182+ /* Right now the CountSet objects below are much larger than they need to be,
183+ for the things we are computing. However, it's not clear that the minimum
184+ information would really put the memory requirement of the program into a
185+ more reasonable range, so keeping all the information seems reasonable. */
188186struct CountSet {
189187 static constexpr auto max_prob_repr = 255.0 ;
190188 // ADS: accepting int16_t because of using -1 for unknown prob vs. 0 prob.
0 commit comments