1- /* clean-hairpins: a program for identifying and removing hairping
2- * reads from paired-end WGBS or RRBS reads.
1+ /* clean-hairpins: a program for identifying and removing hairpin
2+ * reads from paired-end WGBS/PBAT/RPBAT or RRBS reads.
33 *
4- * Copyright (C) 2018-2022 Andrew D. Smith, Liz Ji and Jenny Qu
4+ * Copyright (C) 2024 Andrew D. Smith
55 *
6- * Authors: Liz Ji and Andrew D. Smith
6+ * Author: Andrew D. Smith
77 *
88 * This program is free software: you can redistribute it and/or
99 * modify it under the terms of the GNU General Public License as
@@ -121,7 +121,7 @@ similar_letters_bisulfite_tc_and_ag(const char a, const char b) {
121121
122122// compare two reads to detect the overlapped region
123123static double
124- similarity_both_bisulfite_convsersions (const string &s1, const string &s2) {
124+ similarity_both_bisulfite_conversions (const string &s1, const string &s2) {
125125 const size_t lim = min (size (s1), size (s2));
126126
127127 uint32_t total_letters = 0 ;
@@ -263,7 +263,7 @@ clean_hairpin::analyze_reads(const string &outfile1, const string &outfile2,
263263
264264 // See if inverted duplicates emerge
265265 const double percent_match =
266- similarity_both_bisulfite_convsersions (r1.seq , r2.seq );
266+ similarity_both_bisulfite_conversions (r1.seq , r2.seq );
267267
268268 // ADS: need a bitter way to get this bin identifier
269269 ++hist[floor (percent_match * n_hist_bins)];
@@ -302,7 +302,7 @@ clean_hairpin::analyze_reads(bgzf_file &in1, bgzf_file &in2,
302302 ++hps.n_good_reads ;
303303
304304 const double percent_match =
305- similarity_both_bisulfite_convsersions (r1.seq , r2.seq );
305+ similarity_both_bisulfite_conversions (r1.seq , r2.seq );
306306 ++hist[floor (percent_match * n_hist_bins)];
307307
308308 if (percent_match > cutoff) {
@@ -318,7 +318,7 @@ clean_hairpin::analyze_reads(bgzf_file &in1, bgzf_file &in2,
318318int
319319main_clean_hairpins (int argc, const char **argv) {
320320
321- static const string description = " fix and stat invdup/hairping reads" ;
321+ static const string description = " fix and stat invdup/hairpin reads" ;
322322
323323 try {
324324
0 commit comments