blachlylab / fade

Fragmentase Artifact Detection and Elimination
MIT License
11 stars 3 forks source link

negative values for 'read count', 'Clipped %', etc #29

Open rhalperin opened 2 years ago

rhalperin commented 2 years ago

I am seeing a few bams were all of the values printed to STDERR from fade out are negative. For example:

read count: -1875829843
Clipped %:  -0.297902
% With Supplementary alns:  -0.00544583
Artifact rate:  -0.00885916
% With Supplementary alns and artifacts:    -0.00333769
Artifact rate left only:    -0.00442794
Artifact rate right only:   -0.00443122

Unfortunately, this issue doesn't replicate on the same bam subset to a region, so I'm not able to share data to reproduce the issue.

charlesgregory commented 2 years ago

I check to see if there is potential for those stats to overflow or somehow become negative.