Improve tag handling.

mapdamage/rescale.py

@@ -275,11 +275,12 @@ def rescale_qual_read(bam, read, ref, corr_prob,subs, debug = False,direction="b
     read.qual = new_qual
     # truncate this to 5 digits
     number_of_rescaled_bases = float("%.5f" % number_of_rescaled_bases)
-    # check if the read has a MR tag 
-    for tag in read.tags:
-        if tag[0] == "MR":
-            raise SystemExit("Read: %s already has a MR tag, can't rescale" % (read))
-    read.tags = read.tags + [("MR:f",number_of_rescaled_bases)]
+
+    if read.has_tag("MR"):
+        raise SystemExit("Read: %s already has a MR tag, can't rescale" % read)
+
+    read.set_tag("MR", number_of_rescaled_bases, 'f')
+
     return read
 
 

mapdamage/version.py

@@ -2,5 +2,4 @@
 try:
     from mapdamage._version import __version__
 except ImportError:
-    __version__ = "2.0.7"
-
+    __version__ = "2.0.8"

setup.py

@@ -55,7 +55,7 @@ def run(self):
 setup(
     cmdclass={'install': compileInstall},
     name='mapdamage',
-    version='2.0.7',
+    version='2.0.8',
     author='Aurélien Ginolhac, Mikkel Schubert, Hákon Jónsson',
     author_email='MSchubert@snm.ku.dk, jonsson.hakon@gmail.com',
     packages=['mapdamage'],