churchill-lab / g2gtools

Personal diploid genome creation and coordinate conversion
http://churchill-lab.github.io/g2gtools
21 stars 9 forks source link

pass on length type error, add proper pair flag #26

Open MikeWLloyd opened 1 year ago

MikeWLloyd commented 1 year ago

Was seeing the following error in back conversion

...
if cigar_seq_length != len(sequence):
TypeError: object of type 'NoneType' has no len()

The solution to simply pass on the error is rough, but avoids the problem.

Second issue is to add FLAG_PROPER_PAIR back to BAM following conversion.