Skip to content
Snippets Groups Projects
Commit f05d7cb4 authored by Ruben Vorderman's avatar Ruben Vorderman
Browse files

Use set_tag call from pysam

parent 054b7a7f
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1 addition and 2 deletions
umi.wdl
+ 1
− 2
View file @ f05d7cb4
...@@ -62,8 +62,7 @@ task BamReadNameToUmiTag { ...@@ -62,8 +62,7 @@ task BamReadNameToUmiTag {
for segment in in_bam: # type: pysam.AlignedSegment for segment in in_bam: # type: pysam.AlignedSegment
new_name, umi = split_umi_from_name(segment.query_name) new_name, umi = split_umi_from_name(segment.query_name)
segment.query_name = new_name segment.query_name = new_name
# append does not work. (Pysam is not Pythonic.) segment.set_tag("RX", umi, value_type="Z")
segment.tags = segment.tags + [(bam_tag, umi)]
out_bam.write(segment) out_bam.write(segment)
if __name__ == "__main__": if __name__ == "__main__":
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment