kineticsTools
kineticsTools copied to clipboard
ipdSummary doesn't report individual values that are used for tMean calculation
I am looking for a way to extract individual values that are used to calculate tMean. For example, if there are 5 reads aligned in a certain region, each value of tMean should be calculated from 5 values. I need all those individual values.
Alternatively, how could I get IPDs with respect to my reads instead of my reference? That way, for each read, I would have vector of IPD values in the length of the read (this of course would have to account for errors, especially small insertions in the reads that do not contribute to the tMean calculation).
I will appreciate your help.