diff --git a/CLAM/preprocessor.py b/CLAM/preprocessor.py index 503c2c7..1511ad9 100755 --- a/CLAM/preprocessor.py +++ b/CLAM/preprocessor.py @@ -116,6 +116,7 @@ def filter_bam_multihits(filename, max_tags, max_hits, out_dir, read_tagger_meth if read_tag==-1: continue read.tags += [('RT', read_tag)] ## add the tag + tagged_read = pysam.AlignedSegment() tagged_read.query_name = read.query_name tagged_read.query_sequence = 'N' @@ -239,8 +240,8 @@ def filter_bam_maxtags(obam_fn, ibam_fn, max_tags=1): last_pos = read.positions[0] stack = [read] for new_alignment in new_alignment_list: - new_alignment.query_sequence = '*' - new_alignment.query_qualities = '0' + #new_alignment.query_sequence = '*' + #new_alignment.query_qualities = '0' _ = obam.write(new_alignment) else: stack.append(read) @@ -248,8 +249,8 @@ def filter_bam_maxtags(obam_fn, ibam_fn, max_tags=1): output_counter += len(new_alignment_list) last_pos = read.positions[0] for new_alignment in new_alignment_list: - new_alignment.query_sequence = '*' - new_alignment.query_qualities = '0' + #new_alignment.query_sequence = '*' + #new_alignment.query_qualities = '0' _ = obam.write(new_alignment) ibam.close() obam.close()