search

getzlab / rnaseqc

Fast, efficient RNA-Seq metrics for quality control and process optimization
Other
146 stars 19 forks source link

Faulty estimates of Read Length #77

Open karl-az opened 1 year ago

karl-az commented 1 year ago

I have noticed that RNA-seQC at times underestimates the Read Length

This is the section measuring read length:

https://github.com/getzlab/rnaseqc/blob/c426536132b7f90e108707c5e4da206097e31ad7/src/RNASeQC.cpp#L274-L277

The comparison is on alignmentSize, but the assignment is with alignment.Length(). Could this lead to unfortunate errors? Especially if alignment.Length() sometimes is less than alignmentSize

kachulis commented 1 year ago

I believe this is addressed by #78