SequenceWithQualities does not handle qualstr=None

2013-06-12T06:35:38.316000Z

Reading a BAM with a "*" in the quality causes SequenceWithQualities to throw an exception: object of type 'NoneType' has no len() in init() on the line:

if len( seq ) != len( qualstr ):

According to pysam:

http://www.cgat.org/~andreas/documentation/pysam/api.html#pysam.AlignedRead

The "qual" field can be None if not present (ie a "*" in the SAM)

Thus the fix should be (I think) from:

if read.qual != "*":
         seq = SequenceWithQualities( read.seq, read.qname, read.qual )

to:

if read.qual is not None and read.qual != "*":
         seq = SequenceWithQualities( read.seq, read.qname, read.qual )

Recent changes to 12: SequenceWithQualities does not handle qualstr=None

SequenceWithQualities does not handle qualstr=None