💊 fix corner case where alignment has no aQual score

c4effa6d · Puriney · c26b7880 · c4effa6d
Commit c4effa6d authored Oct 09, 2018 by Puriney
Hide whitespace changes
Inline Side-by-side

Showing with 1 addition and 1 deletion

celseq2/count_umi.py celseq2/count_umi.py +1 -1

No files found.
--- a/celseq2/count_umi.py
+++ b/celseq2/count_umi.py
@@ -63,7 +63,7 @@ def count_umi(sam_fpath, features, stranded='yes',
        except KeyError:
            pass

-        if aln.aQual < accept_aln_qual_min:
+        if not aln.aQual or aln.aQual < accept_aln_qual_min:
            aln_cnt["_low_map_qual"] += 1
            continue