Count tail lengths — count_tails • tailquant

Count tail lengths observed at each site, and how many are considered censored.

Usage

count_tails(sited_reads, min_tail, length_trim, must_be_close_to_site)

min_tail: Tail lengths shorter than this are discarded.
length_trim: This many bases are effectively trimmed from the ends of reads.
must_be_close_to_site: Are tail lengths only valid if reads ended very close to the site position. This applies to tail lengths from read 1, but not tail lengths from read 2.

If umi column is present, weight each read by 1/n within UMI groups.

UMI counts are n_event, n_died.

Read counts are n_read_event, n_read_died.